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FOREWORD 


During  the  meeting  of  the  AGARD  Fluid  Dynamics  Panel  in  Rome,  Italy,  September  1974,  one  day  was  de- 
voted to  the  topic  of  Multiple  Flow  of  Solids  and  Gases  to  review  the  activities  in  this  field  in  various 
countries.  Five  survey  papers  were  presented  by  representatives  of  France,  West  Germany,  United  Kingdom, 
Belgium,  and  the  United  States.  Subsequently,  the  undersigned  was  Invited  to  prepare  an  AGARDogrsph  based 
on  these  presentations. 

The  papers  were  prepared  by  the  individual  authors  and  are  collected  here  in  the  order  in  which  they 
were  presented  at  the  meeting.  Minor  updating  and  clarifications  were  incorporated  in  cooperation  with 
these  authors.  The  final  text,  was  established  with  the  concurrence  of  the  Editor,  Ing.  Gen.  A.  Auriol, 
Directeur  de  l'Institut  Franco-Allemand  de  Recherches  de  Saint-Louis,  France. 

Several  topics  are  discussed  in  more  than  one  paper.  Therefore,  cross  references  were  inserted  in  the 
text  to  facilitate  comparison  of  different  points  of  view.  These  identify  the  papers  consecutively  by  A to 
E and  the  appropriate  subsections.  In  addition,  a siiject  index  was  prepared  to  allow  quick  location  of 
specific  topics. 

It  's  a great  pleasure  to  thank  the  Fluid  Dynamics  Panel  of  AGARD  for  the  invitation  to  prepare  this 
material  and  to  acknowledge  the  cooperation  of  the  authors  of  the  papers,  of  Ing.  Gen.  A.  Auriol,  and  of 
Mr.  J.  A.  Lawford  and  Mr.  M.  C.  Fischer,  Executives  of  the  Fluids  Dynamics  Panel.  The  help  of  Textron 
Bell  Aerospace  Company  and  the  State  University  of  New  York  at  Buffalo,  New  York,  in  the  preparation  of 
the  manuscript  is  greatly  appreciated. 
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I 


George  Rudlnger 


Buffalo,  New  York 
September  1976 
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A.  FRENCH  CONTRIBUTION  TO  AERODYNAMICS  OF  GAS-PARTICLE  MIXTURES 


Paul  Kuentzmann 

Chef  de  Groupe  de  Recherche  A l'O.N.E.R.A. 

Office  National  d'Etudes  et  de  Recherches  Aerospa tiales  (ONERA) 
92320  Chafillon,  France 


SUMMARY 

Mixed  gas-particle  flows  were  studied  in  France  mainly  for  particular  applications.  Three  examples, 
concerning  the  aerospace  field,  are  presented;  they  concern  rocket  propulsion,  laser  anemometry  and  capture 
of  droplets.  From  the  main  problems  encountered  in  these  studies,  it  has  been  possible  to  shed  some  light 
on  the  most  important  fundamental  problems. 

A precise  knowledge  of  the  particle  size  distribution  is  essential  in  most  cases,  and  improvements  in 
optical  techniques  are  desirable.  Velocity  measurements , satisfactory  for  small  concentrations,  should  be 
extended  to  larger  ones.  Temperature  measurement  methods  for  both  phases  exist.  Interactions  between 
particles  should  be  better  known. 

TVio-phase  flows,  the  importance  of  which  is  expected  to  increase  in  the  future,  should  be  at  first  the 
subject  of  two  types  of  research:  fundamental  studies  of  elementary  mechanisms  and  development  of  new  ex- 

perimental means. 


1.  INTRODUCTION 


1.1  Ceneral 


Gas-particle  mixtures  belong  to  the  category  of  dispersions,  that  is  of  media  made  of  a cont'njui 
fluid  phase  and  of  another  phase  distributed  among  discrete  volume  elements.  There  exist,  a prio.  , five 
possible  combinations  according  to  the  nature  of  the  two  phases;  they  can  be  classified  by  consider »ng 
first  the  fluid  phase: 

- Liquid  fluid  phase: 

• Solid  dispersed  phase:  suspension 

• Liquid  dispersed  phase:  emulsion 

. Gaseous  dispersed  phase:  foam 

- Gaseous  fluid  phase: 

• Solid  dispersed  phase:  suspension,  smoke,  fluidized  bed 

• Liquid  dispersed  phase:  mist,  fog 

We  shall  concentrate  here  on  aerodynamic  phenomena,  i.e.,  on  dispersions  in  which  the  fluid  phase  is  gas- 
eous. Mixtures  in  which  the.  fluid  phase  is  liquid,  the  dispersed  phase  being  liquid  or  solid,  can  be  treated 
by  the  same  methods  as  the  gas-particle  mixture;  but  the  properties  of  liquid-gas  mixtures  can  be  more  varied. 

A detailed  listing  of  gas-particle  mixtures  can  be  established  according  to  the  two  following  para- 
meters: the  volume  fraction  occupied  by  the  condensed  phase,  and  the  mean  dimensions  of  the  dispersed 

elements.  This  classification  is  given  in  the  following  table. 

Classification  of  the  mixtures  of  gas  and  particles  (1). 


Particle 

State 

Name  of  the 
Mixture 

Partial  Volume 
Occupied  by  the 
Condensed  Phase 

Mean  Size  of 
the  Particles 

Liquid 

Fog 

« i 

< O.l^m 

Cloud 

« i 

> 0.  U*m 

Aerosol 

« i 

\ very  small 

Solid 

Aerosol 

« i 

J droplets  or  particles 

Smoke 

« i 

< O.lum 

Suspension 
of  dust 

« i 

> 0.1;*m 

Fluidized  dust 

- 1 

l^m 

For  high  volume  concentrations,  the  action  of  the  condensed  phase  on  the  gaseous  phase,  and  interactions 
between  particles,  are  important,  so  that  we  are  led  to  analytic  processes  and  to  experimentations  quite 
different  from  those  concerning  low  concentrations  (see  also  E 3.2,  4 and  7). 

For  low  volume  concentrations,  to  which  we  shall  limit  ourselves,  we  must  first  consider  the  mean  di- 
mensions of  the  particles:  it  is  known  that  the  action  of  the  gaseous  phase  on  the  condensed  phase  depends 

strongly  on  this  parameter. 


The  flow  of  a gas-particle  mixture  is  controlled  by  the  exchanges  occurring  between  the  two  phases: 
exchange  of  momentum  due  to  mechanical  entrainment  of  the  particles,  exchange  of  heat  resulting  from  thermal 


nonequilibrium  between  the  phases,  and  exchange  of  mass  if  there  is  evaporation,  condensation  or  com- 
bustion. 

Problems  related  to  mass  transfer  between  phases  will  not  be  considered  here:  we  shall  only  consider 

problems  where  the  condensed  phase  Is  driven  only  by  aerodynamic  mechanisms.  We  shall  not  limit  ourselves 
to  solid  particles,  as  many  applications  concern  liquid  particles,  these  two  classes  of  problems  being 
often  treated  by  similar  means. 

1.2  Outlook  of  National  Work 


There  does  not  actually  exist  in  France  any  university  laboratory  specialized  in  the  field  of  gas 
flows  with  a low  concentration  of  particles.  This  type  of  flow  is  however  studied  In  applied  fields  at  the 
Paris  VI  University  and  at  the  Atmospheric  Physics  and  Pollution  Laboratory  (Prof.  Brlcart).  Ihe  National 
Center  of  Scientific  Research  (CNRS)  does  not  conduct  any  specific  research  in  this  field;  in  its  activity 
report,  multiphase  flows  are  mentioned  in  the  category  of  rheological  flows. 

On  the  other  hand,  in  some  institutes  of  applied  research  and  some  industrial  companies,  some  effort 
is  made,  usually  under  the  pressure  of  specific  problems.  This  situation  makes  it  difficult  to  draw  any 
overall  conclusion  regarding  fundamental  studies;  that  is  why  we  shall  start  from 'technical  problems,  and 
then  attempt  to  emphasize  the  fields  in  which  fundamental  knowledge  appears  Insufficient. 

The  dispersion  of  the  lines  of  effort  and  their  technical  specialization  explains  the  small  number  of 
publications  permitting  a study  of  gas-particle  flows.  Professor  Fortier  gathered  in  a single  volume  the 
texts  of  lectures  given  at  the  Paris  VI  University  within  the  framework  of  the  Higher  Hydrodynamics  and 
Aerodynamics  Certificate  (1).  The  general  equations  of  multi-phase  media  are  established,  and  the  laminar 
flow  of  a fluid  around  a spherical  particle  is  derived.  Applications  considered  concern  mainly  turbulent 
flows  In  long  cylindrical  ducts,  with  a view  to  determining  concentration  variations  and  pressure  drop. 

More  recently,  we  summed  up  the  work  initially  oriented  towards  the  operation  of  rocket  engines  (2).  Equa- 
tions for  both  phases  are  written  for  a suspension  with  general  assumptions  on  liquid  or  solid  particles, 
distributed  according  to  their  diameter,  the  mass  exchanges  being  taken  into  account  and  the  gaseous  medium 
being  multi-reactive.  A limited  number  of  applications  is  given  in  the  field  of  aerodynamics. 

The  applied  problems  we  are  going  to  examine  concern  flows  within  solid-propellant  rocket  motors,  those 
encountered  during  implementation  of  laser  anemometry,  and  droplet  accretion  on  obstacles.  All  these  prob- 
lems belong  to  the  aerospace  field,  and  the  mechanisms  brought  into  play  are  varied  enough  to  permit  a 
synthesis  of  a fundamental  nature. 

Before  reviewing  these  applications,  let  us  mention  that  some  work  is  carried  out  in  France  on  high- 
concentration  mixtures,  in  particular  with  a view  to  study  the  transport  of  powdered  matter  (Experiments  1 
Fluid  Mechanics  Laboratory,  Orsay  University,  Prof.  Fortier)  and  the  phenomena  of  sedimentation  and  particle 
removal  from  a wall  (Aerothermodynamics  Laboratory  of  CNRS,  Meudon-Bellevue , Prof.  Bernard). 


2.  FLOW  OF  COMBUSTION  PRODUCTS  IN  SOLID  PROPELLANT  ROCKET  MOTORS 
2.1  Problem  Statement 


Flow  in  rocket  motors  often  present  a two-phase  character,  as  in  the  case  of  liquid-propellant  motors, 
during  combustion  of  propellant  droplets  in  the  combustor;  this  problem,  which  involves  mass  transfer,  lies 
outside  the  scope  of  this  paper.  In  most  solid-propellant  motors,  a dispersed  metallic  charge  is  intro- 
duced into  the  propellant  in  order  to  increase  its  performance.  The  metal  is  found  in  the  combustion  pro- 
ducts in  the  form  of  tiny  oxide  particles  which  may  account  for  up  to  40%  of  the  mass  flow.  Figure  1 shows 
alumina  particles  picked  up  in  the  jet  of  an  experimental  nozzle. 

The  formation  of  this  condensed  phase  is  not  to  be  discussed  here;  we  are  only  concerned  with  the 
aerodynamic  aspect  of  the  phenomena  occurring  in  the  motor  combustor  and  nozzle.  Practically,  the  influence 
of  particles  is  being  felt  in  the  following  characteristics: 

- Performance  - Thermodynamic  calculations  made  in  order  to  predict  performance  assume  that  kinetic 
and  thermal  equilibrium  of  the  phases  is  achieved  locally,  and  thus  do  not  take  into  account  the  two-phase 
aspect  of  the  flow.  The  experimental  specific  Impulse  of  metallized  propellants  is  further  away  from  its 
theoretical  value  than  for  non-metalllzed  propellants;  it  thus  appears  necessary  to  explain  and  predict 
the  influence  of  the  interaction  between  phases  on  performance.  The  loss  of  specific  impulse  seems  highly 
dependent  on  the  operating  parameters  of  the  motor,  in  particular  on  its  scale;  for  small  motors  it  may 
exceed  5%,  which  Justifies  the  attention  paid  to  this  problem.  Figure  2 gives,  as  an  example,  the  varia- 
tions of  the  specific  Impulse  discrepancy  6 I/I  as  a function  of  the  ratio  between  curvature  and  aperture 
radii  at  the  throat  of  a blconlcal  nozzle. 

- Wall  Interference  - This  kind  of  problem  includes  the  accretion  of  pait  of  the  condensed  phase  within 
the  combustor,  when  the  motor  is  rotated  about  its  longitudinal  axis,  formation  of  an  alumina  layer  on  the 
nozzle,  and  erosion  of  part  of  the  nozzle  by  particle  impact.  Figure  3 shows  the  thickness  e of  the  alumina 
layer  deposited  at  the  distance  s from  the  beginning  of  the  convergent  part  of  a small  test  motor  nozzle 

for  zero  and  10  revolutions  per  second. 

- External  Aerodynamics  - Among  problems  related  to  Jets,  those  concerning  their  radiolectrlc  behavior 
and  erosive  effects  are  of  particular  Interest.  Photographs  of  rocket  exhaust  Jets  (3)  show  that  the  con- 
densed phase  tends  to  concentrate  along  the  Jet  centerline,  but  the  available  photographs  are  not  suitable 
for  reproduction  (see  also  E,  Fig.  13). 

- Combustion  Instabilities  - Por  determining  the  stability  range  of  the  motor,  it  is  mandatory  to 
know  the  acoustic  loss  due  to  the  condensed  phase;  this  loss  is  directly  dependent  on  the  mean  particle 
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This  paper  Is  limited  to  problems  concerning  steady  Internal  aerodynamics.  In  the  combustor  and  the 
nozzle. 

2.2  Work  Already  Performed 

The  problem  of  the  flow  of  a two-phase  mixture  within  the  motor  may  be  approached  theoretically  by 
using  general  equations.  However,  we  should  first  give  a sufficiently  fine  description  of  the  condensed 
phase  and,  in  view  of  the  complexity  of  the  calculation,  define  the  simplifications  that  will  permit  a 
faster  processing  of  Che  equations. 

A difficulty  pertaining  to  Che  problem  concerns  the  characteristics  of  the  condensed  phase.  Consider- 
ing the  experimental  results  obtained  (3),  we  may  state  that  the  particles  are  nearly  spherical  and  very 
broadly  distributed  In  size.  Figure  4 gives  the  distribution  of  Che  particle  diameters  dp  obtained  experi- 
mentally by  analysis  of  a sample  of  particles  picked  up  within  a rocket  motor  Jet.  Moreover,  the  alumina 
particles  are  liquid  during  most  of  the  expansion,  and  their  collisions  entail  an  Increase  and  a local  var- 
iation of  their  size  distribution;  this  effect  Is  reinforced  by  particle  segregation  In  zones  with  high 
velocity  gradient,  which  Implies  also  a non-uniform  distribution  of  the  condensed  phase  In  cosibustlon 
products. 

Calculation  emphasizes  that  the  effects  are  about  proportional  to  the  square  of  the  particle  mean 
diameter;  so  It  Is  essential  to  determine  experimentally  the  particle  size  distribution.  Two  techniques 
are  used  in  France:  capture  and  optical  measurement. 

The  capture  technique  was  developed  at  length  at  ONERA.  The  particles  are  captured  mainly  In  the  Jet, 
far  from  the  outlet  section;  this  measurement  gives  an  Idea  of  the  particle  s.ze  distribution  after  solidi- 
fication and  during  the  whole  motor  operating  time  (tests  never  exceeded  10  seconds).  A careful  statistical 
analysis  of  a sample  of  deposited  particles  permits  a calculation  of  the  various  mean  diameters  of  the  dis- 
tribution law. 

Measurements  of  the  same  type  were  a. so  performed  by  capture  within  the  combustor  and  the  nozzle.  In 
spite  of  the  difficulties  due  to  the  liquid  state  of  the  trapped  particles,  the  results  prove  the  reality 
of  particle  size  Increase.  The  capture  technique  Is  rather  difficult  to  handle;  moreover,  the  statistics 
established  this  way  are  relatively  Inaccurate  at  extreme  values  of  the  diameter,  and  the  wchanlsm  of 
particle  deposition  on  capture  plates  Is  not  well  known. 

Another  technique  was  more  recently  developed  by  a Joint  effort  of  SEP  (Soclete  Europeenne  de  Propulsion) 
and  CERT  (the  ONERA-dependent  research  center  at  Ttoulouse);  it  is  based  on  radiation  scattering  by  a particle 
cloud,  and  rests  on  the  classical  Mle  theory  (4).  First  measurements  concerned  a rocket  motor  Jet  and  pro- 
vided encouraging  results,  emphasizing  in  particular  the  evolution  of  the  particle  distribution  during  siotor 
operation.  Figure  5 gives  a general  view  of  the  Instrumentation,  Fig.  6 an  example  of  results.  This  method 
is  being  Improved,  on  the  one  hand  with  regard  to  calibration  with  particles  of  known  size,  on  the  other  by 
extending  the  measuring  range  and  the  processing  method.  It  is  also  envisaged  -o  measure  the  particle  size 
distribution  inside  a motor  equipped  with  a window. 

The  quality  of  aerodynamic  calculations  depends  on  available  data  on  the  size  distribution  of  the  con- 
densed phase.  In  spite  of  the  progress  still  to  be  made  In  this  field,  many  calculation  schemes  have  al- 
ready been  developed.  Apart  from  a few  exceptions  (3),  size  Increase  and  collisions  between  particles  are 
neglected  and,  according  to  the  additional  slmpllficat ions  made,  the  calculations  are  more  or  less  realistic 
(Fig.  7). 

A first  approach  consists  In  studying  the  movement  of  a particle  of  given  size  in  a reference  flow  de- 
fined by  a two-phase  equilibrium  calculation  of  this  type,  performed  by  the  Aerospatiale  Company.  It  shows 
that  the  particle  trajectories  do  not  noticeably  depart  from  the  flow  streamlines  except  In  the  zones  ot 
high-velocity  gradient:  In  the  vicinity  of  a slope  discontinuity  of  the  duct  meridian  of  rhe  nozzle  and 

in  the  nozzle.  More  complete  flow  schemes  are  at  present  being  developed  in  various  companies  working  on 
propulsion  - Aerospatiale,  SEP,  SNPE  (the  National  Society  for  Powders  and  Explosives)  - for  the  study  of 
the  various  Interesting  zones,  mainly  the  nozzle.  Only  a single  particle  size  Is  at  first  considered,  the 
flow  Is  usually  assumed  laminar,  various  laws  for  drag  and  thermal  exchange  are  considered,  and  radiation 
phenomena  are  ignored. 

Experiments  were  also  performed  by  Aerospatiale  by  simulating  the  flow  of  combustion  products  In  a com- 
bustor by  a gas-particle  mixture  and  by  determining  the  drag  on  a Jet-thrust-measurlng  ground  facility  for 
various  nozzle  shapes  and  various  loading  ratios;  these  tests  are  of  fundamental  Interest,  as  the  experi- 
mental conditions  are  well  known,  though  they  do  not  represent  exactly  the  i-eal  phenomena  (5). 

The  study  of  wall  phenomena  Is  directly  linked  to  aerodynamic  calculations,  which  should  permit  the 
prediction  of  impacts  in  the  various  zones  of  the  motor.  Moreover  It  Is  necessary  to  predict  the  behavior 
of  particles  trapped  in  the  boundary  layer  In  order  to  determine  the  deposits.  These  problems  seldom  have 
been  approached,  and  then  only  from  very  simplified  capture  calculations  In  which  the  trajectory  Is  being 
studied  In  a reference  gaseous  flew. 

2.3  Problems  Still  to  be  Solved 

The  main  cause  of  uncertainty  Is  the  Incomplete  knowledge  of  size  measurement  of  the  condensed  phase 
and  the  mechanisms  of  particle  collisions  and  accretion.  The  calculation  schemes  must  take  these  phenomena 
Into  account  to  be  realistic.  The  role  of  turbulence  on  the  condensed  phase  should  be  clarified.  So  It 
seems  that,  for  the  moment,  the  research  work  should  concentrate  more  on  the  development  of  measuring 
methods  and  on  fundamental  studies  than  on  the  Improvement  of  calculation  methods. 
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3.  LASER  ANEHOCTOY 


Problem  Statement 


Among  aerodynamic  measuring  methods,  the  optical  ones  are  particularly  lntereatlng  as  they  present  the 
advantages  of  not  disturbing  the  flow,  of  being  usable  In  difficult  conditions  and  of  providing  nearly 
point  measurements.  Among  them,  laser  anemossetry  was  the  subject  of  auny  Improvements  during  the  last 
years.  Hie  principle  Is  nw  well  known  and  can  be  explained  either  by  a calculation  of  the  Doppler  frequency 
of  the  light  scattered  by  the  moving  particles,  or  by  the  crossing  of  the  particles  through  the  Interference 
fringes  created  by  Interaction  of  two  Incident  laser  beams.  The  analysis  of  the  frequency  of  the  aignal 
received  by  the  decector  provides  the  velocity  of  the  particles  entrained  by  the  flow. 

In  Its  practical  utilization  In  aerodynamics,  this  method  raises  a problem  of  Interpretation,  aa  we 
have  to  deduce  the  fluid  velocity  from  the  particle  velocity.  These  two  velocities  sre  the  same  If  the 
particles  are  small  enough  to  follow  Instantaneously  the  velocity  tluctue cions  of  the  gaseous  phase.  It 
Is  thus  important  to  define  precisely  the  tolerable  particle  size,  or  to  establish  the  corrections  permit- 
ting the  passage  from  the  measured  velocity  to  the  fluid  velocity  (see  also  C.7,  D.S,  and  E.9). 

With  powerful  enough  lasers  It  Is  possible  to  v;rV  with  particles  that  are  normally  In  suspension  In 
any  flow.  Injection  of  smoke  can  also  be  envisaged. 

3.2  Work  Already  Performed 

Laser  anemometry  was  simultaneously  developed  at  the  Sa'.nt  Louis  r-anco -German  Institute  (I.S.L.)  (6) 
and  at  ONERA  (7);  good  experiments  were  performed  with  both  types  of  lnatnasentatlon.  The  two  crossed-beam 
instruments  are  different  with  regard  to  signal  analysis,  which  Is  performed  by  counters  at  ISL  and  by  a 
spectrum  analysis  at  ONERA.  The  ISL  inatrtaent  permits  the  elimination  of  spurious  signals  and  seems  to 
be  better  adapted  to  the  measurement  of  rather  low  turbulence  rates  during  relatively  short  times.  The 
ONERA  Instrument  offers  a better  space  resolution  thanks  to  a smaller  measuring  volume;  high  turbulence 
rates  can  be  measured,  but  only  over  relatively  long  clams. 

The  first  ISL  experiments  concerned  the  turbulence  rate  In  a small  free  jet  at  low  speed;  the  rates 
measured  for  the  axial  velocity  fluctuations  In  the  jet  agree  well  with  determinacicns  done  with  the  hoc 
wire  technique-,  as  seen  in  Fig.  8.  The  velocity  field  atound  a two-dimensional  wing  profile  was  then  es- 
tablished, but  without  comparison  with  other  methods,  lr>  order  to  determine  the  possibilities  of  the  Instru- 
ment; measurements  of  the  two  velocity  components  In  the  leading-edge  area  were  obtained.  Results  of  the 
exploration  of  a thin  boundary  layer  compare  favorably  with  the  Blasius  theoretical  profile.  During  these 
various  experiments,  T**L  made  an  effort  to  know  the  size  distribution  of  the  particles  and  natural  dust  In 
suspension  In  the  experimental  Installations:  the  ms  In  method  used  Is  that  of  very  fine  sieves,  obtained 

by  chemical  treatment  of  brass  whiskers,  able  to  remove  particles  1 urn  in  diameter;  it  is  not  however  sure 
that  the  biggest  particles  are  collected  by  the  sieve,  and  that  they  are  not  responsible  for  the  difficul- 
ties encountered  In  the  study  of  low-v. loclty  flows. 

The  ONERA  techniques  were  applied  to  the  study  of  a flow  behind  a step  In  the  S8B  wind  tunnel.  The 
mean  velocity  profile  at  mid-step  compares  favorably  with  pitot  amasurements  and  with  those  performed  with 
a hot  wire  In  a similar  Incompressible  flow.  The  turbulence  rate  Is  estsbllshed  up  to  a high  level,  and 
corresponds  approximately  to  that  provided  by  the  hot-wire  method;  the  differences  observed  between  these 
reo  types  of  measurements  are  found  In  the  turbulence  profile  , as  seen  In  Fig.  9.  The  mean  velocity  pro- 
files and  the  turbulence  rate  profiles  In  the  boundary  layer  of  the  supersonic  flow  In  front  of  the  step 
were  also  recorded;  results  for  mean  velocity  are  In  good  agreement  with  the  pitot  meestrements.  During 
the  experiments  the  air  used  In  the  wind  tunnel  was  thoroughly  filtered  before  Its  passage  through  the 
test  section.  A check  on  dust  size  distribution  by  the  whisker-capture  method  confirmed  that  the  particles 
were  of  submicron  size,  but  the  quality  of  Ghe  comparisons  did  not  justify  further  Investigation  by  this 
technique. 

before  the  development  of  lase’-  anemometry,  ONERA  had  worked  on  photographic  techniques  (6).  The  study 
concerned  Internal  aerodynamics  In  the  -nd  part  of  a solid-propellant  motor,  and  used  a two-dimensional 
model,  fed  with  compressed  gas  through  porous  walls  and  with  20  ^m-dla.  particles  at  some  points  of  the 
walls.  Illumination  of  the  test  section  by  a continuous  laser  permitted  the  determination  of  the  mean 
particle  trajectory  (Fig.  10);  this  trajectory  la  well  defined  for  20  pm  particles,  but  we  could  detect  a 
rapid  dispersion  of  1 pm  particles  after  their  emission  at  the  wall.  Illumination  of  the  teat  section  by 
a pulsed  laser,  either  In  Its  natural  or  g-svltched  mode,  permits  In  principle  the  visualization  of  the 
successive  particle  positions,  and  thus  the  calculation  of  the  velocity  components;  the  values  obtained  for 
1 pm  particles  are  highly  dispersed,  which  may  be  explained  by  flow  turbulence. 

3.3  Problems  Still  to  be  Solved 


In  the  application  of  laser  anemometry  to  classical  aerodynamic  studies,  the  flows  Include  very  few 
particles  In  suspension.  The  action  of  the  condensed  phase  on  the  flow  is  negligible,  and  the  main  problem 
Is  to  determine  the  limit  size  under  which  the  particles  follow  the  flow  fluctuations  without  delay.  In 
this  type  of  experimentation,  we  partially  control  the  particle  size  distribution,  for  Instance  by  thorough 
filtering  of  the  gas;  water-vapor  condensation  has  to  be  avoided.  A check  of  particle  size  can  be  made 
either  by  capture  or  by  optical  measurements,  but  the  quality  of  the  results  obtained  Is  more  generally 
established  by  comparison  with  other  physical  measurements,  such  as  hot  wire,  If  the  velocity  is  sufficient, 
or  pressure.  That  is  why  the  aerodynamic  studl.  a conducted  at  ISL  and  ONERA  have  not  yet  shown  the  need 
to  go  further  with  problems  specific  to  the  condensed  phase. 

If  we  attempt  to  apply  laser  anemometry  to  flows  with  noderflte  oartlcle  content,  or  with  panicles  of 
more  than  1 y.m  diameter,  we  must  expect  to  encounter  new  problems  with  regard  to  signal  analysis  and  Inter- 
pretation of  results.  Velocity  measurements  then  cannot  be  separated  from  those  of  particle  size  and  con- 
centration. 


4.  PARTICLE  CAPTURE 


4.1  Problem  Statement 


Particle  capture  was  already  mentioned  at  the  occasion  of  condensed-phase  size  measurement  and  of  flow 
In  rocket  motors.  We  shall  now  examine  activities  related  to  aircraft  Icing  and  to  rain  efrects.  The  flows 
considered  are  lightly  loaded  with  droplets  of  about  20  )im  diameter  for  Icing  and  of  1 mm  diameter  for  rain. 

The  theoretical  study  of  laws  governing  Icing  sheds  some  light  on  several  physical  phenomena.  The 
first  mechanism  concerns  the  capture  itself  of  droplets  by  the  obstacle;  we  then  have  to  determine  their 
movement  within  the  aerodynamic  field  around  the  obstacle,  and  the  similitude  conditions  of  this  problem 
were  established  quite  a long  time  ago  (9).  Hie  second  phenomenon  to  be  studied  is  that  of  water  and  Ice 
distribution  over  the  aerofoil;  lastly  we  have  to  treat  the  problem  of  thermal  equilibrium  of  the  wall  In 
an  Icing  atmosphere.  There  has  been  no  study  in  France  on  the  latter  two  points,  and  the  results  used  are 
American  or  British.  Rather  general  similitude  conditions  are  thus  available,  and  permit  the  reproduction 
in  the  wind  tunnel  of  given  flight  conditions.  The  tests  performed  have  an  Industrial  character  and  aim 
to  be  able  to  fix  experimentally  the  Icing  conditions  at  the  given  values. 

With  regard  to  rain  erosion,  we  are  mainly  concerned  with  capture  conditions,  and  results  remain  quali- 
tative. Again,  the  point  Is  to  be  able  to  reproduce  given  physical  characteristics. 

4.2  Work  Already  Performed 


Icing  and  rain  erosion  are  studied  at  ONERA  at  the  Modane  Wind  Tunnel  Department.  Icing  studies  were 
recently  reactivated  for  new  French  or  European  aircraft.  In  particular  the  Concorde  was  the  occasion  of 
an  important  series  of  tests,  following  doubts  expressed  abroad  (10).  Experimentation  was  carried  out  in 
the  large  Sl-MA  wind  tunnel  on  a 1/6  scale  model  (Fig.  11).  The  application  domain  of  similitude  conditions 
was  extended  from  1/6  to  1/12  scale  by  studies  of  a more  fundamental  character  on  simpler  geometries  (cyl- 
inders, classical  delta  wings).  The  results  obtained  made  it  possible  to  qualify  satisfactorily  the  Icing 
similitude  technique.  Icing  problems  on  aircraft  engines  are  also  carried  out  at  the  French  Engine  Test 
Center  (CEPr). 

High-speed  (Mach  2)  rain  erosion  tests  are  performed  at  the  Landes  Test  Center  (CEL)  on  aeronautical 
structures  and  on  missiles  (11).  The  elements  under  study  are  placed  on  a trolley  moving  on  a single  or 
dual-rail  track,  600  and  400  m long,  and  accelerated  by  solid  propellant  motors  In  a cloud  of  droplets 
created  by  ejectors  (Fig.  12).  Approximate  similitude  conditions  are  used  to  reprodut  ■ damage  observed  in 
flight;  parameters  to  be  adjusted  are  exposure  time,  rain  intensity,  droplet  size  and  trolley  speed.  This 
installation  Is  operational,  and  work  is  in  progress  to  Increase  the  velocity  range. 

4.3  Problems  Still  to  be  Solved 

The  main  practical  problems  concern  droplet  production  and  dimensional  control  to  make  sure  that  the 
test  conditions  are  satisfied. 

In  the  ONERA  Icing  installations,  the  droplets  are  produced  on  a atomization  grid  by  air  break-up 
action;  the  Influence  of  adjustable  parameters  (water  flow  rate,  air  pressure)  was  determined  empirically. 
The  artlflclal-rain  generating  station  at  CEL  was  made  by  the  Bertin  Company,  under  the  direction  of  the 
official  Launching  Range  Instrumentation  Service  (SECT),  and  Is  made  of  atomization  manifolds;  drop  size 
is  controlled  by  water  pressure. 

Qualification  of  mists  msde  this  way  appeared  as  a difficult  problem.  The  measuring  method  used  for 
rain  erosion  Is  that  of  drop  capture  on  an  oily  plate;  the  photograph  of  the  sample  is  then  analyzed  by 
traditional  means  of  counting  by  size  classes.  Icing  was  the  objective  of  various  measurements;  the  tech- 
nique of  capture  on  an  oily  plate,  also  implemented  at  CEPr,  was  first  employed  and  based  on  a statistical 
analysis,  using  a catalogue  developed  by  the  British  firm  Hawker-Si dde ley ; precision  of  the  results  is 
limited  by  the  various  assumptions  which  have  to  be  made  regarding  the  spherical  shape  of  the  collected 
droplets,  absence  of  droplet  accretion  ana  evaporation,  nonmiscibility  of  oil  and  water,  and  existence  of 
a known  ratio  between  the  diameter  of  the  observed  image  and  that  of  the  actual  droplet.  A photometric 
method  was  tried,  but  the  instrument,  initially  developed  by  the  ESRI  Company  and  improved  with  regard  to 
optics  and  electronics,  did  not  give  satisfactory  results.  An  instrument  developed  by  the  Bertin  Company, 
using  droplet  capture  in  a softened  laver,  measured  imprints  left  after  evaporation,  but  these  measurements 
do  not  agree  well  with  those  obtained  with  oily  plates.  Thus,  there  remains  an  uncertainty  about  the  size 
analysis  of  the  mists  produced,  and  new  measuring  methods  are  currently  investigated. 

Icing  and  rain-erosion  tests  of  industrial  character  do  not  seem  to  raise  specific  two-phase  aerody- 
namic problems.  The  influence  of  some  phenomena  related  to  liquid  particles  (deformation,  accretion)  is 
ignored.  On  the  other  hand,  progress  in  techniques  for  measuring  mist  sizes  would  be  highly  appreciated 
(see  also  B 3.5,  C.4  and  5,  and  E.9). 


5.  SUMMING  UP  THE  PROBLEMS 

We  shall  not  attempt  to  present  a synthesis  of  the  problems  that  appeared  during  the  surv^  of  previous 
applications.  Knowledge  of  si.e  distribution  is  always  necessary;  this  problem  being  assumed  solved,  the 
incidence  of  mechanisms  related  to  particles  on  the  results  seems  to  depend  on  the  application  considered; 
that  is  why  we  shall  only  enumerate  the  tie  Ids  where  progress  in  our  knowledge  would  be  useful. 

5,1  Knowledge  of  the  Condensed  Phase 

There  exist  many  classical  means  for  measuring  the  size  distribution  of  a condensed  phase;  let  us  men- 
tion, for  the  record,  filters  (large  particles),  diffusion  batteries  for  pollution  (fine  particles  or  aero- 
sols), suspension  of  particles  in  an  electrically  conductive  liquid  (Coulter  counter),  separation  devices 
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(cyclones).  These  methods  are  not  applicable  to  the  flows  considered,  which  usually  contain  a small  amount 
of  particles  broadly  distributed  in  size.  The  only  methods  used  are  either  direct  capture  methods  or  In- 
direct ones  based  on  the  optical  properties  of  the  suspension. 

5.1.1  Particle  Capture 

According  to  the  nature  of  particles,  the  capture  element  may  vary  widely:  whiskers  (particles  in 

suspension  wind  tunnels  — ISL) , glass  or  quartz  plates  (alumina  particles  of  combustion  product*-  of  metal- 
lized propellants  — ONERA),  oily  plates  or  plates  coated  with  a plastic  film  (liquid  particles  — Berlin 
Co.).  In  all  cases,  exposure  times  should  be  strictly  controlled. 

Statistical  analysis  of  the  sample  requires  powerful  enlarging  means,  when  particles  are  submicroscoplc 
(use  of  a scanning  electron  microscope  — ONERA).  Automatic  processing  of  the  image  Is  often  out  of  the 
question  in  view  of  particle  superposition,  and  manual  processing  Is  always  long  and  difficult. 

Interpretation  of  results  is  possible  only  If  the  movement  of  particles  in  the  aerodynamic  field  is 
known  and  the  capture  efficiency  calculated.  This  problem  requires  the  knowledge  of  the  drag  laws,  little 
studied  for  particles  of  complex  shape  or  for  clusters  (soot).  When  the  particles  are  small  and  the  flow 
turbulent,  the  deposition  mechanism  is  poorly  understood. 

It  is  desirable  to  calibrate  the  collecting  devices  in  conditions  near  those  of  actual  utilization. 

We  then  must  have  at  our  disposal  particles  of  given  shape  and  size  distribution,  and  possibly  of  uniform 
size;  these  materials  are  still  not  available  in  sufficient  quantities,  especially  for  microscopic  sizes. 

The  capture  technique  was  for  a long  time  the  only  one  available,  and  gave  valuable  information  in 
most  applications.  It  is  however  limited  with  regard  to  precision  for  extreme  size  values,  in  particular 
for  very  broad  distributions.  In  spite  of  progress  accomplished,  it  remains  rather  difficult  to  handle 
and  provides  only  mean  measurements  in  time  and/or  space. 

5.1.2  Optical  Methods 

We  can  count  among  these  methods  photographic  techniques,  holography  and  processes  making  use  of  light 
transmission  or  scattering,  in  particular  by  a cloud.  Photographic  techniques  are,  for  example,  used  to 
determine  the  size  of  droplets  formed  by  an  injector;  their  use  remains  limited  to  relatively  large  dia- 
meters not  requiring  too  big  an  enlargement  which  is  incompatible  with  a sufficient  depth  of  field. 

Holography  was  applied  to  the  study  of  an  icing  cloud  at  ISL,  for  application  to  test  chambers  at 
CBPr  (12)  (Fig.  13).  Results  are  promising  regarding  the  particular  problem  envisaged.  The  field  of  appli- 
cation of  this  technique  seems  however  limited  to  particle  diameters  of  more  than  5 pm  and  low  concentra- 
tion if  the  hologram  quality  is  not  to  be  altered,  and  at  velocities  no  higher  than  100  m/sec.  Processing 
is  relatively  slow,  as  for  capture.  The  hologrcphic  method  is  thus  only  interesting  for  certain  fields 
of  application,  and  progress  remains  to  be  needed  to  make  its  practical  application  easier. 

Processes  making  use  of  transmission  and  scattering  properties  of  a suspension  were  the  object  of 
quite  a few  applications  in  ln&istrlal  fields,  which  should  be  mentioned  in  spite  of  their  non-aerodynamic 
character.  Measurements  of  the  optical  absorption  coefficient  of  a mist  were  performed  at  the  Research 
Division  of  EDF  (the  French  Electricity  Board)  with  a view  to  study  two-phase  phenomena  in  steam  turbines; 
preliminary  studies  concerned  the  determination  of  mean  diameters  in  the  0.02-10  pm  range  in  various  ex- 
perimental conditions  (13).  Applications  are  under  way  or  envisaged  at  the  Hydraulic  Binders  Industrial 
Research  Center  (continuous  control  of  cement  mixers)  and  at  the  Coal  Board  Research  Center  (measurement 
of  the  dust  content  in  mine  galleries). 

The  study  of  light  scattering  by  a particle  cloud  should  permit  a fine  particle-size  determination, 
and  was  the  object  of  an  application  in  chemical  propulsion.  First  results  shed  light  on  several  classes 
of  problems,  which  limit  the  accuracy  of  the  results,  lhe  first  category  of  problems  concerns  the  analysis 
of  scattering;  whereas  the  theory  for  single  scattering  cannot  be  questioned,  we  can  express  doubts  on  the 
role  of  multiple  scattering  in  case  of  relatively  high  particle  concentrations  and  on  that  of  non-homogeneity 
of  the  spatial  distribution  of  the  condensed  phase  in  the  flow.  The  diffraction  coefficients  of  the  par- 
ticle material  should  be  precisely  known,  which  is  not  alvays  possible  for  materials  formed  at  high  temper- 
ature, Optical  characteristics  of  the  set-up  should  be  thoroughly  checked,  as  they  control  the  useful  meas- 
uring range.  Lastly,  deriving  size  distribution  from  measurements  is  an  arduous  numerical  problem.  Pro- 
gress may  be  expected  by  making  for  instance,  simultaneous  measurements  of  scattering  and  attenuation,  or 
by  working  with  several  wavelengths.  Thus,  this  method  is  in  full  development  and  it  is  still  too  early 
to  define  its  limits. 

5.1.3  Synthesis 

There  is  no  universal  method  for  determining  size  distribution  of  a condensed  phase,  in  view  of  the 
extreme  variety  of  the  particles,  of  their  size  range  and  of  their  concentration  in  the  flow.  All  exist- 
ing methods  are  rather  difficult  to  handle  and  do  not  offer  the  expected  accuracy.  In  these  conditions, 
capture  and  optical  processes  should  be  considered  for  the  moment  as  complementary  means.  Only  optical 
methods  present  an  interesting  potential  for  development;  progress  to  be  expected  in  the  knowledge  of  the 
condensed  phase  rests  on  the  research  efforts  in  the  field  of  applied  optics. 

5.2  Velocity  Measurements 

Regarding  velocity  determination,  the  two  main  problems  concern  the  measurement  itself  and  its  use 
in  aerodynamic  studies. 

The  only  methods  envisaged  here  are  the  optical  methods  developed  in  classical  aerodynamics.  Photo- 
graphic techniques,  an  example  of  which  is  given  in  Section  3,  seem  to  be  the  only  ones  applicable  to  par- 
ticles with  a diameter  of  the  order  of  10  ^m  and  relatively  low  concentration  in  the  flow:  laser  anemometry 
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requires  very  low  concentrations  for  signal  analysis. 

Regarding  interpretation,  for  tests  performed  for  the  aerodynamic  study  of  homogeneous  flows  in  which 
the  particles  are  being  used  only  as  tracers,  perturbation  by  the  condensed  phase  is  to  be  avoided;  the 
point  is,  in  this  case,  to  know  if  the  particles  are  small  enough  to  follow  all  the  fluctuations  of  the 
gaseous  phase.  Ihis  problem  is  theoretically  solved  for  spherical  particles  and  a turbulent,  homogeneous 
and  isotropic  flow. 

The  study  of  real  two-phase  flows,  i.e.,  with  a particle  content  sufficient  to  act  on  the  gaseous  flow, 
has  not  been  experimentally  approached.  The  application  of  laser  anemometry  to  mix  tures  with  a high  con- 
centration of  particles  will  probably  be  difficult,  and  the  interpretation  of  measurements  will  only  be  ob- 
tained if  the  local  concentration  of  the  condensed  phase  is  determined  and  if  the  various  forces  acting  on 
the  particles  are  known. 

The  drag  laws  currently  used  in  France  are  taken  from  the  literature  and  are  convenient  for  most  appli- 
cations where  particles  are  solid  and  spherical.  Complementary  fundamental  research  would  be  useful  in 
the  following  fields: 

— behavior  of  an  isolated  particle  in  a viscous  fluid  flow  with  velocity  gradient  (Aerothermics 
Laboratory,  CNRS), 

--  drag  laws  for  liquid  particles,  taking  into  account  deformations  due  to  high  velocity  gradients, 
and  Internal  circulation, 

--  effects  of  the  surface  state  on  the  solid  particle  drag, 

--  study  of  unsteady  entrainment,  with  a view  to  applications  to  turbulent  flows, 

— drag  laws  of  non-spherical  particles  and  of  clusters. 

Velocity  measurements  were,  up  to  now,  only  treated  for  practically  homogeneous  flows.  An  important 
effort  remains  to  be  made  for  two-phase  flows  in  both  the  experimental  and  fundamental  fields  (see  also 
C.8.1). 


5.3  Temperature  Measurements 

Temperature  measurements  are  especially  interesting  for  very  hot  flows;  their  interpretation  raises  a 
particular  problem  when  concentration  and  temperature  gradient  are  high  enough  for  differences  to  appear 
between  the  temperatures  of  the  two  phases.  In  particular,  this  is  the  case  in  the  expansion  of  combustion 
products  in  rocket  motors. 

Gas  temperature  Tg  and  particle  mean  temperature  Tp  have  been  measured  at  ONERA  in  a rocket  Jet  (14). 
Gas  temperature  is  determined  by  absorption  of  two  monochromatic  radiations  in  the  near  infrarej  range. 
Temperature  of  the  condensed  phase  is  deduced  from  the  Jet  emission  characteristics,  assumptions  having 
to  be  made  regarding  homogeneity  of  the  particle  distribution  within  the  flow  and  the  particle  size  dis- 
tribution. Figure  14  gives  an  example  of  the  pressure  p and  particle  temperature  obtained  as  functions 
of  the  test  time  t;  temperature  differences  beLween  the  phases  are  emphasized. 

Analysis  of  the  phenomena  implies  the  knowledge  of  thermal  exchange  laws  between  gas  and  particles. 

The  main  laws  used  result  from  a compilation  of  the  literature,  and  complements  should  be  established  for 
unsteady  mechanisms  and  non-spherical  particles. 

Temperature  and  velocity  measurement  problems  are  different  as  it  is  possible  to  measure  the  gas  tem- 
perature directly.  The  Interest  in  temperature  measurement  is  only  obvious  in  particular  applications. 

Let  us  lastly  remark  that  the  influence  of  thermal  phenomena  on  two-phase  flows  is  usually  smaller  than 
that  of  dynamic  phenomena  (see  also  E.6.1). 

5.4  Interaction  Between  Particles 


This  problem  concerns  only  two-phase  flows  with  rather  high  particle  content.  They  are  however  of 
g ■«:  practical  interest,  in  particular  for  two  of  the  applications  previously  mentioned:  rocket  motor 

internal  aerodynamics  and  droplet  capture. 

When  the  particles  are  solid,  their  collision  tends,  on  the  one  hand,  to  limit  the  velocity  differ- 
ences between  particles  (see  also  E,  Fig.  8)  and,  on  the  other,  to  ensure  their  diffusion  within  the  flow. 
From  a theoretical  point  of  view,  the  movement  of  each  class  of  particle  should  be  studied  by  examining 
the  momentum  exchanges  between  that  class  and  all  other  particles.  This  interaction  term  is  established 
by  a statistical  analysis  of  the  impacts,  and  this  calculation  implies  a good  knowledge  of  the  elementary 
collision  mechanism,  a knowledge  not  yet  acquired  in  a sufficiently  general  manner. 

The  problem  becomes  more  difficult,  but  also  more  important,  when  we  consider  liquid  particles.  Col- 
lisions entail  accretions  and  possibly  breakdowns,  appearing  as  a continuous  evolution  of  the  condensed 
phase,  particularly  unfavorable  for  aerodynamic  predictions.  The  elementary  mechanisms  of  coalescence  are 
little  known,  and  we  can  predict  neither  the  collision  efficiency  nor  the  size  of  the  particles  formed 
when  accretion  is  followed  by  breakdown. 

Interactions  between  particles  seem  to  have  to  be  studied  first  at  the  level  of  elementary  phenomena, 
which  will  probably  entail  implementation  of  new  facilities.  Macroscopic  equations  for  both  phases  will 
then  have  to  be  written  (see  also  E.5).  Radiative  exchanges,  in  particular,  that  have  not  been  mentioned, 
should  be  predicted  with  the  available  optical  theories  (see  also  B.2.4  and  E.3,3). 


6.  CONCLUSIONS 


Aerodynamic  problems  of  gas-particle  mixtures  were  mainly  approached  in  France  in  connection  with 
particular  applications,  which  explains  the  relative  paucity  of  works  of  fundamental  character  and  the 
wide  diversity  of  the  techniques  implemented,  theoretical  approaches  remained  rather  general,  so  tfa  t 
fundamental  phenomena  are  imperfectly  known.  Flow  analyses  should  be  complemented  in  the  fields  already 
studied,  and  expanded  to  high-concentration  suspensions.  Progress  will  depend  on  development  of  experi- 
mental means,  and  mainly  of  optical  methods  whose  performance  should  be  improved  and  implementation  sim- 
plified. 

Applications  present  a variable  degree  of  success  depending  on  the  nature  of  the  particles,  their  size 
distribution  and  their  concentration;  the  most  difficult  problems  are  raised  by  liquid  particles,  broad 
distributions,  high  concentrations  and  high  temperatures.  In  the  present  situation,  new  applications  seem 
to  appear  in  the  field  of  energy  exchanges:  to  raise  the  efficiency  of  thermal  machines  and  of  transfor- 

mation cycles,  fluids  are  sought  that  possess  extreme  physical  properties,  in  particular  with  regard  to 
transformation  temperatures  and  heats.  We  could  thus  imagine  using  "two-phase  fluids"  made  for  example 
of  capsules  of  coated  mineral  salts  in  suspension  within  a fluid  with  low  freezing  point  and  low  viscosity. 
It  will  be  possible  to  optimize  these  fluids  with  a composite  structure,  whose  behavior  will  be  controlled 
at  all  levels  of  the  cycle  by  the  laws  of  two-phase  flows. 
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SUMMARY 

A general  review  of  German  activities  in  the  field  of  gas  flows  with  solid  particles  is  given.  Fun- 
damental research  on  force  and  heat-transfer  laws  and  propagation  of  sound  and  shock  waves  is  discussed 
as  well  as  more  applied  work  on  generation,  conveying,  separation  and  measurement  of  particles. 


1.  INTRODUCTION 

1.1  Technical  Motivation  for  Research  in  the  Field  of  Gas  Flows  with  Solid  Particles 


Two-phase  flows  are  encountered  in  many  engineering  problems  and  in  nature.  In  the  following  only 
mixtures  of  a gas  (or  a liquid)  and  solid  particles  are  considered.  Such  flow  problems  comprise  the 
conveying  of  solid  particles,  dcsty  gas  flow,  separation  of  solid  particles,  measurement  techniques,  and 
the  sampling  of  micrometeorites.  For  the  understanding  of  such  processes,  the  particle  paths  in  a Held 
of  viscous  forces,  the  conductive  and  radiative  heat  exchange,  the  attenuation  of  sound  and  shock  wa-es 
must  be  known.  Fundamental  research  tends  to  clear  up  these  problems  on  single  particles  and  on  partic'e 
clusters.  Including  Interaction  effects.  In  the  following,  activities  in  the  mentioned  fields  are  re- 
viewed, which  have  been  developed  in  Germany  during  recent  years. 


2.  FUNDAMENTAL  RESEARCH 

2.1  Forces  Acting  Upon  Microscopic  Particles 

The  drag  of  small  particles  depends  on  Reynolds  number,  Mach  number,  and  heat-transfer  rate.  It  is 
rather  accurately  known  only  for  limiting  cases  as  discussed  by  Wuest  (1).  For  high  Mach  numbers,  the 
Newtonian  drag  law  with  a quadratic  dependence  on  velocity  is  valid,  for  high  Reynolds  numbers,  too,  a 

quadratic  law  describes  the  turbulent  regime.  For  low  Reynolds  numbers,  however,  the  drag  is  Inversely 

proportional  to  Reynolds  number  in  the  Stokes  regime,  but  inversely  proportional  to  the  Mach  number  in  the 
subsonic  Knudsen  regime.  The  limit  between  both  regimes  is  roughly  given  by  the  condition  that  the  particle 
diameter  is  less  than  5 mean  free  paths.  Under  atmospheric  conditions  this  corresponds  to  a particle  size 
of  0.3  pm.  Of  course,  there  are  smooth  transitions  between  the  different  regimes.  For  the  transitional 
regime  between  free-molecular  'low  and  continuum  flow,  measurements  on  uncooled  spherical  particles  were 
undertaken  by  Legge  and  Koppenvallner  (2)  (Fig.  1).  In  the  Knudsen  flow  regime  there  is  a smooth  transi- 
tion from  the  linear  to  the  quadratic  law  with  increasing  molecular  Mach  number  S (Fig.  2). The  molecular 

Mach  number  is  defined  as  the  ratio  S of  the  mass  velocity  to  the  most  probable  velocity  ./2RT,  so  that 
S = J\/2  M,  where  R is  the  gas  constant,  T the  temperature,  Y the  ratio  of  the  specific  heats  and  M the 
Mach  number.  A characteristic  feature  of  this  regime  is  the  dependence  of  the  drag  of  spherical  particles 
on  surface  temperature.  New  measurements  on  cone  drag  in  free-molecular  flow  by  H.  Legge  (3),  which  are 
still  in  progress  (1974),  have  quantitatively  confirmed  the  theoretical  predictions  of  the  dependence  of 
free-molecular  drag  on  surface  temperature.  No  experimental  or  theoretical  data  are  as  yet  available  for 
the  combined  effects  of  velocity  and  surface  temperature  in  the  transitional  regime  (see  also  E.3.1). 

2.2  Particle  Paths  in  a Field  of  Viscous  and  Centrifugal  Forces 

The  density  of  solid  particles  is  generally  much  larger  than  the  density  of  the  ambient  gas.  Because 
of  Inertial  effects,  particles  cannot  follow  immediately  the  motion  of  the  gas  but  show  relaxation.  The 
same  is  true  also  for  temperature  changes.  As  was  mentioned,  in  the  free  molecular  regime  drag  and  surface 
temperature  are  coupled,  and  therefore  velocity  and  temperature  relaxation  must  be  comsldered  together. 

The  same  problem  was  treated  by  Schmltt-v.  Schubert  (4)using  a Laplace  transformation.  In  the  physical 
plane  the  shape  of  the  frozen  shock  was  given  analytically.  For  small  mass  fraction  of  the  particles, 
analytical  solutions  were  obtained  for  all  flow  variables  on  the  surface  of  a wedge.  For  small  ratios  of 
solid  mass  to  gaseous  moss,  the  particle  path  was  considered  by  Wuest  (1)  for  linear  and  quadratic  drag  laws 
In  supersonic  wedge  and  stagnation  point  flow. 

The  total  changes  of  particle  velocity  v and  temperature  Tp  are  ,;iven  by  the  following  formulas, 
where  vg  and  Tg  are  the  gas  velocity  and  temperature,  and  Pp  is  the  density  of  the  particle  material, 
p the  pressure,  d the  particle  diameter,  p and  k the  viscosity  and  thermal  conductivity  of  the  gas,  cD 
the  drag  coefficient,  Re  and  Nu  the  particle  Reynolds  number  and  Nusselt  number,  and  cp  the  specific  heat 
of  the  particles. 
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are  Che  relaxation  times  for  velocity  and  temperature,  and  all  velocities  must  be  Interpreted  as  vectors. 
For  low  particle  Reynolds  number,  the  action  of  pressure  may  be  neglected.  This  is  valid  even  for  the 
passage  of  a particle  through  a shock  wave,  as  was  shown  by  Wuest  (1)  (see  also  C.2.1  and  E.3). 

Perturbed  two-phase  cylindrical-type  flows  and  two-phase  convex-type  corner  flows  were  calculated 
by  Healy  (5,6)  (see  also  E.9).  The  particle  paths  In  rotating  flows  with  additional  axial  flow  were 
Investigated  by  Muhle  (7). 


2.3  Particle  Clusters 

The  motion  of  particle  "lusters  has  been  described  by  Brauer  and  Thiele  (8).  If  a cluster  Is  com- 
posed of  one  species  only,  it  Is  called  a monocluster;  in  the  case  of  several  species,  it  Is  called  a 
polycluster.  In  the  thesis  by  Kaskas  (Technical  University  Berlin  1970)  a formula  for  the  velocity  of  one 
fraction  of  a polycluster  was  derived.  The  formulae  have  been  applied  especially  to  sedimentation 
processes  under  the  Influence  of  gravity,  and  it  was  observed  that  the  sinking  of  larger  particles  may 
produce  a rising  of  smaller  particles  caused  by  the  reverse  flow. 


2.4  Thermal  Exchanges 

For  very  small  particles  a complete  description  of  heat  transfer  Is  possible  for  arbitrary  velocities. 
For  this  case  the  relaxation  time  for  temperature  was  derived  by  Wuest  (1).  In  the  hypersonic  flow  regime, 
the  total  heat  transfer  to  spheres  is  known  for  a rather  extended  range  of  Reynolds  numbers  as  is  shown 
by  Fig.  3,  where  measurements  in  the  hypersonic  low-density  tunnel  at  Gottingen  are  plotted  and  compared 
with  other  data.  The  heat  transfer  In  gas  flows  with  dust  particle  has  been  investigated  in  the  thesis 
by  Hein  (9). 

Radiant  Interchange  among  suspended  particles  and  Its  effect  on  thermal  relaxation  in  gas-particle 
mixtures  was  studied  by  Schneider  (10)  (see  also  E.3. 3).  Apart  from  the  total  radiative  transfer  due  to 
the  contrlbutlona  of  all  particles,  there  can  also  be  a radiant  Interchange  among  particles  of  different 
sizes  within  the  same  local  region.  The  smaller  particles  follow  the  gas  temperature  more  closely  than 
the  larger  ones  and  this  provides  the  possibility  of  radiative  heat  transfer  between  different  particles. 
From  the  basic  equations  describing  this  process  a linear  integro-dlf ferentlal  equation  was  obtained  for 
small  temperature  disturbances  from  an  equilibrium  state.  The  relative  Importance  of  radiation  and  conduc- 
tion is  characterized  by  the  parameter 

W *°T**an 
RG  ’ TR(rn)  ' nk  rn 

with  the  conductive  relaxation  time  given  by 


C 4nr  k 
n 


and  the  radiative  relaxation  time  given  by 


R 16oT*3ct. 


where  r la  the  reference  radius  (most  probable  particle  radius),  o the  Stefan-Boltzroann  constant,  T* 
the  parElcle  temperature  at  equilibrium,  o an  the  nominal  absorption  cross  section  of  a reference  particle, 
m the  mass  of  an  Individual  particle,  and  c the  specific  heat  of  the  solid  particle  material.  To  get 


an  Idea  of  the  magnitude  of  this  parameter  we  obtain  NR(,  ta  1,  If  T* 


2000  K and  r 


From  a practical  point  of  view,  this  temperature  Is  of  a realistic  order  of  magnitude.  The  analysis 
waa  baaed  on  the  assumption  that  the  particles  are  large  In  comparison  to  the  wave  lengths.  Wien's  law 
yields  a wave  length  of  ^max  ” 1.5  pm  for  maximal  Intensity  at  2000K,  so  that  a particle  radius  of  30 
compiles  well  with  the  assumption. 


In  Fig.  4 the  disturbances  of  the  mean  particle-cloud  temperature  and  of  the  gas  temperature  are 
plotted  as  functions  of  dimensionless  time.  With  Increasing  radiation-conduction  parameter  N„„  (N  • 
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no  radiation;  N * 10  strong  radiation)  the  temperature  range  of  Individual  particles  becomes  narrower. 

For  strong  radiation,  the  temperatures  of  particles  of  different  sizes  tend  to  equalize.  For  the  numerical 
calculation,  a distribution  function  was  used  which  contained  a free  parameter  D to  characterize  the  dis- 
persion which  was  symmetrical  about  the  most  probable  particle  radius. 

The  heat  transfer  by  radiation  in  suspensions  of  solid  particles  in  a gas  flow  was  also  considered 
from  a more  technical  view  point  by  Hein  and  Lowes  (11). 

2.5  Propagation  of  Sound  Waves 

A survey  article  on  wave  propagation  in  suspensions  of  solid  particles  in  gas  flow  was  given  by 
Rudlnger  (12),  who  considered  also  German  contributions  to  this  problem  (see  also  E.6.2).  A thorough  inves- 
tigation of  propagation  of  sound  waves  and  formation  on  shocks  is  given  by  Becker  (13)  and  Schmitt-v. 
Schubert  (14-18).  The  lag  between  the  velocities  and  temperatures  of  gas  and  particles  gives  rise  to  a 
special  class  of  relaxation  phenomena,  and  the  particle  flow  may  serve  as  model  fluid  for  studying  more 
complex  real-gas  relaxation  problems.  For  finite  values  of  the  products  cut  and  of  the  relaxation 
times  and  the  angular  frequency,  Schmitt-v.  Schubert  (18)  obtains  complex  equations  ror  the  frequency- 
dependent  sound  velocity  and  damping.  Simple  formulas  are  obtained  for  small  deviations  from  equilibrium 
and  for  the  limiting  cases  of  the  frequency-independent  sound  velocity  b in  frozen  flow  (^T  ^ -•  ®,  -*  “>)  > 

in  semifrozen  flow  (c^  if  (jut  i -•  0,  tm  ^ — °°>  c2  ff  ibt,  "*  ®,  uffj  -*  0)  end  in  full  equilibrium  (c  with 
cur i " 0UT2  “ 0)  and  the  relation  c < (c.  ^ c^)  < b holds,  that  is,  the  sound  velocity  for  abritrary  relaxa- 
tion parameters  lies  between  the  equilibrium  sound  velocity  c and  the  frozen  sound  velocity  b.  The  sound 
velocity  in  gas  flows  with  solid  particles  was  also  investigated  by  Rumpf  and  Gregor  (19)  and  Weber  (20), 
who  found  also  a decrease  with  diminishing  particle  diameter  and  frequency  and  increasing  mass  rate  and 
velocity  difference  between  particle  and  gas.  The  conveying  capacity  is  limited  by  the  critical  velocity, 
which  is  equal  to  the  local  equilibrium  sound  velocity  and,  therefore,  the  mass-flow  rate  of  solid  particles 
cannot  be  increased  beyond  a maximum  value. 

2.6  Formation  of  Shock  Waves 


The  influence  contributed  by  the  particles  contained  in  the  gas  flow  is  qualitatively  shown  in  Fig.  5 
for  different  conditions.  If  we  introduce  the  relaxation  lengths  X,  for  velocity,  X2  for  temperature  and 
a characteristic  length  X , the  relative  value  of  these  lengths  is  determining  the  formation  of  a shock 
wava,  a fully  dispersed  wove,  or  a combination  of  a frozen  shock  wave  with  a dispersed  relaxation  zone. 

The  relaxation  lengths  and  X2  are  given  by 


Xl,2  = Vl,2 


where  ao  is  the  sound  velocity  at  a reference  state  of  the  gas,  and  T ^ 2 are  given  in  2.2.  Let  u be 
the  gas°velocity  at  a large  distance  before  and  u_,  u^  the  equilibrium’ and  the  frozen  velocity  beRind 
the  shock  wave;  c is  the  equilibrium  sound  velocity  and  b the  frozen  sound  velocity  of  the  mixture  of  gas 
and  particles  at  ¥arge  distance  before  the  shock.  The  results  show  that  a gas  with  solid  particles  behaves 
like  a pure  gas  with  two  inherent  relaxation  processes  (see  also  E.6.1). 


3.  APPLIED  RESEARCH 


3.1  Generation  of  Solid  Particles 


For  laboratory  investigations,  the  generation  of  a dusty  flow  is  of  interest.  Fuchs  and  Murashkevich 
describe  a laboratory  powder  disperser  (21).  The  efficiency  of  such  a device  may  be  checked  by  short-time 
micro-holography  (22). 

3.2  Conveying  of  Solid  Particles 

There  is  no  activity  in  Germany  on  two-phase  effects  in  rocket  propulsion.  However,  internal  flows 
with  solid  particles  are  widely  investigated,  and  survev  articles  on  pneumatic  conveying  in  horizontal 
straight  or  curved  tubes  are  given  by  Brauer  and  Schmldt-Traub  (23)  and  Vollheim  (24,25).  The  similarity 
laws,  valid  for  these  phenomena  were  checked  by  Siegel  (26). 

Experiments  of  Krotzsch  (27)  on  the  gas  flow  with  solid  particles  in  vertical  tubes  brought  the 
interesting  result  that  below  a limiting  concentration,  the  particle  velocity  is  determined  only  by 
collisions  with  the  wall,  whereas  collisions  between  particles  are  negligible.  The  behavior  is  very  similar 
to  that  of  rarefied  gas  flow,  and  a mean  free  path  d/6<x/2  should  be  introduced,  where  d is  the  mean  particle 
diameter  ano  c the  volume  concentration  of  the  solid  phase.  The  limiting  concentration  of  the  solid  phase 
then  is  reached  when  the  mean  free  path  is  equal  to  the  tube  diameter  D;  this  yields 

_ _d 

6R/2 

(see  also  C.6).  Another  aspect  of  the  transport  of  solid  particles  is  the  structure  of  turbulence,  which 
was  Investigate'  by  Kazanskij  (28)  for  water  flow  through  pipes  and  the  conveying  noise,  which  was  investi- 
gated by  Flatov  (29).  Noise  increases  with  increasing  gas  velocity  and  with  increasing  mass  rate  of  solid 
particles;  the  noise  level  is  lower  with  smaller  diameter  of  the  particles. 

3.3  Separation  of  Solid  Particles 

Although  industrial  dust  separators  are  outside  the  scope  of  this  review,  some  results  on  '•yclone 
separators  are  also  of  aerodynamic  Interest.  Since  the  early  work  of  Prandtl  on  cyclone  separators,  their 
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flow  behavior  has  often  been  Investigated.  For  the  cleaning  of  chimney  gas,  the  cyclones  are  now  replaced 
by  electrostatic  or  other  filters.  However,  they  are  further  used  In  chemical  engineering  for  separation 
of  high  particle  loads  at  relatively  low  Reynolds  numbers.  Separation  takes  place  In  the  Inlet  bend  and  In 
the  swirl  flow  of  the  Immersion  tube.  Calculation  of  optimal  cyclones  leads  to  extremely  long  tubes  If  the 
wall  shear  stress  Is  not  taken  Into  account  (30)  and  according  to  Gloger  and  Lukas  (31)  the  efficiency 
may,  therefore,  be  Improved  by  applying  boundary-layer  suction. 

3.4  Measurement  Techniques  for  Gas  Flows  With  Solid  Particles 

3.4.1  Analysis  of  Particle  Size 

A survey  on  methods  for  the  measurement  of  particle  sizes  is  given  b;  Leschonski  (32).  In  the  sedimen- 
tation analysis,  the  settling  velocity  in  a suitable  liquid  is  compared  with  that  of  an  equivalent  sphere. 
With  particle  sizes  below  1 pm,  Brownian  motion  influences  the  motion  in  the  gravitation  field  and,  there- 
fore, sedimentation  in  a centrifugal  field  is  preferred.  The  precision  of  the  sedimentation  /jnalysls  in 
a centrifugal  field  was  investigated  by  Alex  (33).  Other  methods  are  sifting,  optical  measurement,  and 
specific  surface  measurement.  The  application  of  Laser  light  sources  allows  not  only  to  measure  the  size 
of  particles  but  also  their  concentration  (34). 

3.4.2  Measurement  of  Mass  Rate  of  Solid  Particles  in  a Gas  Flow 


For  the  measurement  of  mass  flow  rate  of  solid  particles  in  a gas  flow  numerous  methods  are  used,  which 
range  from  photoelectric  dust-density  meters  to  instruments  which  measure  absorption  of  radioactive  radia- 
tion (35)  or  scattering  of  light  (36).  In  nearly  all  cases,  the  measuring  system  Is  in  a bypass  of  the  main 
flow  and  according  to  Fig.  6,  the  particle  concentration  is  measured  too  high  if  the  suction  velocity 
In  the  probe  Is  lower  than  the  freestream  velocity  W0  or  too  low  In  the  opposite  case.  The  probe  diameter 
and  the  diameter  of  the  limiting  stream  tubes  of  the  gas  D and  of  the  particles  Dg  are  indicated  in 
the  figure.  The  possible  errors  are  calculated  by  Bohnet  (37).  Zenker  (38)  has  verified  the  calculations 
by  comprehensive  measurements  (see  also  C.8.2). 

3.5  Sampling  of  Mlcrometeorltes 

In  some  satellite  or  high-altitude  rocket  experiments,  cosmic  dust  was  collected  and  analyzed. 

Extremely  small  particles  do  not  hit  the  collecting  surface  but  follow  the  streamline,  whereas  larger 
particles  are  less  deviated  due  to  inertial  forces.  An  analysis  was  given  by  Wuest  (1)  for  some  two-dimen- 
sional flow  configurations  such  as  supersonic  wedge  flow  and  viscous  hypersonic  stagnation  flow  around  a 
sampling  plate.  A representative  result  is  given  in  Fig.  7,  where  Tw  is  the  wall  temperature  and  Tf 
the  recovery  temperature  of  the  heat-insulated  plate,  and  pp  are  the  undisturbed  densities  of  the  gas  and 
of  the  particles  at  a large  distance  from  the  plate,  d is  the  particle  diameter,  and  c«  the  particle  drag 
coefficient.  The  ratio  of  the  width  of  the  particle  stream  tube  which  impacts  on  the  plate  is  plotted  as 
function  of  the  parameter  a defined  in  the  figure.  These  results  show  that  the  deviation  function  is 
greatly  Influenced  by  the  temperature  of  the  plate  if  the  flow  is  free-molecular  or  near  free-molecular. 


4.  CONCLUSIONS 

Main  Fundamental  Problems  to  be  Solved  in  the  Future 

1.  Drag  force  and  heat  transfer  on  spheres  are  known  in  different  regimes.  Additional  measurements 
are  needed  in  transitional  low  density  regime  for  varying  surface  temperatures. 

2.  In  the  low-Reynolds  number  regime  (near  free  molecular  flow)  the  heat  transfer  is  dependent  on 
surface  temperature.  No  experimental  data  are  available  on  this  behavior. 

3.  Only  scarce  measurements  on  nonspherical  particles  are  available. 

4.  The  behavior  of  small  particles  in  turbulent  flow  may  depend  on  the  ratio  of  particle  diameter  and 
turbulence  length  scale.  No  measured  data  on  these  phenomena  are  available. 


5.  Gas  flow  with  solid  particles  may  be  largely  influenced  by  Che  ratio  of  particle  mean  free  path  to 
a characteristic  length.  Only  few  measurements  are  known  where  this  influence  has  been  recognized.  Funda- 
mental experiments  are  needed  to  get  a clear  insight.  Theoretical  approximations  should  follow  the  lines 
of  the  BKG-model  in  kinetic  theory  (39).  This  may  give  an  adequate  description  of  interaction  problems  in 
dense  particle  flow. 

6.  An  experimental  investigation  of  relaxation  phenomena  is  needed.  This  may  comprise  sound  and  shock 
propagation  measurement  or  a direct  measurement  of  the  velocity  and  temperature  lag  of  the  particles. 

7.  In  high-altitude  micrometeorite  problems,  the  particles  flow  in  an  ionized  gas.  In  view  of  such 
problems,  particle  flow  should  also  be  studied  in  an  ionized  gas. 
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Figure  4 Radiant  interchange  among  suspended  particles,  according  to 
W.  Schneider  (10). 
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Figure  7 Sampling  of  micrometeorites.  Deviation  function  for  hypersonic 
free  molecular  flow  around  a plate,  according  to  W.  Wuest  (1). 
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SUMMARY 

This  paper  represents  a survy  of  research  currently  being  carried  out  In  the  United  Kingdom.  Under- 
standing of  solid-gaseous  flows  is  much  less  complete  than  in  the  case  of  fluid  flows  partly  because  of  the 
limitations  Imposed  by  the  available  measuring  techniques.  Several  such  techniques  are  disucssed.  The 
complex  nature  of  gas-particle  flows  has  led  to  an  attempt  to  obtain  a generalized  empirical  solution  which 
Is  briefly  described.  Also  included  are  discussions  on  particle  deposition,  entrainment,  and  erosion  caused 
by  the  impact  of  airborne  particles. 
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LIST  OF  1YMB0LS 

Semidepth  of  flat  plate  or  radius  of  cylinder 

Cross-sectional  area  of  particle 

Drag  coefficient 

Particle  diameter 

Pipe  diameter 

Electrical  permittivity  of  particle  material 

Gravitational  acceleration 

Electric  current 

Leakage  current 

Wall  current 

Roughness 

Pipe  length 

Acceleration  length 

Mass  of  particle 

Mass  flow  rate 

Mass  flow  ratio,  (Mg/Mj) 

Pressure  drop  over  length  1 of  pipe 
Pressure  drop  due  to  acceleration 
Electrical  resistance  per  unit  length  of  pipe 
Reynolds  number 
Time 

Velocity  component  In  x direction 
Fluid  velocity  in  undisturbed  region 
Friction  velocity 

Velocity  component  In  the  y direction 
Nominal  fluid  velocity,  (M^/O ^nD^/4) 

Component  of  deposition  velocity  due  to  interparticle  electrical  repulsion 

Component  of  deposition  velocity  due  to  turbulence 

Settling  velocity 

Minimum  transport  velocity 

Pracedins  mg*  yank 
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Suffixes 


Particle  coordinates 

Distance  of  grazing  streamline  from  axis  of  undisturbed  flow 
Time  delay 

Coefficient  of  restitution 

Component  of  deposition  velocity  due  to  interaction  of  turbulence  and  electrical  components 
Coll,sion  efficiency 
Pipe  inclination 

2 

Friction  factor,  O^/pjV  /8) 

Dynamic  viscosity 
Kinematic  viscosity 
Density 

Electrical  charge  density 
Density  ratio,  (ps/of) 

Stokes  number  (see  Eq.  6) 

Wall  shear  stress 

Electrical  potential,  Function  of  ... 

Fluid 

Solid 

Suspension 


1.  INTRODUCTION 

While  a considerable  volume  of  research  work  is  currently  being  carried  out  in  the  United  Kingdom 
in  the  field  of  solid-gaseous  suspension  flows,  by  far  the  major  effort  is  related  to  applications  in  the 
chemical,  food  processing  and  the  materials  handling  industries.  Hjwever,  a number  of  these  researches  may 
be  of  interest  to  workers  in  the  aeronautical  industry,  and  in  thi;  paper  the  author  has  selected  a number 
of  researches  which  could  have  relevance  in  the  aeronautical  fielj. 

Now  although  the  paper  is  primarily  concerned  with  solid-gr.seous  flows,  the  author  also  discusses 
solid-liquid  suspension  flows  where  there  appears  to  be  a simi’arity  between  gaseous  and  liquid  flows  and 
where  the  use  of  a liquid  in  lieu  of  a gas  may  be  helpful  in  mxiel  studies. 

It  will  be  appreciated  that  our  understanding  of  the  mechanism  of  solid-gaseous  flows  is  much  less 
complete  than  in  the  case  of  fluid  flows,  and  this  stems  partly  from  the  limitations  imposed  by  the  available 
measuring  techniques.  However,  the  relatively  new  technique  of  laser-Doppler  anemometry  appears  to  offer 
the  prospect  of  measuring  both  the  particle  and  fluid  velocity  distributions  and  the  particle  concentration 
distribution,  and  such  information  is  likely  to  lead  to  a greater  physical  understanding  of  suspension  flows. 
This  technique  and  its  application  in  this  field  will  also  be  discussed  in  the  paper.  Where  bulk  flow 
behavior  is  to  be  examined,  the  technique  of  cross  correlation  has  been  used  very  successfully  following 
recent  research  and  development  work  carried  out  in  the  U.K. 

While  many  investigators  have  formulated  mathematical  models  for  the  turbulent  fLow  of  a suspension 
of  solids  under  conditions  of  dilute  concentration  and  including  compressiblity  and  heat  transfer  effects, 
the  complex  nature  of  such  flows  has  led  the  writer  to  attempt  to  obtain  a generalized  empirical  solution 
using  the  technique  of  dimensional  analysis,  for  the  friction  factor  and  other  flow  parameters. 

The  problems  of  particle  deposition  and  entrainment  and  the  erosion  caused  by  the  impact  of  airborne 
particles  such  as  sand  would  appear  to  be  of  some  interest  to  aeronautical  engineers,  and  these  too  are 
discussed. 

Inevitably,  in  a review  paper  of  this  nature,  the  writer  has  to  be  selective  and  it  is  hoped  that  the 
selection  of  items  for  disucssion  will  prove  to  be  of  interest  to  those  working  in  the  field  of  aeronautics. 

2.  PARTICLE  DEPOSITION 


2.1  Collision  Prediction 


In  the  design  of  impingement  and  inertia  filters  for  the  removal  of  particles  from  an  air  stream,  it 
is  necessary  to  be  able  to  predict  the  trajectories  of  the  particles  in  the  region  of  an  obstruction  placed 
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in  the  path  of  the  suspension  flow.  On  the  other  hand,  there  are  a number  of  practically  important  flow 
situations  in  which  it  is  undesirable  for  the  particles  to  collide  with  an  obstruction.  For  example,  in 
turbomachinery  applications  the  collision  of  airborne  particles  could  lead  to  the  deposition  of  particles 
on  the  blade  surface  with  the  possibility  of  impairing  the  aerodynamic  performance,  or  the  collision,  in 
the  case  of  sand  or  other  abrasive  particles,  could  lead  to  serious  erosion  (see  also  E.9). 


In  an  attempt  to  be  able  to  predict  the  conditions  for  collision,  Morsi  and  Alexander  (1-4)  have 
studied  the  trajectories  of  fine  airborne  particles,  1 < d < 100  pm,  in  a uniform  flow  field  around  various 
body  shapes  placed  in  a suspension  flow.  They  determined  the  trajectories  analytically,  using  the  velocity 
field  for  a potential  flow  around  various  body  shapes  for  the  description  of  the  fluid  velocities,  and 
combining  these  velocities  with  the  equations  of  motion  for  the  particles  given  by 

m(dug/dt)  « CpPfCUf  - us)2  A/2 

(1) 

m(dvs/dc)  “ - vg)2  A/2 

(2) 

where 

u = dx/dt 
s 

(3) 

Vfi  ■ dy/dt 

(4) 

and 

Cp  = Kj/Re  + K2/Re2  + Kj  (5) 

In  Eq.  (5),  the  Reynolds  number  Is  defined  in  terms  of  the  x-component  of  the  relative  velocity,  thus 
Re  ■ (uf  - ug)d/v,  or  in  terms  of  the  y-component  of  the  relative  velocity,  Re  « (vj  - vg)d/v.  The 
constants  K.,  K.  , Kj  are  selected  to  give  a good  fit  over  the  range  of  Re.  This  analysis  represents  an 
approximation,  because  it  is  based  on  different  Reynolds  numbers  for  the  x-  and  y-directions . It  should 
not  introduce  large  errors  as  long  as  the  component  of  the  relative-velocity  vector  is  much  smaller  in  one 
coordinate  direction  than  in  the  other.  This  assumption  apparently  was  implied  by  Alexander  and  Moral  but 
not  explicitly  stated.  It  decouples  Eqs.  (1)  and  (2)  and  thus  makes  them  more  tractable.  A rigorous 
treatment  would  have  to  compute  the  drag  force  in  the  direction  of  the  relative-velocity  vector  and  then 
use  the  x-  and  y-components  of  this  force  in  Eqs.  (1)  and  (2)  (see  also  E.3.1).  Numerical  solution  of  the 
equations  made  it  possible  to  determine  the  grazing  trajectories  and  hence  the  collision  efficiency  Tl  ; 
the  latter  is  defined  by  y0/a  for  a flat  plate  or  a circulat  cylinder,  where  a is  the  semi-depth  of  tfie 
flat  plate  or  the  radius  of  the  cylinder,  and  y is  the  perpendicular  distance  from  the  grazing  trajectory 
in  the  undisturbed  region  to  the  axis  of  symmetry.  Typical  trajectories  for  a flat  plate  and  for  a cylinder 
are  shown  in  Figs.  1 and  2 for  different  Stokes  numbers,  where  the  Stokes  number  is  defined  by 

o = (2/9)  (Ua/v)  (p*)  (d/a)2  (6) 


with 


p*  - p ahf 

It  should  be  noted  that  in  the  case  of  the  flow  over  an  inclined  flat  plate  or  an  aerofoil,  the  distortion 
of  the  trajectories  and  the  asymmetry  of  the  flow  means  that  the  collision  efficiency  has  to  be  defined  in 
terms  of  the  ratio  of  the  vertical  area  of  the  powder  flow  which  collides  with  the  flat  plate  or  aerofoil 
and  the  vertical  projected  area  of  the  flat  plate  or  aerofoil. 

In  the  work  described,  it  will  be  apparent  that  the  gravitational  force  on  the  particle  has  been 
neglected,  and  while  no  great  error  is  Introduced  in  the  case  of  particles  having  low  settling  velocities, 
this  is  unlikely  to  be  so  for  larger  or  more  dense  particles,  and  a term  (mg)  should  be  Included  in  Eq.  (2). 
Alexander  and  Moral  have  in  fact  included  the  latter  and  have  shown  that  gravitational  effects  cause  a 
vertical  displacement  of  the  trajectories. 

The  work  would  appear  to  provide  a useful  basis  for  further  studies  into  the  deposition  of  particles 
on  aerofoils  and  other  body  shapes  and  for  erosion  and  Impact  studies  (see  also  A. 4,  B.3.5,  E.3  and  9). 

2.2  Deposition  Velocity  of  Charged  Particles 

Work  in  the  field  of  Electro-Gas  Dynamics  (EGD)  is  currently  being  studied  by  Musgrove  (5,6),  and  one 
of  the  interesting  problems  examined  is  concerned  with  the  drift  velocity  of  the  charged  particles  towards 
the  walls  of  the  EGD  generator  duct.  This  study  arose  because  it  was  thought  that  the  low  efficiency 
obtained  when  an  EGD  generator  was  operated  at  atmospheric  pressure  was  largely  due  to  parasitic  fields 
caused  by  the  deposition  of  charge  on  the  walls  of  the  duct. 
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In  this  study,  silicon  carbide  particles  having  a mean  diameter  of  4.5  were  transported  through  a 
hlgh-soda^glaas  duct  of  50  mm  In  diameter  (Pig.  3a)  at  a velocity  of  54  m/s;  the  particle  concentration  was 
0.30  kg/m  , corresponding  to  a number  density  of  1012/m3.  The  resistivity  was  maintained  constant  by 
temperature  control  at  a value  of  304  GQ/m.  Under  these  conditions,  Husgrove  showed  that  the  deposition 
velocity,  V , could  be  determined  from  measurements  of  the  electrical  potential  distribution  along  the  duct; 
he  obtalnedua  value  of  0.22  m/s. 

In  his  analysis,  he  assumed  that  the  deposition  velocity  Is  made  up  of  three  components:  a component 

arising  from  the  fluid  turbulence,  Vt,  a component  due  to  Inter-particle  electrical  repulsion,  Ve,  and  a com- 
ponent T]VtVe  due  to  the  Interaction  between  the  turbulence  and  electrical  components.  The  value  of  T) 
must  be  determined  experimentally.  He  further  assumed  that  the  radial  field  Is  proportional  to  the  charge 
density,  and  thus  Vfi  “ (i/i  )V  , where  the  suffix  (o)  refers  to  conditions  at  x - 0,  at  the  upstream  end 
of  the  conducting  duct  section.  Now  it  will  be  clear  from  Fig.  3b  that  the  sum  of  the  current  1 through 
the  duct  due  to  the  flow  of  the  charged  particles  and  the  wall  current  due  to  the  transfer  of  the  charge 

to  the  duct  wall  by  the  deposited  particles  will  be  equal  to  the  leakage  current  1.,  where 

L 


-(di/dx)dx 


2TtrP  v dx 
e D 


' "r2peVf 


-di/dx  = 2(i/r)  (VD/Pf) 


Now  from  the  above  assumptions  we  have 


VD  - vtfi  +*/i0>Veo(l/Vt  +TI>] 


and  substituting  for  Vp  from  Eq.  (9),  we  obtain 


-(1/1)  (di/dx)  - (2/4:)  (Vt/V£)  fl  + (i/i0>Veo(1/Vt  +T1)] 


In  this  equation,  the  current  1 may  be  determined  from  the  potential  gradient  and  the  current  1.  using 
the  equation 

1 - 1£  - (1/R)  (dp/ dx)  (12 


where  d<p/dx  is  obtained  by  numerical  or  graphical  differentiation.  Musgrove  thus  solved  Eq.  (11)  by  plotting 
- (1/1)  dl/dx)  against  i/i0,  obtaining  2Vt/rV£  from  the  intercept  corresponding  to  1/ i0  » 0,  and  2VtVeo(l/Vt+ 
Tl)/rVf,  from  the  slope  of  the  straight  line.  If  T|  is  assumed  to  be  zero  then  Vt  and  Veo  are  determined 
directly,  and  hence  Vp.  However,  if  the  experiment  Is  repeated  at  a different  fluid  velocity  while  main- 
taining the  current  1Q  and  hence  Veo  constant  at  the  same  value  for  each  velocity,  then  Tt  may  be  determined 
and  the  assumption  Tl  ■ 0 Is  unnecessary. 

In  cases  where  the  current  flow  is  small,  the  turbulent  component  of  the  deposition  velocity  dominates 
and  In  this  case  Musgrove  showed  that 


<p  - (RIq/KL)  [(1  - e_Kx)/x  - (1  - e'KL)/L] 


where  K ■ 2Vp/rVf  and  L Is  the  length  of  the  conducting  tube.  Clearly,  a single  measurement  of  che  potential 
<p,  corresponding  to  a given  value  of  x suffices  to  determine  Vp.  He  also  showed  that  when  Vp  « rVf , the 
potential  distribution  Is  parabolic,  and  that  In  this  case, 

V_  - 4q>  V./l  RL2  (14) 

D max  f o ' 


The  determination  of  Vp  using  Eqs.  (13)  or  (14)  is  preferable  to  that  based  on  Eq.  (11),  since  the 
latter  requires  that  the  electrical  potential  should  be  differentiated  twice  and  this  Introduces  serious 
computational  errors. 

While  this  technique  of  determining  the  deposition  velocity  was  developed  In  respect  of  EGD,  It  is 
likely  to  find  application  In  a number  of  other  fields  and  could  be  used  for  a wide  range  of  experimental 
conditions. 
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3.  PARTICLE  REMOVAL 

Having  discussed  at  some  length  the  work  concerned  with  the  deposition  of  particles,  we  now  turn  our 
attention  to  some  recent  work  carried  out  by  Cleaver  and  Yates  at  Liverpool  University  (7)  on  the  problem 
of  particle  detachment  of  submicron-sized  particles.  In  this  work,  they  develop  a criterion  for  the  removal 
of  small  particles  which  is  based  on  the  non-steady  behavior  of  the  viscous  sub-layer  of  a turbulent 
boundary  layer.  They  suggest  that  the  turbulent  bursts  which  continually  disrupt  the  viscous  layer,  which 
are  not  unlike  miniature  tornadoes,  may  be  sufficient  to  cause  instantaneous  lift  forces  sufficient  to 
detach  a particle. 

In  arriving  at  their  criterion  for  particle  removal,  they  assume  that  the  force  responsible  for  the 
removal  is  the  hydrodynamic  lift  force  and  use  the  equation  due  to  Saffman  (8),  which  shows  that  the  lift 
is  proportional  to  p fv2(dU*/v)  , where  U*  is  the  friction  velocity.  Since  previous  work  has  shown  that 
the  adhesive  force  is  proportional  to  the  particle  diameter,  and  since  the  particles  will  only  be  removed 
if  the  ratio  of  the  lift  force  to  the  adhesive  force  is  greater  than  unity,  they  arrive  at  the  relationship, 

(Pfv2/<>)  (dU*/v)3  > Const.  (15) 

Since  U*  “ (tq/v)^,  then  for  a given  fluid  the  criterion  for  particle  removal  reduces  to 

T d4/3  > B (16) 

o 


where  T is  the  wall  shear  stress,  d is  the  particle  diameter,  and  B is  a constant.  It  is  this  criterion 
that  Cleaver  proposes  for  the  removal  of  particles,  and  although  he  does  not  give  values  of  B for  particles 
of  different  shape  or  for  different  types  of  adhesive  force,  he  uses  the  experimental  dat^  of  Vlsser  (9) 
for  carbon  black  on  a cellophane  substrate  to  show  that  T0d^' J has  the  same  order  of  magnitude  for  different- 
sized particles,  a value  of  337  x 10  6 g-cm1' ^le^^tor  particles  of  0.21  i*m  in  diameter.  The  consistent 
units  to  be  used  in  Eq.  (16)  are  therefore  g.cm/s  for  Tq  and  cm  for  d. 

To  assess  the  rate  at  which  particles  are  removed,  he  assumes  that  only  one  burst  of  turbulence 
occurs  at  any  one  time  over  an  area  equal  to  630(v/U*)  x 135(v/U*),  and  that  the  bursts  repeat  over  a time 
interval  of  75  V/U*  , the  location  of  the  bursts  being  random.  He  also  assumes  that  the  area  of  the 
burst  is  equal  to  (rr/4)(20  v/U*)2,  and  that  a proportion  a of  this  area  is  cleaned.  By  dividing  the  control 
area  into  a number  of  sites,  m,  where 


m - (630  v/u*)  x (135  v/U*)/0.7854i20  V/U*)2 


(17) 


he  points  out  that  after  n bursts,  where 


n - t(U*2/75v) 


(18) 


the  probable  percentage  of  cleaned  surface  is  given  by 


S(t)  = 1 - (1  - or/270)  (U*2/7 5v)  t 


(19) 


where  a is  assumed  to  be  of  the  order  of  1/100. 


4.  EROSION  STUDIES 

During  the  transport  of  abrasive  materials,  such  as  sand  and  carborundum  in  shot  blasting  operations, 
wear  at  pipe  bends  and  in  the  pipelines  presents  a serious  problem.  Arundel  (10)  and  his  co-workers  at 
Liverpool  Polytechnic  have  been  carrying  out  experimental  studies  to  assess  the  relative  importance  of 
solids  loading,  air  velocity  and  the  angle  of  impact.  The  results  of  this  work  indicate  that,  foi  veloci- 
ties less  than  100  m/s,  Llnatex  rubber  appears  to  have  a high  resistance  to  wear,  though  at  higher  veloci- 
ties mild  steel  and  tool  steel  are  significantly  better.  The  wear  rate  for  mild  steel  is  given  by  the 
equation 


q * 40  x 10  *(M*)^‘3  v3'3  mm3 /minute 


(20) 


where  the  impact  angle  is  in  the  range  from  20°  to  30°,  the  mass  flow  ratio  M*  lies  between  0 and  8,  and 
the  air  velocity  V is  155  m/s. 

In  their  future  work  they  intend  to  study  the  wear  rate  in  bends  in  pipes  having  diameters  in  the 
range  20  < D < 40  mn  where  the  bend  radius  will  be  varied  up  to  a maximum  of  6D.  The  particle  velocities, 
trajectories  and  concentration  in  the  bend  will  also  be  studied  by  laser-Doppler  anemometry.  Ultrasonic 
techniques  will  be  used  to  locate  the  regions  of  maximum  wear  (see  also  A. 4,  B.3.5,  and  E.9). 
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INGRESS  OF  FOREIGN  BODIES  INTO  AIRCRAFT  ENGINE  INTAKES 


The  problem  of  the  Ingress  of  foreign  bodies,  such  as  birds,  hailstones,  Ice,  etc.  In  the  atmosphere, 
or  of  stones,  grass,  sand,  etc.  on  the  ground,  or  manufacturing  debris  left  In  the  aircraft  following 
manufacture  or  maintenance.  Is  a problem  that  Is  being  given  serious  consideration  In  the  U.K.  (11). 

A number  of  actions  may  be  taken  to  avoid  the  problem,  by  avoiding  areas  of  high  bird  concentration, 
good  housekeeping  In  respect  of  manufacture  and  maintenance,  and  by  avoiding  ln-llne  ahead  take  off. 

Although  such  actions  are  helpful,  further  research  and  development  Is  required.  Studies  of  the  flow 
conditions  required  for  the  pick-up  of  particles  and  of  the  trajectories  of  particles  in  the  vicinity  of  the 
aircraft  in  general  and  the  engine  intakes  in  particular  should  be  undertaken  In  order  that  the  conditions 
for  the  ingress  of  foreign  bodies  may  be  avoided  or  that  preventive  action  may  be  taken  by  careful  design 
and  location  of  the  engine  Intakes.  A considerable  amount  of  work  has  already  been  carried  out  In  relation 
to  pick-up  velocities,  particularly  In  the  case  of  small  particles  of  sand.  This  was  studied  extensively 
some  two  decades  ago  by  Bagnold  (12).  More  recently  Bagnold  (13)  has  been  studying  the  nature  of  saltation 
of  particles  In  water,  and  this  too  could  be  helpful  even  though  the  flow  conditions  are  not  directly 
comparable.  During  the  commissioning  of  power  stations  the  pipework  has  to  be  purged  to  remove  debris 
Introduced  during  erection  which  could  be  carried  Into  the  turbines.  A recent  study  of  this  problem  by 
Cranfleld  and  Lawrence  (14)  gives  some  useful  Information  on  pick-up  velocities,  and  the  analysis  provides 
a basis  for  the  further  study  of  the  pick-up  of  large  particles  close  to  aero-engine  Intakes. 

Design  studies  should  be  undertaken  in  parallel  with  the  research  work,  to  examine  the  feasibility  of 
modifying  the  design  and  location  of  the  engine  Intakes,  and  of  providing  guards,  separators,  and  filters 
upstream  of  the  Intake  (see  also  A. 4,  B.3.5,  and  E.9). 

6.  EMPIRICAL  CORRELATIONS  FOR  SUSPENSION  FLOWS 

Though  many  workers  have  developed  empirical  correlations  for  the  various  flow  characteristics  of 
suspensions,  with  few  exceptions  these  workers  have  confined  their  studies  to  a limited  range  of  the  var- 
iables involved.  It  thus  appeared  to  the  author,  that  there  would  be  merit  in  attempting  to  generalize 
the  equations  for  suspension  flows,  using  for  this  purpose  the  technique  of  dimensional  analysis.  A 
detailed  account  of  this  treatment  is  given  by  Rose  and  Duckworth  (15)  who  deal  with  the  flow  of  suspensions 
of  particles  in  liquids  and  gases  in  pipelines  at  various  inclinations  from  the  horizontal  to  the  vertical 
position.  Equations  are  obtained  for  the  length  of  pipe  required  for  the  particles  to  accelerate  from 
rest  to  the  velocity  corresponding  to  the  conditions  of  established  flow,  for  the  pressure  drop  required 
to  accelerate  the  particles  to  this  velocity,  for  the  pressure  gradient  required  to  maintain  the  flow  of  as 
suspension  under  established  flow  conditions  and  for  the  minimum  fluid  velocity  required  for  the  transport 
of  the  particles  in  the  form  of  a suspension.  They  show  that,  to  an  acceptable  degree  of  accuracy,  the 
pressure  drop  for  established  flow  may  be  computed  from  the  pressure  drop  In  the  pipe  when  conveying  the 
fluid  phase  only  plus  a pressure  drop  which  depends  upon  the  characteristics  of  the  solid  phase.  Data 
relating  to  pipe  lines  varying  in  diameter  from  30  to  400  mm,  conveying  solids  of  specific  gravities  of 
1.5  to  11.1  In  liquids  and  gases,  are  satisfactorily  correlated  by  means  of  the  equations  developed.  These 
equations,  some  of  which  have  been  subsequently  modified  by  the  author,  are  summarized  below  (see  also 
B.3.2). 


dP, 


sA 


la/d 


/ (p  fV  /2) 


<{(Ms/pfg*D5/2)  x (D/d)*  x (p*)*J 


M*  x ^ (V2/gdp*2)  x a92(9) 


(21) 


(22) 


dpm(p fV2/2) (L/D)  - q£PfVD/p,k/D]  + [^(M*)  x ^(d/D)  x ^(e)  x ^(P*)  x F<P5(0) 
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Vg/Vf  - 0.5  x ^(V^gdP*2)  x A<P2 (9 ) 


Ap  -Apm  + APgA  + PfV  /2 


V /V 

C 00 


ytp^d/D)  x vT>2(0)  x M* 
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(24) 


(25) 


(26) 


The  function  Aqj, 


and  ip  In  these  equations  refer  to  correlations  for  acceleration,  friction  and 


velocity  and  can  be  found  in  Rif.  15.  As  an  example,  the  relationship  between  the  solids  friction  factor 


...  — -ound  In  R^f. — r , r 

function  jff>g(V  /gD)  and  the  Froude  number  VZ/gD  Is  shown  in  Fig.  4.  It  should  be  noted  that,  although 
Eq.  (26)  Is  of  the  same  form  as  that  used  In  Ref.  15,  the  equation  Is  based  on  subsequent  experimental  work 
carried  out  by  the  author  for  the  transport  of  particles  In  air.  The  equation  covers  the  transport  of 
uniformly  sized  particles  of  glass,  sand,  flyash,  sodium  bicarbonate,  aluminum  silicate,  polyester  and 
seeds.  In  the  size  range  0.006  < d/D  < 0.1  and  for  pipe  Inclinations  In  the  range  0 < 9 < 90°. 
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It  should  also  be  noted  that  the  above  equations  are  restricted  to  particles  of  sensibly  spherical 
shape  having  a narrow  size  distribution.  Current  work  Is  concerned  with  the  study  of  the  Influence  of  the 
particle  shape  factor,  the  spread  of  particle  sizes,  particle  cohesiveness,  and  the  Influence  of  electro- 
static effects.  In  a recent  series  of  experiments  (16)  by  the  author,  In  connection  with  the  Influence  of 
electrostatic  charge  level,  on  the  pressure  gradient,  it  was  found  that  the  solids  friction  factor  Xg 
showed  a significant  Increase  as  the  dimensionless  parameter  l2/eM  V3  was  Increased,  though  further  work  Is 
required  before  definitive  conclusions  can  be  made  in  respect  of  the  correlation. 


7.  APPLICATION  OF  LASER-DOPPLER  ANEMOMETRY  TO  SOLID-GASEOUS  FLOWS 

The  laser-Doppler  technique  Is  particularly  appropriate  for  the  measurement  of  two  component  flows,  in 
that  It  does  not  Interfere  with  the  flow  and  makes  the  measurement  of  local  velocities  and  concentrations 
possible  (see  also  A. 3,  D.5,  and  E.9).  Up  to  the  present  time  only  a limited  number  of  preliminary  Investi- 
gations have  beet,  carried  out  In  the  U-K. , using  single  and  two-mode  particle-size  distributions  flowing 
in  vertical  ducts  (17,18).  This  work  has  indicated  that  it  is  not  necessary  to  use  mono-sized  particles, 
since  a signal  is  obtained  from  each  particle.  It  is  intended  to  extend  the  present  work  to  include  the 
measurement  of  the  particle  sizes,  velocity  and  concentration  distribution.  Future  work  will  examine  the 
radial  distribution  of  the  axial  and  radial  velocities  of  the  particles  and  the  fluid  carrier,  and  the 
particle  and  fluid  Reynolds  stresses,  in  an  attempt  to  refine  the  models  for  the  prediction  of  the  shear 
stress  and  pressure  gradients  in  suspension  flows.  However,  the  technique  appears  to  be  limited  to  the  study 
of  relatively  dilute  concentrations  of  the  order  of  1.5%  by  volume  and  15:1  by  mass. 

At  Southampton  University  (17)  flow  fields  in  flames  as  well  as  in  cold  flows  have  been  seeded  with 
either  solid  particles  or  aerosols,  in  order  to  make  velocity  measurements.  Instruments  have  been  developed 
which  make  simultaneous  measurements  of  all  three  components  of  velocity,  in  real  time,  in  hot  and  cold 
turbulent  flows.  Intensities  , spectra,  scales  and  Reynolds  stresses  of  the  turbulence  have  been  deduced. 
Preliminary  studies  of  the  effects  of  particle  loading  on  the  turbulent  flow  have  been  made,  but  the  main 
interest  has  been  in  conditions  where  the  particles  do  not  significantly  influence  the  flow. 

Crossed-beam  correlation  studies  have  also  been  carried  out  at  Southampton,  in  which  measurements  of  the 
covariance  of  the  signals  from  two  crossed  laser  beams  in  a two-phase  turbulent  flow  have  yielded  statisti- 
cal information  about  the  turbulence  in  Che  vicinity  of  the  point  of  intersection  of  the  two  beams. 


8.  MEASURING  TECHNIQUES  FOR  VELOCITY,  VOLUME  FLOW  AND  MASS  FLOW 

8.1  Velocity  and  Volume  Flow  Measurement 

In  industrial  applications,  particularly  in  the  chemical  and  food  processing  Industries,  we  frequently 
require  to  monitor  the  mean  velocity  and  the  volume  flux  of  the  suspension,  and  it  is  clearly  desirable  to 
avoid  interference  between  the  instruments  used  and  the  suspension  flow.  For  this  purpose,  cross-correla- 
tion techniques  of  measurement  have  been  studied  extensively  in  the  U.K.  by  Beck  and  his  colleagues  at 
Bradford  University  (19). 

In  these  studies,  two  capacitance  transducers,  which  consist  of  an  electrode  fitting  flush  with  L.ie 
internal  surface  of  the  duct  and  insulated  from  it,  are  located  at  stations  A and  B separated  by  a distance 
L along  the  length  of  the  duct  (see  Fig.  5).  During  the  passage  of  the  suspension,  the  capacitance  varies 
and  produces  a signal  m(t)  at  the  upstream  transducer,  and  a signal  n(t)  at  the  downstream  transducer. 

The  time  delay  t between  the  random  signals  produced  by  the  transducers  is  related  to  the  volume  flux 
through  the  duct.  To  determine  it,  and  hence  the  volune  flux,  an  adjustable  time  delay  8 is  inserted  into 
the  output  from  the  upstream  transducer;  the  product  of  the  delayed  upstream  signal  and  the  downstream  signal, 
m(t  - B)  ■ n(t),  is  integrated  over  a time  to  yield  a cross-correlation  function  for  a range  of  values  of 
B,  as  shown  in  Fig.  5.  To  determine  T,  use  is  made  of  the  fact  that  the  cross-correlation  function  will 
be  a maximum  when  the  delay  of  the  upstream  signal  8 is  equal  to  the  time  delay  t of  the  random  signals 
over  the  distance  L between  the  two  transducers.  Once  T is  determined,  it  is  a simple  matter  to  determine 
either  the  suspension  velocity  or  the  volume  flux  since, 

Suspension  Velocity  “ L/t 

Volume  Flux  ■ (Volume  of  Pipe  Over  Length  L)(t 

This  technique  has  been  used  successfully  for  the  on-line  control  of  chemical  process  plants  and  for 
the  metering  of  dust-laden  gas  streams  (see  also  A. 5. 2). 

8.2  Mass  Flow  Measurement 

Since  the  mass  flux  of  solids  carried  in  suspension  in  a carrier  fluid  affects  the  turbulence  level  of 
the  flow,  and  since  the  turbulence  level  can  be  monitored  by  a capacitance  transducer,  it  follows  that  the 
mass  flow  of  the  solids  can  be  measured  by  means  of  A/C  voltmeter  connected  to  the  transducer  (20). 

Although  this  method  of  mass  flow  measurement  requires  calibration,  it  has  been  used  quite  satisfactorily 
for  the  flow  of  a wide  range  of  powders  and  granular  materials,  including  cement,  flour,  wheat  and  barley 
(see  also  B.3.4.2). 

The  techniques  can  be  used  for  the  measurement  of  nongravlmetric  quantities,  such  as  electrostatic 
charge,  and  it  has  been  suggested  by  the  Bradford  workers,  that  they  could  be  used  for  measuring  the  amount 
of  charge  being  carried  in  an  aircraft  refueling  pipe.  In  the  latter  application,  abnormally  high  charge 
conditions  exist,  and  these  can  lead  to  explosions  which  could  be  avoided  if  the  above  techniq""*  were  used 
for  detection. 
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9.  CONCLUSION 

In  Chls  paper  a number  of  researches  in  Che  field  of  multiple  flows  of  solids  in  gases  which  have 
relevance  to  aeronautical  applications  have  been  briefly  outlined,  though  with  few  exceptions  the  motiva- 
tion for  the  research  has  not  come  from  the  aircraft  Industry.  A number  of  other  researches  have  been 
discussed  which  are  primarily  of  fundamental  interest  and  where  the  techniques  used  are  likely  to  find 
general  applications. 

The  work  on  particle  deposition  and  collision  has  application  to  the  build-up  of  particulate  matter 
on  compressor  and  turbine  blading,  particularly  when  this  is  considered  in  conjunction  with  the  interesting 
work  on  aerodynamic  removal  by  Cleaver. 

While  the  prime  objective  of  aircraft  designers  must  be  to  avoid  the  ingress  of  foreign  bodies  such 
as  sand  into  the  engine,  and  the  work  of  Bagnold  and  others  should  be  helpful  in  this  respect,  nevertheless, 
it  is  useful  to  know  something  of  the  factors  effecting  erosion  due  to  particle  impact,  and  the  work  of 
Arundel,  et  al.  is  likely  to  be  helpful. 

Laser-Doppler  and  cross-correlation  techniques  are  widely  used  in  the  aeronautical  field,  and  it  will 
be  interesting  to  see  how  these  are  applied  to  the  field  of  two-component  flows.  The  suggestion  by  Beck 
that  a simple  capacitance  transducer  may  be  used  to  monitor  electrostatic  charge  levels  in  aircraft 
fuel  lines  is  of  particular  interest. 

In  conclusion,  it  is  hoped  that  this  brief  survey  of  work  in  the  U.K.  proves  to  be  useful  in  stimulating 
discussion  and  further  work  in  the  field  of  multiple  flows  of  solids  in  gases. 
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Figure  1 Particle  trajectories  in  the  vicinity 
of  a flat  plate  perpendicular  to  the 
direction  of  flow  for  two  values  of  the 
Stokes  number  (4).  (Dimensions  are  in 
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Figure  2 Particle  trajectories  in  the  vicinity 
of  an  inclined  flat  plate  (15°)  and  of  a 
circular  cylinder  for  two  values  of  the 
Stokes  number.  (Dimensions  are  in  mm) 
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Figure  5 Cross-conrelation  technique  used 
for  metering  of  suspension  flows  (19) 
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SUMMARY 

Research  done  at  the  VKI  is  briefly  reviewed.  The  topics  discussed  include  low-speed  and  high-speed 
gas-particle  flow.  Theoretical  and  experimental  aspects  are  discussed  be  well  as  some  work  on  Instrumen- 
tation. In  addition,  observations  in  a fluidized  bed  are  briefly  described,  and  development  of  a laser- 
Doppler  velocimeter  is  outlined. 
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LIST  OF  SYMBOLS 

particle  diameter 
static  pressure 
total  pressure 
particle  mass  flow 
gas  mass  flow 
particle  velocity 
gas  velocity 
longitudinal  position 
transversal  position 

pressure  differential  measured  in  a mixture 

pressure  differential  due  to  gas  only 

loading  ratio  • Q /Q 
P 8 

density  of  particle  phase  in  the  mixture 


1.  INTRODUCTION 

This  paper  is  a survey  of  the  basic  work  done  at  the  VKI  on  gas-particle  flows  ^nce  1970.  It  deals 
successively  with  low-speed  flows,  high-speed  flows,  fluidization  and  the  development  of  a laser-Doppler 
velocimeter. 


2.  LOW-SPEED  GAS-PARTICLE  FLOWS 

In  a flowing  suspension,  the  loading  ratio,  or  ratio  of  the  mass  flow  of  solids  to  that  of  the  gas, 
is  a very  important  parameter.  Buchlln  and  Rlethmulier  (1)  improved  an  existing  technique  of  instantaneous 
measurements  of  the  loading  ratio  in  a venturi  tube.  They  analyzed  the  theoretical  conditions  for  an 
optimum  geometry  and  developed  an  analog  data-processlng  system  which  allcws  a direct  reading  of  the  loading 
ratio.  It  was  experimentally  and  theoretically  demonstrated  that  the  ratio  of  parti  tie  and  gas  velocities, 
upstream  of  the  venturi  tube,  had  a strong  effect  on  its  calibration  factor.  Th.s  latter  can  be  reasonably 
well  predicted  as  demonstrated  by  the  good  agreement  found  between  experimental  and  theoretical  results. 
Typical  venturi-tube  calibration  curves  for  two  inlet  velocity  ratios  are  presented  in  Fig.  1,  where  curves 
a and  b refer  to  partlcle-to-gas  velocity  ratios  of  0.6  and  0.4.  Rlethmulier  and  Glnoux  (2)  demonstrated 
that  a detailed  analysis  of  gas-particle  flows  was  possible  by  using  a laser-Doppler  velocimeter  (see 
instrumentation  section).  Nouvellon  and  Rlethmulier  (3)  showed  experimentally  that  the  behavior  of  parti- 
cles entrained  by  a gas  flow  was  strongly  dependent  on  particle-wall  and  particle-particle  interactions. 

The  analysis  of  velocity  and  concentration  measurements,  performed  in  a horizontal  pipe,  led  to  the  con- 
clusion that  a bulk-flow  approach  was  not  sufficient  to  describe  a flowing  suspension.  A study  of  the 
kinetics  of  the  particles  was  thus  started  by  Bouton  and  Rlethmulier  (4).  Particle  trajectories  and 
velocity  distribution  functions  were  determined  with  a laser-Doppler  velocimeter  in  a suspension  flowing 
in  a complex  geometry.  Particle  velocity  and  concentration  profiles  measured  in  a horizontal  converging 
duct  are  represented  in  Fig.  2.  The  particle  diameter  for  these  experiments  is  d ■ 480  * 30  p.m,  and  the 
loading  ratio  is  71  * 1.98.  Axial  distance  x is  measured  from  an  arbitrary  origin  with  x ■ 100  cm  being 
located  at  the  inlet  of  the  test  section  (see  inset  of  upper  part  of  the  figure).  For  the  conditions  of 
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Fig.  2,  the  probability  density  function  of  the  particle  velocity,  measured  at  one  location  in  the  converg- 
ing duct,  is  plotted  in  Fig.  3.  From  the  analysis  of  these  results,  a physical  interpretation  of  the  flow 
of  a suspension  in  a complex  geometry  was  propos€d.  Particle-particle  interactions  were  shown  to  be  of 
small  importance  compared  to  the  wall-particle  interaction.  A computer  program,  based  on  a "semi-kinetic" 
model  of  the  flow  has  been  written.  The  results  of  the  calculation  were  in  good  agreement  with  the  experi- 
mental data. 


3.  HIGH  SPEED  GAS-PARTICLE  FLOWS 

Theoretical  and  experimental  investigations  of  gas-particle  nozzle  flows  have  been  performed  at  the 
von  Karrnan  Institute  since  1970.  Emphasis  was  given  to  the  effects  of  high  loading  ratios.  A one-dimen- 
sional analysis  has  been  programmed  for  an  IBM  1130  which  was  used  by  Dumortier  and  Ginoux  (5)  to  perform 
a dimensionless  parametric  study  of  gas-particle  nozzle  flows  in  which  velocity,  temperature,  and  pressure 
variations  were  computed.  Transonic  solutions  for  large  loading  ratios  have  also  been  examined.  They 
provided  the  initial  conditions  for  the  analysis  of  the  gas-particle  flow  In  an  axisymme trie  nozzle  which 
was  performed  by  Johnson  (6).  His  program  yielded  a numerical  solution  for  the  supersonic  part  of  the 
nozzle  in  which  the  partial  differential  equations  were  solved  by  the  method  oi  characteristics. 

Riethmuller  and  Thiry  (7)  measured  mass  flows  and  static  pressure  distributions  in  nozzles  for  high 
loading  ratios.  The  flow  was  fed  to  the  nozzle  from  a constant-pressure  reservoir,  and  the  diameter  of 
the  particles  was  250  ^m.  The  converging  part  of  the  nozzle  was  ellipsoidal  and  the  diverging  part 
conical  (1.5°  half-angle).  Typical  results  are  shown  in  Fig.  4,  where  the  mass  flow  of  the  gas  has  been 
plotted  against  the  mass  flow  of  the  particles.  Only  data  for  the  last  part  of  the  nozzle  are  shown, 
since  no  measurements  were  made  upstream  of  this  portion. 

The  analysis  was  based  on  simultaneous  integration  of  the  continuity,  momentum,  and  energy  equations 
for  both  the  gas  and  the  particles  (see  also  E.5  and  7).  Newtonian  drag  of  isolated  spheres  was  assumed. 

The  drag  coefficient  depends  on  both  Reynolds  number  and  Mach  number  but,  effectively,  varied  only  between 
0.4  and  0.6.  Each  particle  was  assumed  to  have  uniform  internal  temperature,  and  the  heat-transfer  equation 
was  written  for  isolated  spheres  with  forced  convection  only.  The  calculations  were  started  at  the  inlet 
section  of  the  nozzle  for  a prescribed  initial  particle  velocity,  but  it  was  found  that  the  results  are 
quite  insensitive  to  this  assumption.  A comparison  of  the  measured  static  pressure  distributions  with 
theoretical  results  is  presented  in  Fig.  5.  The  predicted  effect  of  the  particles  on  the  pressure  distri- 
bution is  Oood  agreement  with  the  experimental  results. 

After  a facility,  especially  designed  for  high-speed  gas-particle  flows  with  high  loading  ratios,  was 
completed,  an  investigation  of  a gas-particle  flow  over  an  obstacle  wac  performed  by  Bouffinier  and  Rieth- 
rauller  (8).  Impact  forces  of  a gas-particle  jet  on  a target  as  well  as  localized  impacts  on  a small  probe 
were  determined.  These  experimental  results  have  been  helpful  in  providing  an  insight  into  the  behavior  of 
a high-speed  gas-particle  jet  over  an  obstacle. 


4.  FLUIDIZATION 

An  investigation  of  the  hydrodynamic  behavior  of  a two-d imens ional  bed  in  the  presence  of  a cyl indrica I 
obstacle  was  performed  by  Ginoux,  De  Geyter,  and  Kilkis  (9). 

At  superficial  velocities  below  minimum  fluidization  conditions,  local  fluidization  was  photograph- 
ically observed  near  the  obstacle.  Information  on  the  dynamics  of  the  bubble  were  deduced  from  the  pressure 
variations  measured  on  the  surface  of  the  cylinder. 

At  minimum  fluidization  velocity,  a study  was  made  of  the  interaction  of  a single  bubble  with  the 
obstacle.  These  results  lead  to  a better  understanding  of  the  heat- transfer  mechanism  in  a fluidized  bed. 


5.  LASER-DOPPLER  VELOCIMETER 

The  development  of  a laser-Doppler  velocimeter  was  started  at  the  von  Karrnan  Institute  in  1972. 

Barker,  Riethmuller  and  Ginoux  (10)  set  up  a "fringe"  system  for  the  measurement  of  the  velocity  of 
large  solid  particles  (up  to  500  nm)  in  a two-phase  flow.  The  analysis  of  the  special  requirements  for 
large-particle  measurements  was  further  improved  by  Riethmuller  (11)  who  proposed  a new  model  of  the  quality- 
factor  distribution  curve  for  large  particles.  The  Doppler  signal  obtained  from  the  photodetector  is 
characterized  by  its  quality  factor  which  is  the  ratio  of  the  amplitude  of  the  Doppler  frequency  and  the 
overall  amplitude  of  the  signal.  This  quality  factor  greatly  depends  on  the  ratio  between  particle  dia- 
meter and  fringe  spacing.  It  was  generally  accepted  that  it  decreased  rapidly  for  large  values  of  the 
ratio,  but  a certain  number  of  experiments  led  to  the  conclusion  that  for  very  large  ratios  of  particle 
diameter- to -fringe  spacing,  the  quality  factor  should  rise  again  and  then  decrease  further.  This  departure 
from  the  previously  accepted  curves  should  be  dependent  on  the  ratio  of  the  particle  diameter  to  the  laser 
light  wavelength.  An  anticipated  quality-factor  distribution  curve  was  thus  proposed,  which  explained  why 
it  was  possible  to  measure  the  velocity  of  large  particles  with  a partic  e diameter  to  fringe  spacing  ratio 
larger  than  30.  Durst  and  Eliasson  (12)  have  confirmed  this  hypothesis. 

In  the  meantime,  an  electronic  Doppler  processor  was  developed  by  the  electronics  laboratory  of  the  von 
Karrnan  Institute.  This  enabled  many  investigations  to  be  undertaken  with  the  laser-Doppler  velocimeter. 
Riethmuller  and  Olivari  (13)  have  performed  velocity  measurtuients  in  wet-steam  turbine  cascades.  Measure- 
ments in  centrifugal  pumps  and  in  two-phase  flows  have  also  been  made  with  this  equipment.  A more  recent 
analysis  of  the  application  of  the  laser-Doppler  velocimeter  to  unsteady  flow  measurements  has  been  per- 
formed by  Riethmuller  (14).  Figure  6 shows  a view  of  the  VKI  laser  velocimeter  set  up  for  measurements 
in  a free  jet  (see  also  A. 3,  C.7  and  E.9). 


6.  CONCLUSION 


Work  on  gas-particle  flow  has  been  in  progress  at  the  VKI  only  since  1970.  The  preceding  sections 
describe  the  status  of  this  program  at  the  present  time,  and  work  in  this  area  is  continuing. 
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E.  FUNDAMENTALS  AND  APPLICATIONS  OF  GAS-PARTICLE  FLOW 
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Buffalo,  New  York 


SUMMARY 

Tills  survey  deals  with  flows  of  a gas  in  which  small, rigid,  and  permanent  particles  are  suspended. 
Particle  concentrations  range  from  so  low  that  the  particles  do  not  affect  the  gas  flow  and  can  be  treated 
as  single  particles  to  so  high  that  the  particles  occupy  an  apprecialbe  volume  fraction  of  the  mixture. 
However,  dense-phase  transport  and  fluidized  beds  are  not  discussed.  A number  of  important  applications 
are  outlined.  The  discussion  covers  the  dynanics  of  single  particles  in  continuum  and  low-density  flow, 
the  thermodynamics  of  gas-particle  mixtures,  and  the  basic  flow  equations  for  one -dimensional  flow.  Wave 
propagation,  nozzle  flow  and  gas-particle  Jets  are  discussed  in  some  detail.  Additional  examples  of  ana- 
lytical and  experimental  results  are  given  to  illustrate  important  situations.  Finally,  a brief  discussion 
is  given  on  problems  for  which  satisfactory  solutions  have  not  yet  been  obtained.  Emphasis  is  placed  on 
work  performed  in  the  United  States. 
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LIST  OF  SYMBOLS 
cross-sectional  area  of  duct 
speed  of  sound  in  the  gas  phase 

frozen  and  equilibrium  speed  of  sound  in  the  gas-particle  mixture 

drag  coefficient,  Eq.  (4) 

friction  coefficient  for  a flat  plate 

transformed  friction  coefficient 

specific  heat  of  particle  material 

specific  heat  of  gas  at  constant  pressjre  and  constant  volume 
particle  diameter 
duct  diameter 

internal  energy  of  the  gas,  the  oarticles  and  the  mixture 
effective  force  between  two  particle  clouds,  Eq.  (37) 
correction  to  the  Stokes  drag 

friction  factor  for  a pipe  for  flow  of  pure  gas  and  of  gas  with  solid  particles 

Cunningham  correction  for  the  drag  coefficient 

enthalpy  of  the  gas,  the  particles , a^d  the  mixture 

heat-transfer  coefficient,  Eq.  (15) 

momentum-flux  ratio,  Eq.  (47) 

Knudsen  number 
velocity  ratio 

complex  wave-propagation  parameter 
thermal  conductivity  of  gas  phase 
length  of  nozzle 

characteristic  length  for  the  interaction  between  two  particle  clouds,  Eq.  (38) 
Mach  number  of  the  gas  flow 
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Mach  number  of  the  gas-particle  mixture,  based  on  the  frozen  and  equilibrium 
speed  of  sound 

mass  flow  of  gas  and  particle  phase 
Nusselt  number 
Prandtl  number 
pressure 

reservoir  pressure 

gas  constant  for  gas  phase  and  equilibrium  mixture 

Reynolds  number,  defined  by  Eq.  (5)  unless  otherwise  noted 

Reynolds  number  of  gas  flow  In  a pipe 

gas  and  particle  temperature 

reservoir  temperature 

time 

velocity  of  flat  plate 

gas  and  particle  velocity  components  in  direction  x 
gas  and  particle  velocity  components  in  direction  y 
coordinate  in  the  main  flow  direction 
maximum  lateral  penetration 

coordinate  in  the  lateral  direction  to  the  main  flow 
thickness  of  particle-free  l^yc*  near  flat  plate 
equivalence  parameter  for  gas-particle  jets,  Eq.  (45) 

ratio  of  the  specific  heats  for  gas-phase  and  equilibrium  gas-particle  mixture 

relative  specific  heat 

particle  volume  fraction 

loading  ratio 

molecular  mean  free  path 

gas  viscosity 

gas  density 

density  of  particle  material 
average  density  of  gas-particle  mixture 
concentration  of  the  gas  and  particle  phases 
velocity  relaxation  time,  Eq.  (9) 
temperature  relaxation  time,  Eq.  (17) 
particle  mass  fraction,  Eq.  (19) 

correction  factor  for  to  allow  for  finite  particle  volume  fraction,  Eq.  (14) 
circular  frequency  of  sound  waves 
and  asterisks  (*)  are  defined  as  needed 


1.  INTRODUCTION 

Gas-particle  flows  represent  one  segment  of  the  general  class  of  two-phase  flows,  that  is,  flows  of 
mixtures  of  solid  and  liquid,  solid  and  gaseous,  or  liquid  and  gaseous  constituents.  Depending  on  the 
phases  involved  and  their  distribution  within  the  flow,  a large  variety  of  flow  patterns  can  be  observed 
The  present  discussion  deals  only  with  flows  in  which  small  particles  are  suspended  in  a gas.  This  more 


57 

limited  topic  still  covers  an  area  that  Is  far  too  extensive  for  a survey  such  as  this,  and  further  re- 
strictions must  be  made.  Some  of  the  choice:  were  Indicated  In  the  initial  Invitation  to  prepare  this 
review, and  others  were  influenced  by  the  writer's  personal  experience. 

Viscous  Interaction  and  convective  heat  transfer  between  a perfect  gas  and  rigid  permanent  particles 
will  be  discussed,  while  mass  transfer  by  evaporation,  condensation  or  chemical  reaction,  heat  transfer 
between  the  flowing  mixture  and  the  walls  of  a duct,  as  well  as  breakup  of  liquid  droplets  will  not  be 
considered.  Also  excluded  are  effects  of  electrostatically  charged  particles  (see  C.2,2).  The  particle 
concentrations  considered  vary  from  extremely  low,  so  that  the  gas  flow  is  not  affected  by  the  presence 
of  the  particles,  to  sufficiently  high,  that  not  only  modifications  cf  the  gas  flow  become  important,  but 
also  the  volume  occupied  by  the  particles  must  be  taken  into  account.  However,  the  special  problems  asso- 
ciated with  fluidized  beds  and  pipelines  for  conveying  of  powdered  materials  will,  not  be  specifically  dis- 
cussed (see  also  B.3.2  and  C.6).  Experimental  techniques  will  not  be  presented  as  such,  but  some  experi- 
mental data  are  included  to  illustrate  specific  flow  properties. 

Even  with  these  restrictions,  there  are  many  flows  of  practical  importance  which  require  sound  under- 
standing. An  enormous  literature  is  available,  as  indicated,  for  example,  by  the  number  of  references  in 
the  frequently  cited  books  by  Fuchs  (1)  and  Soo  (2),  each  of  which  lists  some  900  references.  The  older 
survey  by  Torobin  and  Gauvln  (3)  also  contains  nearly  500  references.  Thus,  the  references  given  here 
must  be  considered  as  representative,  and  no  attempt  was  made  to  be  all-inclusive.  Some  basic  aspects  of 
the  mathematical  analysis  are  briefly  outlined  in  the  present  review,  but  reference  to  the  original  liter- 
ature should  be  made  for  further  details.  In  accordance  with  the  invitation  to  prepare  this  review,  em- 
phasis is  placed  on  work  performed  in  the  United  States. 

Since  pr^.cical  applications  provide  the  major  motivation  for  the  study  of  gas-particle  flow,  the  more 
important  applications  are  indicated  in  the  next  section.  Subsequent  sections  deal  with  the  interaction 
of  a gas  with  individual  suspended  particles  and  with  particle  clouds,  the  thermodynamic  properties  of 
equilibrium  gas-particle  mixtures,  and  several  specific  flows  These  are  followed  by  a section  in  which 
numerical  or  experimental  examples  are  given  to  illustrate  the  flow  properties,  and  the  final  section  in- 
dicates some,  as  yet  unresolved,  problems. 


2.  APPLICATIONS  OF  GAS-PARTICLE  FLOW 

This  section  provides  a brief  overview  of  some  important  applications  of  gas-particle  flow.  Illustra- 
tive numerical  or  experimental  examples  are  presented  in  the  subsequent  sections. 

The  particle  concentration  in  flows  of  practical  interest  may  vary  from  extremely  low,  so  that  the 
gas  flow  is  not  affected  by  the  presence  of  the  particles,  to  sufficiently  high,  that  not  only  the  gas 
flow  is  modified  but  that  the  particles  occupy  a significant  fraction  of  the  volume  of  the  mixture. 

If  the  gas  flow  is  not  affected  by  the  particles,  the  latter  may  be  used  as  tracers  to  determine  local 
gas  motion.  An  old  and  well-known  technique  is  to  make  the  particles  (smoke)  visible  by  suitable  illumi- 
nation. A more  recent  development  is  the  laser-Doppler  technique,  e.g.,  Stevenson  and  Thompson  (4),  in 
which  the  particle  velocity  is  inferred  from  the  Doppler  frequency  shift  of  scattered  laser  light.  Because 
of  their  inertia,  the  ability  of  particles  to  follow  rapid  changes  of  the  flow  velocity  may  be  inadequate, 
and  it  becomes  important  to  analyze  the  accuracy  with  which  particles  of  a given  size  can  be  used  as  tracers 
in  a specific  flow.  An  example  of  this  application  is  given  in  Section  9,  Sufficiently  small  particles 
can  follow  even  turbulent  fluctuations,  and  the  laser-Doppler  technique  is  thus  suitable  for  turbulence 
studies  without  Interference  with  the  flow.  The  motion  of  suspended  particles  also  is  of  interest  for 
studies  of  pollution  by  particulate  matter  from  smoke  stacks  and  other  sources,  e.g..  Stern  (5),  Soo  (6). 

A change  in  the  conditions  at  some  point  of  a flow  is  "signaled"  to  other  points  by  waves  which  pro- 
duce the  necessary  flow  adjustments.  It  is  thus  important  to  understand  how  the  propagation  characteristics 
of  these  waves  depend  on  the  composition  of  the  medium  and  the  properties  of  the  phases  themselves.  Energy 
dissipation  associated  with  particle  slip  with  respect  to  the  gas  may  be  useful  for  attenuation  of  blast 
waves,  Kalavsky  and  Rudin  (7), or  noise, Temkin  and  Dobbins  (8,9).  In  his  review  of  pressure  waves,  Rudinger 
(10)  distinguishes  three  types:  (1)  shock  waves  which  produce  a transition  from  one  flow  condition  to 
another  of  higher  pressure;  (2)  sound  waves  which  are  periodic  disturbances  of  small  amplitude  and  fixed 
frequency;  (3)  large-amplitude  waves  of  arbitrary  wave  form  which  may  produce  a net  rise  or  drop  of  the 
pressure.  The  properties  of  these  waves  are  summarized  in  Section  6. 

As  long  as  the  particles  do  not  affect  the  gas  flow,  analysis  of  the  gas  flow  is  uncoupled  from  that 
of  the  particle  motion.  If  the  particle  concentration  becomes  high  enough  to  affect  the  gas  flow,  the 
analysis  becomes  much  more  complicated,  because  both  phases  of  the  flow  must  be  considered  at  the  same 
time. 

Many  solid  rocket  propellants  contain  a fine  metallic  powder.  These  particles  provide  heat  release 
by  combustion,  but  their  main  purpose  is  to  suppress  pressure  oscillations  in  the  combustion  chamber,  e.g., 
Horton  and  McGle  (11),  Dobbins  and  Temkin  (12).  After  combustion,  such  additives  may  form  condensed  solid 
particles  (e.g.,  aluminum  oxide)  which  must  be  exhausted  with  the  rest  of  the  combustion  products.  These 
particles  may  cause  additional  erosion  in  the  rocket  nozzle  by  Impact  on  the  wall  and  Increased  heat  trans- 
fer. Furthermore,  the  particles  lag  during  the  high  flow  accelerations  in  the  nozzle,  and  the  resultant 
viscous  drag  on  the  gas  flow  leads  to  some  degradation  of  performance.  Nozzle  flows  are  discussed  at  some 
detail  in  Section  7. 

In  a rocket  nozzle,  the  particles  may  represent  30-407,  of  the  total  mass  flow.  Much  higher  concen- 
trations, 907.  or  more,  of  the  total  flow  may  be  utilized  to  transport  granulate  materials  over  long  dis- 
tances through  conveying  pipelines.  This  application  has  developed  into  an  Important  technology  of  its 
own,  e.g.,  Spencer,  et  al.,  (13),  (Xtfen  (14),  Thornton  (15);  although  conveying  will  not  be  specifically 
discussed,  some  of  the  analytical  procedures  described  in  the  following  sections  are  relevant  to  it. 
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For  vehicles  flying  In  the  atmosphere  through  dust  or  rain  clouds.  It  Is  Important  to  know  the  Impact 
velocity  with  which  small  particles  hit  the  vehicle  surface.  Studies  of  particle  impact  go  back  to  Langmuir 
and  Blodgett  (16)  who  Investigated  the  Impact  of  water  droplets  on  cylinders  representing  the  leading  por- 
tion of  airplane  wings.  At  high  altitudes,  the  air  ccn  no  longer  be  treated  as  a continuum,  because  the 
molecular  mean-free-path  becomes  comparable  with  or  even  larger  than,  the  particle  dimensions.  However, 
at  hypersonic  speeds,  the  compression  by  the  leading  shock  may  raise  the  density  sufficiently  to  allow 
the  flow  around  the  nose  of  the  body  to  be  treated  as  a continuum.  On  this  basis,  Probsteln  and  Fassio 
(17)  computed  the  collection  efficiency  of  slender  wedges,  cones,  cylinders  and  spheres. 

Erosion  by  suspended  particles  is  important  also  for  turbomachinery  in  which  particles  are  deliberately 
included  in  the  flow  for  the  purpose  of  punning,  or  are  present  as  a result  of  ingestion  from  the  atmosphere, 
Spencer,  et  al.  (13),  Tabakoff,  et  al.  (18-22),  Clevenger  and  Tabakoff  (23),  Grant  and  Tabakoff  (24).  An  il- 
lustration of  particle  trajectories  in  a turbine  cascade  is  given  in  Section  9 (see  also  A. 4 and  C.4). 

The  inability  of  particles  to  follow  rapid  flow  changes  is  customarily  used  to  separate  particles 
from  a gas  flow  or  to  classify  particles  into  different  size  fractions,  Langmuir  and  Blodgett  (16),  Levin 
(25),  Golovin  and  Putnam  (26),  Lapple  (27).  Cyclones,  multiple-plate  impactors,  and  other  devices  are 
commercially  available  in  many  forms.  Since  particles  of  different  sizes  respond  to  flow  changes  at  dif- 
ferent rates,  collisions  between  particles  may  be  induced  which  can  modify  the  flow,  Marble  (28,29),  or 
may  be  of  interest  for  their  own  sake,  for  instance,  coalescence  of  cloud  droplets  under  the  influence 
of  thunder,  Temkin  (30). 

The  converse  problem  arises  if  particles  are  to  be  retained  in  a vessel  through  which  gas  is  flowing. 

In  a fluidized  bed,  gravity  provides  the  force  necessary  to  hold  the  particles  against  the  drag  of  an  up- 
ward directed  gas  flow.  Such  beds  have  become  of  extreme  importance  for  many  chemical  processes,  e.g.. 

Leva  (31),  Zenz  and  Othmer  (32).  An  extension  of  the  gravitational  fluidized  bed  is  a rotating  bed  in 
which  the  centrifugal  field  takes  the  place  of  gravity.  Such  schemes  have  been  proposed  to  contain  the 
active  material  of  a nuclear  reactor  and,  at  the  same  time,  extract  the  heat  energy  by  the  flowing  gas, 

Hatch, et.  al.  (33),  Jackomis  and  von  Ohain  (34),  Thrman  and  Hasinger  (35),  Ravets,  et  al.  (36).  As  pointed 
out  earlier,  special  problems  associated  with  fluidized  beds  are  not  discussed  in  this  review. 

A suggestion  that  gas  flows  of  extremely  high  total  enthalpy  might  be  produced  by  using  gas-particle 
flows  was  made  by  Johnson  and  von  Ohain  (37).  They  proposed  to  accelerate  a stream  of  solid  particles  to 
the  high  velocity  that  can  be  obtained  with  a heated  light  gas  (hydrogen  or  helium)  and,  after  separation 
of  the  particles  from  this  gas,  to  use  the  particles  for  the  acceleration  of  a working  fluid  (air)  to  a 
high  velocity.  Ihis  technique  would  produce  the  final  velocity  without  going  through  phases  of  extremely 
high  static  pressures  or  temperatures.  The  technical  feasibility  of  this  scheme  has  not  yet  been  demon- 
strated. 

The  foregoing  potential  application  would  involve  extensive  use  of  gas-particle  Jets.  Understanding 
of  such  Jets  is  important  also  for  atomizing  systems,  rocket  exhausts,  injection  of  powdered  fuels  into  a 
combustion  chamber  by  means  of  a carrier  gas  (or  liquid),  and  water-augmented  air  cushion  vehicles.  In  the 
last  mentioned  application,  it  was  proposed  to  obtain  thrust  augmentation  by  injecting  the  water  spray, 
produced  by  the  escaping  flow  from  the  air  cushion,  into  the  propulsive  duct,  Hope-Gill,  et  al.  (38);  in 
effect,  this  application  would  involve  injection  into  a cross  flow  through  a relatively  long  and  narrow 
longitudinal  slot. 

Although  a great  deal  is  known  about  gas  J ets , e.g. , Abramovich  (39),  Margason  (40),  Skifstad  (41), 
Kamotanl  and  Greber  (42),  only  little  information  is  available  on  gas-particle  Jets.  Abramovich  (39)  at 
first  treated  a two-phase  Jet  as  if  it  were  a gas  Jet  of  increased  density.  Subsequently  (43),  he  reported 
that  two  Jets  of  equal  average  density  injected  into  a quiescent  gas  develop  differently.  The  gas-particle 
Jet  seems  narrower  and  longer  ranged.  He  ascribed  this  different  behavior  to  the  effect  of  the  particles 
on  the  turbulence  of  the  flow. 

An  attempt  was  made  by  Soo  (2)  to  analyze  particle  diffusion  for  a Jet  issuing  into  a stagnant  gas. 

He  concluded  (on  page  359)  that  "diffusion  of  particles  is  completely  similar  to  the  diffusion  of  momentum 
of  the  fluid."  This  finding  seems  to  contradict  his  statement  on  page  83  that  the  particle  diffusivity 
is  only  1 to  107.  of  the  eddy  diffusivity. 

A few  experiments  on  lateral  injection  of  a gas-particle  Jet  into  a supersonic  cross  flow  were  reported 
by  Edelman,  et  al.  (44).  These  investigators  used  low  particle  flow  rates,  mostly  below  507.  of  the  Jet  gas 
flow.  Particle  trajectories  were  computed  for  two  extreme  cases.  In  one  limit,  the  particles  were  assumed 
to  be  in  equilibrium  with  the  gas;  in  the  other  limit,  the  gas-phase  trajectory  was  assumed  unaffected  by 
the  particles,  and  the  particle  trajectory  was  computed  for  various  drag  coefficients.  The  results  differed 
remarkedly  for  different  assumptions.  Experimentally,  they  found  that  the  particles  remain  "locked"  with 
the  gas-phase  stream,  and  concluded  that  "reasonable  estimates  of  Jet  penetration  can  be  obtained  using 
existing  empirical  relations."  Salzman  (45)  injected  gas-particle  Jets  laterally  into  an  air  stream  and 
determined  both  the  centerline  of  the  trajectory  and  Jet  spreading  in  the  plane  of  the  trajectory  with  an 
Isokinetic  sampling  probe.  The  turbulence  level  of  the  cross  stream  had  little  influence  on  the  Jet  tra- 
jectory but  had  a marked  effect  on  dispersion.  He  observed  some  separation  of  the  centerline  of  the  particle 
stream  from  that  of  the  carrier  gas.  Rudinger  (46)  also  studied  lateral  injection  into  a flow  and  observed 
complete  separation  of  the  particles  from  the  carrier  gas.  These  effects  are  closely  related  to  the  size  of 
the  injected  particles  and  are  further  discussed  n.  Section  8. 

Almost  all  flows  of  technical  importance  are  turbulent.  It  is  therefore  of  considerable  significance 
that  the  friction  coefficient  for  gas-particle  flow  may  be  considerably  lower  than  that  for  pure  gas  flow. 

The  existence  of  this  effect  has  been  known  for  some  years,  and  has  variously  been  ascribed  to  electro- 
static effects,  Richardson  and  McCleman  (47),  reduction  of  the  gas  viscosity,  Sproull  (48),  thickening  of 
the  laminar  sublayer,  Peskln  and  Dwyer  (49),  and  general  turbulence  suppression,  Soo  (2).  These  and  other 
related  studies  were  reviewed  and  extended  by  Kane  and  Pfeffer  (50)  and  Pfeffer  and  Kane  (51)  who  concluded 
that  the  thickening  of  the  laminar  sublayer  is  the  predominant  effect.  With  particle  sizes  between  10  and 
60  pm  and  particle  flow  rates  up  to  four  times  the  gas  flow  rate,  they  found  a reduction  of  the  friction 
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pressure  drop  by  about  40%  under  optimum  conditions.  A few  of  their  results  are  discussed  in  Section  9. 


3.  GAS-PARTICLE  INTERACTION 
3.1  Viscous  Drag  of  Single  Particles 

Particles  suspended  in  a gas  flow  move  under  the  influence  of  forces  exerted  on  them  by  the  gas,  but 
if  the  particle  concentration  is  low  enough,  the  reaction  of  these  forces  affects  the  gas  flow  only  to  a 
negligible  extent.  Interactions  of  particles  with  wakes  of  other  particles  then  also  are  of  no  consequence. 
This  condition  is  satisfied  if  the  average  spacing  between  particles  is  about  50  or  mare  particle  diameters. 
The  particles  then  represent  not  more  than  about  1%  of  the  mass  of  the  mixture  (Section  4). 


The  equation  of  motion  for  a spherical  particle  in  a gas  flow,  which  Itself  need  not  be  constant, 
is  given,  e.g,  Hinze  (52),  Soo  (2),  Rudinger  (53), by 
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where  the  substantial  derivatives  must  be  taken  along  the  trajectories  of  the  gas  and  of  the  particles,  so 
that 
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Lateral  forces  on  the  particles  may  develop  as  a result  of  particle  rotation,  slip  in  the  shear  flow 
of  a boundary  layer,  ar d other  effects.  Such  forces  will  not  be  discussed  here, but  some  reference  to  them 
is  made  in  Sections  3.2  and  9. 


The  left-hand  side  of  Eq.  (1)  reprrsents  the  product  of  mass  and  acceleration,  and  the  right-hand  side 
Indicates  various  forces  that  act  on  the  particle.  The  viscous  drag  is  discussed  below.  The  second  term 
shows  the  effect  of  the  pressure  gradient  in  the  gas  for  which  the  gas  acceleration  can  be  substituted 
Hinze,  (52),  because  of  the  momentum  equation 
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If  the  gas  flow  is  affected  by  the  particles,  Eq.  (3)  does  not  hold  (Section  5),  and  the  substitution  then 
is  not  permissible.  The  third  term  indicates  the  force  needed  to  accelerate  the  added  mass  which,  for 
spherical  particles,  equals  one-half  of  the  mass  of  the  displaced  fluid.  The  integral  term,  also  known  as 
Basset  force,  accounts  for  deviation  of  the  flow  pattern  around  the  particle  from  that  for  steady  flow; 
it  represents  the  influence  of  the  history  of  the  motion  on  the  instantaneous  force.  Under  the  usual  con- 
ditions, in  which  p„  is  about  three  orders  of  magnitude  greater  than  0,  only  the  viscous  drag  needs  to  be 
considered,  Hinze  (52),  Fuchs  (1),  Hjelmfelt  and  Mockros  (54).  However,  the  effect  of  the  history  of  the 
motion  may  affect  the  particle  trajectory  under  extreme  conditions.  Lewis  and  Gauvin  (55)  injected  particles 
into  the  flow  of  an  argon  plasma  torch  and  determined  their  motion  by  streak  photography.  They  measured 
the  residence  time  between  leaving  the  torch  exit  and  arrival  at  a distance  of  178  mm.  These  results 
were  compared  with  calculated  residence  times  based  on  tne  best  available  drag  data,  Beard  and  Pruppacher 
(56)  (see  Table  1).  Under  these  extreme  conditions,  the  effect  of  acceleration  was  to  increase  the  resi- 
dence time  by  about  10%  for  30-^m  particles  and  30%  for  140-p.m  particles.  Their  preliminary  recommendation 
was  to  use  the  acceleration  terms  with  correction  factors  previously  determined  by  Odar  and  Hamilton  (57) 
for  accelerated  motion  of  spheres  in  a liquid. 

Viscous  drag  is  a complicated  function  of  the  flow  conditions  which,  in  general,  can  be  described 
only  by  empirical  relationships  based  on  experimental  data.  The  drag  force  for  a spherical  particle  can 
be  written  as 

2 

Viscous  Drag  = CD  • p | u - upj(u  - u ) (4) 


which  defines  the  force  in  terms  of  a drag  coefficient,  the  frontal  area  of  the  particle,  and  the  dynamic 
head  of  the  flow  relative  to  the  particle.  The  square  of  the  relative  velocity  is  written  in  this  form 
to  Insure  the  correct  sign  of  the  force.  In  continuum  flow,  the  drag  coefficient  is  a function  only  of 
the  particle  Reynolds  number 


(5) 


For  low  Reynolds  numbers,  the  well-known  classical  Stokes  solution 
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applies,  Hie  drag  then  takes  the  simple  form 
Stokes  Drag  = 3nDp.(u  - u^) 


(7) 


For  higher  Reynolds  numbers,  extensive  experimental  data  are  available.  The  relationship  between  C 
and  Re  is  shown  in  Fig.  1,  e.g.,  Schlichting  (58),  and  is  usually  referred  to  as  "standard  drag  coefficient." 
The  agreement  of  data  from  different  sources  appears  quite  good,  but  note  the  discussion  of  Table  1 below. 
Line  1 in  the  figure  represents  Stokes  theory  of  Eq.  (6)  and  line  2 the  extension  derived  by  Oseen.  Many 
empirical  fits  to  the  standard  drag  curve  have  been  proposed,  e.g,  Torobin  and  Gauvin,  Part  I (3),  Fuchs 
(1),  Soo  (2).  The  equation  of  motion  for  a single  particle  then  becomes 
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and  f is  a correction  factor  to  the  Stokes  drag.  Since  f is  a function  of  the  Reynolds  number,  Eq.  (8)  can, 
in  general,  be  integrated  only  numerically.  For  Stokes  drag  (f  = 1)  and  for  a constant  gas  velocity,  Eq. 

(8)  indicates  that  the  velocity  difference  u - up  decays  exponentially  with  a relaxation  time  tv  given  by 
(9).  For  larger  Reynolds  numbers,  rv  merely  becomes  a convenient  qualitative  measure  of  velocity  re- 
laxation (see  also  B 2.2  and  C.2.1). 

The  prpsent  survey  deals  with  small  particles,  usually  smaller  than  100  ^m.  The  Reynolds  number  then 
rarely  exceeds  several  hundred.  For  this  range,  the  relationship 
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represents  a good  fit  to  the  standard  drag  curve.  It  was  first  proposed  by  Klyachko  in  1934  (l).  Apparently 
independently,  it  was  later  suggested  by  Putnam  (59)  who  pointed  out  that  it  may  lead  to  integrable  equa- 
tions, For  most  applications,  these  data  are  sufficiently  accurate,  but  situations  may  arise  when  a higher 
accuracy  is  needed,  for  instance,  the  study  by  Lewis  and  Gauvin  (55)  referred  to  earlier.  Maxworthy  (60) 
pointed  out  that  for  very  low  Reynolds  numbers,  it  is  "impossible  to  be  sure  of  the  drag  to  better  than 
+ 20%,  even  if  one  is  very  generous  with  data  that  scatters  between  the  Stokes  drag  and  values  above  the 
Oseen  drag."  He,  and  subsequently  Beard  and  Pruppacher  (56),  performed  more  accurate  measurement s.  Hie 
latter  are  considered  by  Lewis  and  Gauvin  the  most  precise  ones  currently  available.  A comparison  of  the 
standard  drag  coefficient,  as  given  in  a table  by  Fuchs  (1)  with  Eqs.  (4)  and  (10),  and  some  of  the  data 
by  Beard  and  Pruppacher  is  shown  in  Table  1. 


Table  1 

Comparison  of  Particle  Drag  Coefficients 


Reynolds 

Number 

Standard 

Drag  Coefficient 
Eq.  (10)  Stokes 

B6.P* 

0.01 

2400 

2420 

2400 

0.1 

243 

249 

240 

242.7 

1 

26.9 

28.0 

24.0 

26.45 

10 

4.33 

4.26 

2.40 

4.149 

100 

1.09 

1.10 

0.24 

1.073 

500 

0.568 

0.552 

0.048 

1000 

0.469 

0.424 

0.024 

2000-10000 

0.40 

* Beard  and  Pruppacher  (56) 


This  table  indicates  that  Eq.  (10)  represents  the  standard  drag  coefficient  to  better  than  5%  for 

Reynolds  Numbers  up  to  about  500  and  to  better  than  10%  even  at  Re  * 1000.  Hie  agreement  with  the  data  by 

Beard  and  Pruppacher  is  not  quite  so  good.  In  contrast,  Stokes  drag  is  a good  approximation  only  for 

Reynolds  numbers  less  than  one.  Nevertheless,  it  is  frequently  used  for  qualitative  calculations,  because 

it  leads  to  particularly  convenient  relationships. 

Since  the  viscous  drag  and  the  gas  and  particle  velocities  are  vectors,  Eq.  (8)  represents  only  one 
component  of  a vector  equation,  and  analogous  equations  could  be  written  for  the  other  components  (see 
also  C.2.1).  Often,  flows  are  treated  as  quasi-one-dimensional,  and  the  question  of  other  components  does 
not  arise.  It  should  be  pointed  out,  however,  that  the  drag  vector  has  the  direction  of  the  vector  dif- 
ference of  the  gas  and  particle  velocities  and  that  the  absolute  value  of  this  difference  must  be  used  in 
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calculating  the  Reynolds  number  and  the  drag  force.  Thus,  if  u,  up,  v and  vp  represent  the  components  of 
the  gas  and  particle  velocities  in  two  dimensions,  the  absolute  value)  u - up | in  Eqs.  (4)  and  (5)  must  be 
replaced  by  Ku  - u„)2  + (v  - v )7]^.  Consequently,  the  same  value  of  f appears  in  all  component  equations 
and  the  gas  velocity  in  one  direction  affects  the  component  of  the  drag  in  a direction  at  right  angles 
to  it.  Only  for  a constant  value  of  f,  as  for  Stokes  drag,  do  the  component  equations  become  uncoupled 
and  can  be  solved  Independently.  This  coupling  of  the  equations  can  have  significant  consequences  as 
demonstrated  in  Section  3.2  and  9. 

Throughout  the  preceding  discussion  it  is  implied  that  the  gas  flow  can  be  treated  as  a continuum.  If 
the  particles  are  extremely  small,  or  the  gas  density  is  quite  low,  the  molecular  mean  free  path  in  the  gas, 
X,  may  become  comparable  with  the  particle  dimensions,  and  it  then  becomes  necessary  to  apply  a correction 
to  the  drag  coefficient.  This  correction,  sometimes  called  "Cunningham  correction"  can  be  expressed  in 
terms  of  the  Knudsen  number  which  is  defined  as  Kn  = X/D.  To  correct  the  drag  coefficient,  the  latter  is 
divided  by  G which  is  given  by,  e.g.,  Fuchs  (1), 

G - 1 + Kn  (2.492  + 0.84  e-1'74^)  (Ha) 

It  follows  from  kinetic  theory  that  the  Knudsen  number  is  a function  of  the  ratio  of  the  Mach  number 
of  the  flow  to  the  particle  Reynolds  number,  Schaaf  and  Chambr^  (61).  Experimentally  determined  drag 
coefficients  therefore  can  also  be  presented  in  terms  of  Mach  number  and  Reynolds  number.  Carlson  and 
Hoglund  (62),  after  reviewing  available  data  proposed  the  relationship 

G - 1 + (M/Re)(3.82  + 1.28  e‘1,25  Re/M)  (lib) 

Beckwith  and  Bushnell  (63)  also  derived  correlations  to  represent  the  then  available  data.  They 
suggested  relationships  with  constants  that  are  changed  for  various  ranges  of  the  Mach  and  Reynolds  num- 
bers. Details  of  this  correlation  and  a table  of  the  required  constants  may  be  found  in  a paper  by  Yanta 
(64).  The  most  recent  measurements  appear  to  be  the  ballistic-range  data  by  Bailey  and  Hiatt  (65)  which 
cover  the  entire  range  from  continuum  to  near- free-molecular  flow,  free-strean  Mach  numbers  from  0.1  to 
20,  and  Reynolds  numbers  from  0.2  to  10^.  A sunmary  of  their  data  is  presented  in  Fig.  2.  Korkan,  et  al. 
(66)  modified  the  Beckwith  and  Bushnell  correlation  to  allow  for  different  gas  and  particle  temperatures 
and  to  improve  the  fit  of  the  experimental  data.  Still  further  improvements  are  described  in  a recent 
paper  by  Henderson  (67)  (see  also  B.2.1). 


3.2  Viscous  Drag  of  Particle  Clouds 

If  the  particle  concentration  is  sufficiently  high,  the  reaction  of  the  sum  of  the  drag  on  all  par- 
ticles becomes  large  enough  to  affect  the  gas  flow.  Furthermore  the  force  on  each  particle  may  be  affected 
by  the  presence  of  other  particles.  Direct  collisions  between  particles  of  equal  size  are  infrequent,  ex- 
cept at  the  high  concentrations  of  fluidized  beds  or  dense-phase  conveying.  However,  suspensions  of  uni- 
formly sized  particles  are  rare,  and  collisions  between  particles  of  different  size  then  may  result,  be- 
cause they  follow  different  trajectories.  Also  the  volume  of  a particle  wake  may  be  two  or  three  orders 
of  magnitude  larger  than  the  volume  of  the  particle  itself.  Thus,  particle-wake  interactions  may  be  sig- 
nificant, even  if  the  volume  occupied  by  the  particles  is  negligible.  The  resultant  uncertainty  of  the 
drag  coefficient  was  emphasized  by  Torobin  and  Gauvin  (3)  and  Hoglund  (68). 

Since  no  generally  valid  methods  are  available  to  evaluate  deviations  from  the  standard  drag  coeffi- 
cient, it  becomes  necessary  to  derive  effective  drag  coefficients  from  experimental  observations.  The 
results  of  several  such  studies  were  collected  and  discussed  by  Torobin  and  Gauvin,  Part  V (3), and  by 
Soo  (2).  The  data  collected  by  Soo  are  shown  in  Fig.  3,  and  the  contrast  between  this  figure  and  Fig.  1 
is  striking.  While  some  of  the  scatter  of  the  data  may  be  the  result  of  experimental  errors,  the  major 
trends  are  probably  real.  Torobin  and  Gauvin  stated  that  "the  use  of  the  standard  drag  curve  for  drag 
estimation  to  the  second  figure  ...  outside  the  Stokesian  region  is  not  valid." 

Ingebo  (69)  measured  the  acceleration  of  spray  droplets  injected  into  a wind  tunnel  and  found 

CD  " o?84  (6  < Re  < 400)  (12) 

shown  as  line  14  in  Fig.  3.  He  believed  that  the  deviations  from  the  standard  drag  were  caused  by  the 
acceleration  of  the  particles. 

Crowe,  et  al.  (70)  used  high-speed  movies,  and  Selberg  and  Nicholls  (71)  employed  multiple-flash 
photography,  to  observe  the  particle  motion  in  a shock  tube.  Both  investigations  yielded  higher  than 
standard  drag  coefficients.  The  difference  was  tentatively  ascribed  to  surface  roughness  of  the  particles. 

Rudlnger  (72,73)  used  streak  photography  and  light  scattering  to  determine  particle  motion  in  a shock 
tube.  Both  techniques  led  to  the  same  result,  shown  as  line  R48  in  Fig.  3, 
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for  29-t*m  particles  and  values  about  607.  higher  for  62-um  particles.  Other  features  of  the  results  are 
that  Cq  approaches  the  standard  drag  curve  for  higher  Reynolds  numbers  and  that  it  is  Independent  of  the 
particle  loading  ratio  over  the  range  of  the  experiments  (0.05  - 0.36).  To  explain  these  results,  turbu- 
lence, wall  effects,  electrostatic  forces,  particle  rotation,  or  asymmetric  vortex  shedding  were  consid- 
ered. With  the  possible  exception  of  turbulence,  none  of  these  effects  should  be  expected  to  have  a large 
direct  effect  on  the  drag  coefficient,  but  all  of  them  tend  to  produce  longitudinal  and  lateral  pertur- 
bations of  the  particle  trajectories.  Such  perturbations  are  too  small  to  be  observed,  and  the  effective 
drag  coefficient  is  thus  derived  from  the  average  trajectory.  Even  if  the  standard  drag  coefficient  is 
applicable  to  the  microscopic  particle  motion,  the  apparent,  or  effective,  drag  coefficient  derived  from 
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the  average  motion  is  different  because  of  the  nonlinearity  of  the  equation  of  motion  outside  the  rangc- 
of  Stokes  drag.  On  the  basis  of  this  consideration,  Rudinger  (73)  could  qualitatively  explain  all  fea- 
tures of  the  shock-tube  data. 


The  affect  of  particle  slip  appears  in  the  gas  initially  only  in  the  particle  wakes.  Since  mixing 
of  these  wakes  with  the  rest  of  the  gas  is  not  instantaneous,  neighboring  particles  may  not  be  immediately 
affected.  This  mixing  delay  also  may  influence  the  effective  drag  coefficient,  Rudinger  (74).  Ihe  effect 
of  a particle  cloud  on  the  drag  of  individual  particles  must  depend  also  on  the  parLjcle  volume  fraction 
e.  A statistical  analy  is  of  Tam  (75)  for  small  slip  velocities  yields  a correction  factor  for  Stokes 
drag  given  by 


Y 


4 + 3<  + 3 (8e  - 3c2)^ 
(2  - 3e)2 


(14) 


For  the  concentrations  used  in  the  foregoing  shock  tube  experiments,  the  part.cle  volume  fraction  is  less 
than  10“3  and  this  correction  factor  then  yields  less  than  1.1.  Therefore,  it  cannot  explain  the  obser- 
vations, but  the  fact  that  the  derivative  dY/de  approaches  infinity  as  c approaches  zero  indicates  that 
even  small  particle  concentrations  may  have  a significant  effect  on  the  drag  coefficient. 


All  phenomena  mentioned  in  the  foregoing  may  affect  the  particle  motion  in  a manner  that  is  unknown 
at  this  time  and  also  may  depend  on  the  particular  flow  system.  Ihey  contribute  to  the  uncertainty  of 
the  drag  coefficient,  but  their  consequences  are  implied  in  effective  drag  coefficients  tn«t  art  derived 
from  experimental  observations.  This  reasoning  may  help  to  explain  the  wide  scatter  of  the  experimental 
data  shown  in  Fig.  1 (see  also  C.6). 


3.3  Gaj-Particle  Heat  Transfer 


In  compressible  flow,  any  velocity  change  is  accompanied  by  a temperature  change.  Thus,  both  the 
velocity  and  the  temperature  of  suspended  particles  can  lag  behind  the  corresponding  gas  values.  The 
rate  at  which  the  particle  temperature  can  follow  the  gas  temperature  is  described  by  the  equation 
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(15) 


where  the  heat- transfer  coefficient  h can  be  conveniently  expressed  in  terms  of  the  Nusselt  no  * 

Nu  * hD/k.  For  pure  heat  conduction,  the  value  Nu  * 2 applies,  but  additional  heat  transfer  b>  ernvec- 
tion  takes  place  if  the  particle  moves  with  respect  to  the  gas.  The  Nusselt  number  then  becomes  a func- 
tion of  the  particle  Reynolds  number.  A frequently  used  relationship  is,  e.g.  , Knudsen  and  Katz  '/6) 

Nu  = 2.0  + 0.6  Pr1/3  Rel/2  (16) 

where  Pr  = uCp/k  Is  the  Prandtl  number.  If  the  particle  temperature  were  high  enough  for  radiative  heat 
transfir  to  become  important,  a more  elaborate  heat-transfer  coefficient  would  be  needed,  Simmons  and 
Spadaro  (77),  Soo  (2)  (see  also  B.2.4). 


If  Nu  « 2 is  assumed  and  T is  constant,  Eq. 
decreases  exponentially  with  a time  constant 


(15)  indicates  that  the  temperature  difference  T - Tp 
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where  the  last  form  is  obtained  by  introducing  Eq.  (9)  and  the  definition  of  the  Prandtl  number  (see  also 
B.2,2).  Since  6 * c/cp  is  of  order  one,  and  the  Prandtl  number  of  many  gases  is  close  to  2/3,  the  time 
constants  for  velocity  and  temperature  relaxation  often  are  approximately  equal,  Marble  (78),  Rudinger 
(53). 


As  in  the  case  of  velocity  relaxation,  t-j.  becomes  merely  a convenient  qualitative  indication  of  the 
rate  of  temperature  relaxation  if  the  particle  moves  with  respect  to  the  gas  (see  also  B.2.2). 

The  foregoing  relationships  imply  that  the  temperature  is  uniform  throughout  the  particle,  but  heat 
flux  into,  or  out  of,  the  particle  requires  internal  temperature  gradients.  Uniform  temperature  can  be 
assumed  if  the  rate  of  temperature  equalization  within  the  particle  is  fast  compared  with  the  rate  of 
heat  transfer  between  the  particle  and  the  gas.  Rudinger  (53)  has  shown  that  this  condition  is  almost 
always  satisfied,  except  for  particles  of  very  low  thermal  conductivity  immersed  in  a gas  of  very  high 
thermal  conductivity,  such  as  magnesium  oxide  particles  ir  hydrogen. 

So  far,  only  heat  exchange  with  single  particles  has  been  considered.  In  the  case  of  particle  clouds, 
uncertainties  arise  for  the  same  reasons  as  for  the  drag  coefficient.  The  tew  available  experimental  data 
were  summarized  by  Hoglund  (68)  with  the  statement  that  a "line  through  these  data  lies  about  a factor  of 
5 below  the  single-particle  data  at  Re  10,  crosses  the  single-particle  curve  at  Re  50  and  rises  to 
about  twice  the  single-particle  curve  at  Re  100," 


4.  THERMODYNAMIC  PROPERTIES  OF  GAS- PARTICLE  MIXTURES 

Before  discussing  the  analysis  of  gas-particle  flow,  it  is  useful  to  consider  the  general  thermodynamic 
properties  of  such  mixtures  which  enter  into  all  flow  calculations.  Derivation  of  these  quantities  can  be 
found  scattered  through  many  publications,  e.g.,  Soo  (79),  Hoglund  (68),  Marble  (29,78),  Kliegel  (80), 
Krlebel  (81),  Rudinger  (53,82),  Kralko  and  Sternin  (83).  In  deriving  these  properties,  careful  distinc- 
tion must  be  made  between  the  density  of  the  particle  material  pp  and  the  amount  of  particle  material  in  a 
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unit  volume  of  the  mixture  - the  particle  concentration  - denoted  by  op.  Similarly,  there  la  a distinc- 
tion between  the  gas  density  p and  the  gas  concentration  a.  The  densities  and  concentrations  are  related 
by 


o - «p  and  c » (1  - *)p  (18) 

P P 

where  « is  the  volume  fraction  occupied  by  the  particles.  Hie  mass  fraction  of  the  particles  therefore 
is  given  by 


Combination  of  these  equations  yields 


_f_  „ _B  £_  = _2_  fi_ 

!-•  ° °p  pp 


°M  " 0 + CTp  ’ 

An  alternative  form  of  this  equation  Is 


t ‘ (W 

(Although  all  publications  make  or  Imply  the  distinction  between  density  and  concentration,  the  terminology 
and  symbols  used  are  not  uniform.  Care  is  therefore  needed  in  the  interpretation  of  a particular  nomen- 
clature.) 

If  the  mixture  is  flowing  through  a duct  of  cross  section  A,  the  local  mass  flows  are  given  by  m - ouA 
and  iip  = apupA.  These  flow  rates  are  constant  only  if  the  flow  is  steady.  Their  ratio,  known  as  the  load- 
ing ratio,  then  becomes 
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It  is  clear  from  these  equations  that  t)  is  constant  for  steady  flow,  but  tp  is  constant  only  if  up/u 
is  also  constant.  Furthermore,  * is  constant  only  if  ap  is  constant,  because  the  particles  may  be  consid- 
ered as  incompressible. 

Under  usual  conditions,  the  density  ratio  p/pp  is  of  the  order  of  10”3,  Ihe  particle  volume  fraction 
thus  becomes  significant  only  for  rather  high  particle  loadings  and  can  be  neglected  for  low  or  moderate 
loadings.  Note  that  vanishing  particle  volume  implies  mathematically  that  p goes  to  infinity  and  ( to 
zero  in  such  a manner  that  the  product  (pp  remains  equal  to  op. 

Ihe  average  distance  between  neighboring  particles  is  of  the  order  D/s*^.  For  particle  diameters 
less  than  0.1  mm,  mass  fractions  as  low  as  0.03,  and  a density  ratio  p/pp  of  the  order  of  10“^  the  average 
particle  spacing  would  be  less  than  about  3mm.  Consequently,  both  the  particles  and  their  average  separa- 
tions generally  are  quite  small  compared  with  the  dimensions  of  typical  flow  fields.  Thus,  it  seems  natu- 
ral to  treat  the  particles  as  species  of  heavy  molecules  with  a molecular  weight  that  is  several  orders  of 
magnitude  larger  than  that  of  the  gas,  even  for  particles  as  small  as  0.01  p,m.  For  any  mixture  composition 
of  practical  interest,  the  number  density  of  the  particles  thus  is  insignificant  in  comparison  with  that 
of  the  gas  molecules,  and  the  contribution  of  the  particles  to  the  pressure  of  the  mixture  is  negligible. 
The  pressure  of  a gas-particle  mixture  therefore  is  that  of  the  gas  phase  alone  which,  for  a perfect  gas, 
is  given  by 

p - pRT  (24) 

To  obtain  an  equation  of  state  for  the  mixture,  the  gas  density  in  Eq.  (24)  must  be  expressed  in  terms 
of  the  density  of  the  mixture,  and  one  obtains,  with  the  aid  of  Eq.  (21) 


0M(l-<p) 


where  “ (l-tp)R  is  the  gas  constant  of  the  mixture.  Clearly,  Eq.  (25)  is  not  the  equation  of  state  of  a 
perfect  gas  because  of  the  density  term  in  the  denominator  which  is  directly  related  to  the  particle  volume 
fraction.  A mixture  of  solid  particles  with  a perfect  gas  may  be  treated  as  a perfect  gas  with  modified 
thermodynamic  properties  only  if  the  particle  volume  can  be  neglected.  This  result  is  a direct  consequence 
of  the  fact  that  the  particle  volume  does  not  change  during  compression  of  the  mixture. 

The  internal  energy  is  given  by  the  properly  weighted  average  of  the  internal  energy  of  the  consti- 
tuents as 


(1-<P)E  + q£ 


“ (1-<P>CVT  + <PcTp 


The  enthalpy  then  follows  from  Its  definition  and  Eq.  (22)  as 

* Em  + <c„T  + ? ) + <P(=T  + f"  ) 

M M pM  v p P Pp 

(27) 

■ (!*<p)cpT  + <P(cTp  + ) “ (1"V)H  + tpHp 

Thus,  the  enthalpy  of  the  particle  phase  depends  on  the  gas  pressure.  Since  «p/pp  is  of  order  «,  according 
to  Eq.  (21),  this  contribution  is  omitted  if  the  particle  volume  is  neglected. 

The  ratio  of  the  specific  heats  of  an  equilibrium  mixture  (u  = Up)  is  obtained  as  the  ratio  of  the 
properly  weighted  specific  heats  at  constant  pressure  and  constant  volume.  It  is  most  conveniently  ex- 
pressed in  terms  of  the  loading  ratio  which,  for  equilibrium  flow,  is  related  to  the  particle  mass  frac- 
tion by  11  • tp/  (l-<p) , and  one  obtains 


VM 


1+yT|6 


(28) 


This  result  shows  that  the  specific  heat  ratio  of  a mixture  is  always  smaller  than  that  of  the  gas  phase 
and  approaches  one  as  the  loading  ratio  is  increased,  A specific-heat  ratio  of  one  implies  isothermal 
flow.  In  this  case,  the  heat  capacity  of  the  particles  is  so  large  compared  with  that  of  the  gas  that 
temperature  changes  in  the  gas  during  compression  or  expansion  can  be  compensated  by  heat  transfer  to  or 
from  the  particles  without  significantly  changing  their  temperature.  Some  analyses  of  heavily  loaded  gas- 
particle  flow  can  be  considerably  simplified  by  assuming  isothermal  flow  (Section  7). 

The  last  property  to  be  discussed  here  is  the  speed  of  sound.  As  in  the  case  of  relaxing  gas  flows, 
one  must  distinguish  between  the  frozen  speed  of  sound  aj  and  the  equilibrium  speed  ae.  In  frozen  flow, 
no  viscous  interaction  or  heat  transfer  between  the  gas  and  the  particles  takes  place,  and  aj  thus  becomes 
equal  to  the  speed  of  sound  in  the  gas  phase 

af2  - a2  = yRT  (29) 


Hie  equilibrium  speed  of  sound  has  been  derived  in  different  ways  by  various  investigators  with  or 
without  taking  the  particle  volume  into  account,  but  all  results  can  be  reduced  to 

= (1-<P)[1  + 6a/  (!l-m)1 (30) 

^ 3 J (l-c)2[l  + v6<p/U-<p)] 

If  e is  neglected,  ae  approaches  zero  raonotonically  as  tp  is  increased,  but  at  sufficiently  high  par- 
ticle concentrations  the  influence  of  t becomes  significant,  and  ae  reaches  a minimum  value.  Still  higher 
concentrations  lead  to  a rapid  increase  of  ae  which  reaches  infinity  at  tp  * 1,  the  speed  of  sound  in  the 
incompressible  solid.  Ihis  behavior  is  demonstrated  in  Fig.  4 for  y = 1.4,  a density  ratio  o/Pp  = 0.01 
and  the  typical  value  6=1.  Parts  of  the  curves  are  shown  also  for  the  rather  extreme  values  6 = 0.05  and 
6 = 10.  This  figure  demonstrates  the  large  effect  of  the  particle  loading  on  ae  which  is  only  slightly 
modified  by  the  particle  volume  in  the  range  of  moderate  particle  loadings  considered  here. 


5.  GENERAL  FLOW  EQUATIONS 

The  motion  of  Isolated  particles,  which  do  not  affect  the  gas  flow,  is  discussed  in  Section  3.  To 
describe  flows  with  higher  particle  concentrations  requires  that  the  velocity,  concentration  and  tempera- 
ture of  the  gas  and  of  the  particles  be  determined  as  functions  of  the  position  and  time  coordinates.  The 
required  six  equations  are  represented  by  the  continuity,  momentum,  and  energy  equations  for  each  phase. 
por  three-dimensional  flows,  these  equations  can  be  written  in  tensor  notation.  Marble  (29,78),  or  in 
wctor  form,  Kralko  and  Sternln  (83),  Soo  (2).  Numerical  solutions  for  general  three-dimensional  flows 
would  require  major  computing  efforts,  and  investigations  usually  are  limited  to  cases  of  one-dimensional 
flow. 


Several  assumptions  must  be  made  to  derive  the  basic  flow  equations,  some  of  which  already  have  been 
mentioned  earlier: 

1.  The  gas  obeys  the  perfect  - gas  low,  and  the  specific  heats  are  constant; 

2.  The  particles  are  spherical,  of  uniform  size  (unless  specially  noted),  and  incompressible;  their 
specific  heat  is  constant,  and  the  temperature  is  uniform  within  each  particle; 

3.  The  particles  are  uniformly  distributed  over  the  cross-section  of  a duct,  and  their  size  and 
average  spacing  are  small  compared  with  the  cross-sectional  dimensions  of  the  duct; 

4.  Duct  flows  are  treated  as  one-dimensional,  so  that  changes  of  the  cross-sectional  area  must  be 
sufficiently  gradual.  Boundary- layer  effects  and  heat  exchange  with  the  walls  are  not  con- 
sidered (unless  specially  noted).  However,  viscosity  and  thermal  conductivity  of  the  gas  enter 
into  the  calculation  of  particle  drag  and  heat  transfer  between  the  particles  and  the  gas; 

5.  Both  the  drag  coefficient  and  the  Nusselt  number  must  be  prescribed  as  functions  of  the  Reynolds 
number  and,  if  necessary,  also  of  ths  Mach  number  or  Knudsen  number; 

6.  The  particles  do  not  contribute  to  the  pressure  and  do  not  collide  with  each  other  (unless 
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specially  noted) ; 

7,  In  the  flows  discussed  here,  no  external  forces  (such  as  gravity  or  electrostatic  fields)  or  mass 
transfer  between  the  gas  and  the  particles  are  considered. 

Derivation  of  the  flow  equations  follows  the  same  general  approach  as  for  a pure  gas.  Numerous  refer- 
ences can  be  found  for  one -dimensional  steady  flow  with  or  without  taking  the  particle  volume  into  account, 
e.g.,  Soo  (2,79),  Bailey  et  al.  (84),  Kliegel  (80),  Marble  (29,78),  Hassan  (85),  Rudinger  (86).  These 
equations  can  be  extended  to  nonsteady  flow,  Rudinger  and  Chang  (87),  Mlgdal  and  Agosta  (88),  Rudinger 
(53),  Wallis  (89). 


The  continuity  equations  for  the  gas  and  for  the  particles  are 


and 
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Hie  momentum  and  energy  equations  for  the  gas  and  for  the  particles  contain  gas-particle  interaction 
terms  with  opposite  signs.  It  is  therefore  convenient  to  eliminate  these  terms  once  by  adding  the  cor- 
responding equations.  The  momentum  equation  for  the  entire  mixture  thus  becomes 
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where  the  substantial  derivatives  have  the  meaning  previously  indicated  by  Eqs.  (2).  The  effect  of  the 
particles  is  not  directly  "felt"  over  the  entire  cross-section  of  the  duct  but  appears  initially  only  in 
the  particle  wakes.  Soo  (2,90)  suggested  that  the  momentum  term  for  the  particles  in  Eq.  (33)  therefore 
should  be  multiplied  by  an  "effectiveness  correction"  with  a value  between  zero  and  one.  Rudinger  (74) 
pointed  out  that  a violation  of  Newton's  laws  would  result  if  this  factor  were  assumed  constant,  because 
it  implies  that  the  viscous  drag  acts  on  the  particles  while  its  reaction  would  not  affect  the  entire 
gas;  he  suggested  that  the  delay  in  the  mixing  of  the  particle  wakes  with  the  rest  of  the  gas  might  be 
taken  into  account  by  letting  Soo's  effectiveness  correction  go  to  one  as  the  particles  approach  equili- 
brium with  the  gas.  Because  of  the  uncertainty  of  the  mixing  rate  of  the  wakes,  the  consequences  of  this 
delay  are  not  known,  but  they  are  implied  in  any  experimentally  determined  drag  coefficient  (section  3.2). 

The  energy  equation  of  the  mixture  is  given  by 
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The  remaining  two  equations  represent  the  gas-particle  interaction.  Particle  mc'.ion  is  described  by 
Eq.  (1).  Omitting  the  history  integral  and  terms  of  order  p/pp,  and  substituting  Eq.  (33)  for  the  pressure 
gradient,  yields  an  equation  analogous  to  Eq.  (8)  for  a single  particle 
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The  e-term  in  this  equation  appears,  because  the  pressure  gradient  now  depends  on  the  acceleration  of  both 
phases. 

The  relationship  of  gas-particle  heat  transfer,  after  substitution  of  Eq.  (17)  and  the  definition  of 
the  Nusselt  number  into  Eq.  (15),  becomes 
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The  system  of  equations  to  be  solved  is  thus  given  by  Eqs.  (31)- (36)  and  appropriate  assumption  for 
Cj)  and  Nu.  Since  the  gas  pressure  and  the  particle  volume  fraction  appear  as  additional  variables  in 
these  equations,  the  equation  of  state,  Eq.  (24),  and  Eqs.  (18)  are  needed  to  complete  the  system.  In  gen- 
eral, solution  of  these  equations  requires  a considerable  computational  effort.  Fortunately,  simplifica- 
tions often  are  possible:  the  flow  may  be  steady,  the  particle  volume  fraction  may  be  negligible,  or  de- 

viations from  equilibrium  flow  may  be  so  small  that  linearization  in  terms  of  such  deviations  becomes 
permissible.  For  steady  flow  through  a duct  of  constant  cross  section  Eqs.  (31)- (34)  can  be  integrated 
to  yield  algebraic  relationships,  and  only  Eqs.  (35)  and  (36)  remain  differential  equations.  Even  then, 
analytical  solutions  are  rarely  obtainable.  The  analysis  of  several  special  flows  is  indicated  in  the 
following  sections. 


Occasionally,  it  may  be  of  interest  to  assess  the  consequences  if  not  all  particles  are  of  the  same 
size.  The  particulate  phase  then  is  divided  into  groups  each  of  which  contains  only  particles  of  one  size 
having  its  own  relaxation  times  Ty  and  t^.  Equations  corresponding  to  Eqs.  (32),  (35)  and  (36)  then  must 
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be  written  for  each  group,  and  the  particle  terms  In  Eqs.  (33)  and  (34)  must  be  replaced  by  sums  of  the 
corresponding  terms  for  each  particle  size.  The  resultant  extended  system  uf  equations  can  be  solved  In 
the  same  manner  as  the  original  system  at  an  Increased  computational  effort,  Kllegel  (80),  Kriebel  (81), 

Soo  (2). 

Since  particles  of  different  sizes  approach  equilibrium  with  the  gas  as  different  rates,  collisions 
between  particles  may  become  important.  A technique  to  take  such  collisions  Into  account  was  suggested 
by  Marble  (28,29).  He  assumed  e “ 0 and  f « 1 and  derived  the  force  acting  between  two  clouds  of  particles 
having  diameters  Dj  and  O2  as 


(u*  - u*  )/l 


(37) 


In  this  equation,  the  asterisks  indicate  the  velocities  which  the  particles  would  have  in  the  absence  of 
collisions,  and  l represents  a characteristic  length 
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Ulus,  the  particle  acceleration,  given  by  Eq.  (35), must  be  written  for  each  particle  size  and  contains 
the  additional  collision  term  F/op|  for  the  particles  of  size  Df  and  -F/ap2  for  the  particles  of  size  D2. 
TVo  more  equations  of  the  form  of  Eq.  (35),  but  without  the  collision  terms,  must  be  written  for  u*^  and 
u*2  to  account  for  these  new  variables.  An  example  of  this  approach  is  shown  in  section  6.1. 

6.  PRESSURE  WAVES 


6.1  Shock  Waves 


An  analysis  of  shock  waves  in  gas-particle  mixtures  was  first  performed  by  Carrier  (91),  and  a number 
of  analyses  were  published  subsequently,  e.g.,  Soo  (2,79),  Kliegel  (80),  Marble,  (29,78),  Rudlnger  (53,92), 
Kriebel  (81)  (see  also  B.2.6).  If  the  shock  propagates  through  a duct  of  constant  cross  section, the  flow 
is  steady  in  a coordinate  system  fixed  in  the  shock  wave,  and  Eqs.  (31)- (34)  can  be  Integrated  to  yield 
algebraic  relationships.  Ihe  gas-particle  mixture,  generally  assumed  to  be  in  equilibrium,  enters  the 
wave  with  a velocity  that  must  be  greater  than  the  equilibrium  speed  of  sound  aeo  given  by  Eq.  (30),  where 
subscript  zero  Indicates  the  conditions  ahead  of  the  shock.  Hie  equilibrium  shock  Mach  number  is  thus 
defined  as  Me  “ u0/aeo.  the  shock  front  is  traversed  by  the  particles  in  such  a short  time  that  their 
velocity  and  temperature  has  not  changed.  Accordingly,  the  change  of  the  gas  across  the  shock  front  is 
determined  by  the  frozen  Mach  number  Mf  » u0/a0  and  is  given  by  the  well-known  Rankine-Hugonlot  shock 
relationships  for  a pure  gas.  Immediately  behind  the  shock  front,  both  the  velocity  and  temperature  of 
the  partlc  les  thus  are  different  from  the  corresponding  gas  values. 


A great  distance  behind  the  shock,  equilibrium  between  the  phases  is  reestablished  (ue  3 upe  and 
Te  " Tpe).  Hie  basic  equations  then  simplify  sufficiently  that  analytical  relationships  can  be  obtained 
between  the  upstream  and  the  equilibrium  conditions  even  for  a finite  particle  volume,  Rudlnger  (82).  If 
c ■ 0,  the  equations  become  the  same  as  for  a pure  gas  for  a Mach  number  Me  and  a specific-heat  ratio  Ym 
given  by  Eq.  (28).  Relaxation  from  the  frozen  to  the  equilibrium  conditions  must  be  obtained  by  numerical 
Integration  of  the  equations  for  assumed  correlations  for  Cp  and  Nu.  A typical  relaxation  zone  is  shown  in 
Fig.  5 for  the  conditions  given  in  the  caption.  Hie  final  equilibrium  conditions  are  indicated,  and  the 
large  extent  of  the  relaxation  zone  should  be  noted. 


It  is  instructive  to  explore  the  consequences  of  making  different  assumptions  for  Cp  and  Nu.  Results 
for  the  gas  pressure  and  the  particle  velocity  are  shown  in  Fig.  6 for  all  combinations  of  three  drag 
coefficients  and  two  heat  transfer  coefficients.  The  drag  coefficient  has  a marked  effect,  while  the  in- 
fluence of  the  hea t-transfer  coefficient  is  quite  small.  Corresponding  results  for  the  gas  and  particle 
temperatures  are  shown  in  Fig.  7 which  illustrates  that  variations  in  the  relaxation  zone  need  not  be  mono- 
tonic  and  that  heat  transfer  is  more  Important  than  drag,  as  one  should  expect  (see  also  A. 5. 3). 

If  the  particle  volume  cannot  be  neglected,  the  equilibrium  velocity  is  increased.  This  effect  becomes 
significant  for  t0  ” 0.01  and  is  quite  large  for  e0  greater  than  0.05.  In  contrast,  the  equilibrium  pres- 
sure is  decreased  but  the  change  is  almost  insignificant^  Rudlnger,  (53,82). 

Oblique  shocks  may  be  analyzed  by  considering  the  normal  and  tangential  components  as  in  the  case  of 
a pure  gas,  Morgenthaler  (93),  Okauchl  (94).  Note  that  the  equations  for  the  two  components  are  not  un- 
coupled unless  the  factor  f in  Eq.  (35)  is  assumed  constant,  as  pointed  out  in  Section  3.1. 

If  different  particle  sizes  are  considered,  Kliegel  (80),  Kriebel  (81),  Soo  (2),  the  velocity  relaxa- 
tion for  each  particle  size  appears  similar  to  the  one  in  Fig.  5 with  the  smaller  particles  approaching 
equilibrium  with  the  gas  faster  than  the  larger  ones.  Zung,  see  Marble  (29),  performed  calculations  for 
two  particle  sizes  with  allowance  for  collisions  between  the  particles.  This  example  is  shown  in  Fig.  8 
and  demonstrates  that  such  collisions  tend  to  draw  the  two  particle  velocities  together  (see  also  A. 5, 4). 

Shock  transition  requires  Mg  > 1,  but  the  existence  of  a discontinuous  shock  front  implies  Mf  > 1. 
Since  the  equilibrium  speed  of  sound  is  smaller  than  the  frozen  speed  of  sound,  a velocity  range  ae  < u0  < 
ef  • aD  exists  for  which  shock  transition  is  possible  but  a discontinuous  shock  front  is  not.  The  limiting 
value  of  Me  for  which  Mf  « 1 is  shown  in  Fig.  9 as  a function  of  the  particle  mass  fraction  for  three 
values  of  6 and  for  y " 1.4.  Within  this  range,  only  a gradual  transition  from  the  initial  to  the  final 
equilibrium  conditions  is  possible  in  the  form  of  a dispersed  shock  wave.  A special  starting  procedure  for 
the  numerical  integration  is  needed  to  obtain  initial  deviations  from  equilibrium.  Kriebel  (81)  assumed 
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initial  perturbations  that  would  be  produced  by  a weak  compression  wave,  and  Rudinger  (95)  linearized  the 
equations  to  obtain  starting  conditions  in  terms  of  a small  perturbation  of  the  gas  velocity.  Results  for 
dispersed  shocks  were  presented  also  by  Kliegel  (80).  An  example  of  a dispersed  shock  is  shown  in  Fig.  10 
and  comparison  with  Fig.  5 indicates  a much  longer  relaxation  zone. 

It  is  evident  from  Fig.  9 that  dispersed  shocks  can  have  appreciable  strength.  Dispersed  oblique 
shocks  should  be  expected  to  occur  also  in  supersonic  flow  around  sufficiently  slender  bodies.  Because 
of  the  absence  of  a discontinuous  shock  front,  dispersed  shocks  may  not  be  readily  visible  in  schlleren 
photographs,  and  no  observations  of  such  waves,  apparently  have  been  reported  so  far. 

6 ■ 2 Sound  Waves 


The  rate  at  which  particles  can  respond  to  fluctuations  in  a flow  is  indicated  by  the  relaxation  times 
Ty  and  Tj  of  Eqs.  (9)  and  (17).  Interaction  of  sound  waves  with  a cloud  of  particles  therefore  becomes 
most  Important  if  their  period  is  not  too  different  from  the  relaxation  times. 

Flow  changes  in  a sound  field  are  small,  and  linearization  of  the  basic  flow  equations  then  becomes 
possible.  Temkin  and  Dobbins  (8)  used  this  approach  to  derive  periodic  solutions  of  the  form  exp  [i(Kx-uut)] 
where  K « k^  + ik2  is  a complex  propagation  parameter.  The  real  part  k^  = ma  represents  the  dispersion  and 
the  imaginary  part  k2  the  attenuation  of  the  wave.  Some  of  their  results  are  shown  in  Figs.  11  and  12  where 
these  parameters,  in  dimensionless  form,  are  plotted  as  functions  of  idt  for  several  values  of  the  loading 
ratio.  Maximum  attentuation  occurs  near  bt  » 1 and  is  shifted  toward  higher  frequencies  with  Increasing 
loading  ratio.  The  sound  velocity,  which  has  a value  ao  for  1)  = 0,  approaches  the  equilibrium  value  ae, 
given  by  Eq.  (30),  for  low  frequencies,  and  the  frozen  value  af  = aD  at  high  frequencies;  it  varies  most 
rapidly  near  ujt  * 1.  Marble  (29)  performed  a similar  analysis.  He  compared  this  theory  with  an  earlier, 
more  detailed,  analysis  by  Epstein  and  Carhart  (96),  which  considered  the  flow  field  around  individual 
particles,  and  concluded  that  it  is  completely  adequate  for  wavelengths  which  are  long  compared  with  the 
size  of  the  particles  (see  also  B.2,5). 

6.3  Large-Amplitude  Waves 

The  general  flow  equations,  Eqs.  (31)- (36)  are  hyperbolic  and  therefore  can  be  solved  by  the  method 
of  characteristics  for  pressure  waves  of  arbitrary  shape  and  magnitude.  Rudinger  and  Chang  (87)  treated 
the  unknown  gas  and  particle  conditions  as  one  set  of  unknowns,  while  Migdal  and  Agosta  (88)  considered 
the  particle  terms  as  sources  or  sinks  of  drag  and  heat  in  a gas  flow.  There  is  little  practical  differ- 
ence between  these  points  of  view.  Both  studies  assumed  a negligible  particle  volume. 

The  six  simultaneous  partial  differential  equations  have  six  characteristics  in  the  position-time 
plane  which  everywhere  have  the  local  velocities  u + a,  u and  up,  where  Up  is  counted  threefold.  This 
degeneracy  of  the  u.-characteristlcs  was  explained  by  Sauerwein  and  Fendell  (97)  as  a consequence  of  ne- 
glecting the  partial  pressure  of  the  particles.  Flow  changes  along  the  characteristics  are  given  by 
equations  which  contain  all  the  terms  of  a pure-gas  flow,  e.g.,  Rudinger  (98),  with  additional  terms  that 
are  proportional  to  the  local  particle  concentration.  The  most  distinguishing  feature  of  large-amplitude 
waves  is  their  continually  changing  wave  form  with  the  possibility  that  compression  waves  coalesce  to  form 
shock  waves  at  initially  unknown  times  and  locations.  Numerical  solutions  therefore  may  involve  substan- 
tial computational  efforts, 

A different  approach  was  taken  by  Johnson,  et  al.  (99)  who  explored  the  "particle-ifi-cell"  method, 
e.g.,  Amsden  (100).  Their  example  of  shock  propagation  and  reflection  yielded  good  agreement  with  results 
obtained  by  the  method  of  characteristics.  They  concluded  that  the  method  could  be  applied  with  confidence 
to  problems  in  two  and  three  dimensions,  as  had  been  done  for  problems  of  pure-gas  flow. 


7.  NOZZLE  FLOW 

For  a steady  flow  through  a nozzle,  Eqs.  (31),  (32),  and  (34)  can  be  integrated  to  yield  algebraic  re- 
lationships which  must  be  solved  together  with  the  remaining  Eqs.  (33),  (35),  and  (36).  This  system  of 
equations  can  readily  be  written  in  a form  that  is  suitable  for  numerical  integration  for  a prescribed  flow 
at  the  nozzle  inlet  and  a selected  nozzle  shape  A(x),  e.g.,  Bailey,  et  al.  (84),  Hultberg  and  Soo  (101), 

Soo  (2).  In  general,  the  prescribed  inlet  flow  may  not  be  able  to  pass  through  the  nozzle  throat  (choking) 
or  may  recompress.  Instead  of  expand,  after  passing  through  the  throat.  This  behavior  is  a general  property 
of  nonequilibrium  nozzle  flows,  since  the  flow  rate  compatible  with  a throat  area  depends  on  the  relaxation 
processes  in  the  nozzle.  It  is  then  necessary  to  use  some  iteration  scheme  by  changing  either  the  throat 
or  the  inlet  flow  until  compatible  conditions  are  found. 

Because  of  the  lateral  flow  accelerations  in  the  nozzle,  the  particles  tend  to  move  toward  the  nozzle 
wall  in  the  convergent  section  and  stay  near  the  axis  in  the  divergent  section  (see  also  A. 2.1).  A schlleren 
photograph  of  a two-dimensional  nozzle  flow  with  14.4-^m  particles  by  Gilbert  et  al.  (102)  is  shown  in  Fig. 
13.  These  authors  state  that  the  result  with  2.7-pm  particles  is  "remarkably  similar."  It  can  be  seen  that 
the  particle  distribution  is  only  moderately  nonuniform  in  the  convergent  section  near  the  throat  but 
severely  so  in  the  divergent  section.  The  size  of  the  nozzle  in  Fig.  13  was  not  stated,  but  one  should 
expect  the  nonuniformity  to  become  less  severe  as  the  size  of  the  nozzle  is  Increased.  Calculations  based 
on  one-dimensional  llow  through  the  divergent  section  of  a nozzle,  therefore,  may  be  only  qualitatively 
correct. 

To  overcome  the  problem  of  uneven  particle  distribution  in  the  divergent  section,  Bailey,  et  al.  (84) 
determined  the  two-dimensional  gas  flow  without  particles  by  the  method  of  characteristics  for  steady 
supersonic  flow  and  then  computed  particle  trajectories  in  this  flow.  Since  the  effect  of  the  particles 
on  the  flow  is  neglected,  this  method  Is  valid  only  for  relatively  low  particle  loading.  A computing  pro- 
cedure, based  on  the  method  of  characteristics,  in  which  gas-particle  interaction  is  taken  Into  account, 
was  developed  by  Kliegel  and  Nickerson  (103)  and  used  to  compute  the  flow  through  various  nozzles. 
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The  performance  of  a nozzle  Is  reduced  by  the  losses  resulting  from  particle  lag.  Bailey,  et  al.  (84) 
emphasized  that  once  an  appreciable  lag  has  been  established  during  the  flow  accelerations,  the  loss  cannot 
be  recovered  by  lengthening  of  the  nozzle  to  let  the  particles  approach  equilibrium  with  the  gas.  The 
losses  can  be  minimized  by  making  the  area  changes  as  gradual  as  possible,  and  within  a prescribed  overall 
length,  an  optimum  nozzle  shape  can  be  determined.  Optimum  nozzle  contours  for  gas-particle  flow,  appar- 
ently, were  first  derived  by  Marble  (104)  from  the  linearized  flow  equation.  Later,  Hoffman  and  Thompson 
(105)  developed  a numerical  variational  technique  which  was  applied  by  Elsbernd  and  Hof&nan  (106)  to  a 
study  of  the  influence  of  various  design  parameters  on  thrust  and  optimum  nozzle  contour.  Their  results 
are  summarized  In  Table  2. 


Table  2 

Influence  of  design  parameters  on  optimum  thrust  and  nozzle  contours  (+  indicates  a 
significant  Influence,  and  - denotes  relative  unimportance). 


Parameters 

Thrust 

Contour 

Particle  size 

+ 

- 

Particle  loading 

+ 

+ 

Particle  drag  coefficient 

+ 

- 

Particle  heat-transfer  coefficient 

- 

- 

Nozzle  Inlet  angle 

+ 

- 

Radius  of  nozzle  throat 

- 

+ 

Nozzle  size 

+ 

_ 

If  the  flow  of  a gas-particle  mixture  through  a prescribed  nozzle  is  needed,  one  must  integrate  the 
equations  by  one  of  the  Indicated  methods.  More  convenient  techniques  are  available  if  a general  under- 
standing Is  to  be  obtained  rather  than  numerical  answers  to  a specific  problem. 

One  method  Is  to  assume  that  deviations  from  equilibrium  are  small.  If  one  sets  t ■ 0,  u * Up  and 
T * Tp,  the  basic  equations  reduce  to  the  well-known  equations  for  one-dimensional  nozzle  flow  with  the 
thermodynamic  properties  outlined  in  Section  4.  Small  deviations  from  this  flow  then  may  be  analyzed  by 
linearizing  the  equations  In  terms  of  a suitable  small  parameter.  Rannle  (107)  and  Marble  (29,78)  derived 
a number  of  relationships  between  the  perturbations  of  the  flow  variables. 

For  larger  deviations  from  equilibrium,  an  Inverse  method  may  be  used,  based  on  the  assumption  that 
e ■ 0 and  that  f and  Nu  are  constants.  If  the  particle  velocity  Is  prescribed  as  a function  of  x,  Eq.  (34) 
can  be  solved  for  the  gas  velocity.  The  remaining  variables  and  the  nozzle  shape  then  can  be  found  by 
successive  substitution  into  the  other  basic  equations.  This  approach,  in  which  the  nozzle  shape  A(x)  is 
Implied  In  the  Initial  assumption  for  u(x),  is  no  more  arbitrary  than  if  A(x)  is  assumed  and  up(x)  computed. 
For  example,  if  up  increases  linearly  with  x,  then  up/u  = K is  constant,  and  analytical  relationships  are 
obtained  for  all  flow  variables  in  the  same  form  as  the  nozzle  equations  for  one-dimensional  gas  flow  with 
modified  thermodynamic  properties  which  are  functions  of  K.  For  K * 1,  these  become  equal  to  those  given 
In  Section  4.  This  method  has  been  used  by  a number  of  authors,  e.g.,  Soo  (2,79),  Bailey  et  al.  (84), 
Kllegel  (80).  Nozzles  for  which  the  velocity  ratio  Is  constant  have  a shape  similar  to  that  of  typical 
nozzles,  particularly  near  the  throat;  Kllegel  called  them  "constant-fractional-lag  nozzles."  A similar 
analysis  was  performed  by  Hassan  (85)  who  assumed  up  to  be  proportional  to  x%. 

The  foregoing  analyses  apply  for  particle  loadings  for  which  the  particle  volume  can  be  neglected.  If 
the  e-terms  are  retained  In  the  equations,  numerical  Integration  Is  required.  Rudlnger  (86)  performed  such 
calculations  for  the  convergent  section  of  constant-fractional-lag  nozzles.  The  Influence  of  the  particle 
volume  becomes  noticeable  when  < approaches  5%.  In  these  flows,  the  gas  and  particle  tenperatures  vary 
only  slightly  between  the  flow  reservoir  and  the  nozzle  throat,  consistent  witi-  the  statement  In  Section  4 
that  high  loadings  tend  to  produce  Isothermal  flow.  If  one  assumes  T ■ Tp  - conet,  the  energy  equation  and 
the  heat-transfer  equation  are  eliminated,  and  the  remaining  equations  lead  to  simple  algebraic  solutions 
In  terms  of  the  velocity  ratio  K.  A comparison  of  this  approximate  method  with  the  "exact"  solution,  based 
on  numerical  integration  of  the  equations.  Is  sham  In  Fig.  14,  where  the  ratios  of  the  throat  conditions, 
marked  by  asterisks,  are  plotted  as  functions  of  the  loading  ratio  and  the  particle  volume  fraction  at  the 
nozzle  Inlet.  Results  of  the  approximate  analysis  with  t ■ 0 are  Indicated  by  the  broken  lines.  Subscripts 
r and  o Indicate  condition  at  the  reservoir  and  at  the  nozzle  inlet.  It  Is  Interesting  to  note  that  a 
range  of  T)  exists,  roughly  between  5 and  50,  for  which  the  Isothermal  flow  with  t » 0 agrees  with  the 
"exact"  solution  within  about  3%.  The  agreement  becomes  poorer  at  lower  loading  ratios,  where  the  flow 
should  not  be  treated  as  Isothermal,  and  at  higher  loading  ratios,  where  the  particle  volume  should  not  be 
neglected.  The  approximate  theory  with  * » 0 leads  to  the  particularly  simple  relationships 
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The  simplifications  that  result  from  the  assumption  of  isothermal  flow  may  have  other  useful  applica- 
tions. Rudinger  (108)  used  it  to  analyze  the  addition  of  heated  particles  to  a gas  flow. 


8.  GAS -PARTICLE  JETS 

The  limited  information  available  on  gas-particle  jets  is  indicated  in  Section  2.  If  such  jets  behaved 
like  equivalent  gas  jets,  the  large  literature  on  gas  jets  would  be  directly  applicable.  To  be  equivalent, 
a gas  Jet  should  have  not  only  the  same  average  density  but  also  the  same  mass  flow  and  momentum  flux. 

Particles  injected  into  a gas  flow  must  travel  a considerable  distance  to  reach  the  gas  velocity,  and 
they  may  impact  on  the  walls  of  the  Jet  tube  and  thereby  loose  kinetic  energy.  In  general,  one  should  not 
assume,  therefore,  that  the  gas  and  particle  velocities  are  equal.  Let  and  uj  be  the  concentration  and 
velocity  of  the  gas  Jet  which  carries  particles  at  a loading  ratio  7],  concentration  op  and  velocity  up; 
denote  by  02  and  U2  the  corresponding  quantities  of  a heavier  gas  without  particles.  Equal  average  density 
then  requires 
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where  Eq.  (23)  has  been  used.  The  flow  rates  are  equal  if  p2U2  * ajuj(l  + 1)).  After  substitution  for 
P2/al  from  Eq.  (42),  the  ratio  of  the  two  gas  velocities  becomes 
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Similarly,  equal  momentum  flux 
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and  it  follows  that 
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It  is  always  possible,  at  least  in  principle,  to  select  a gas  which  satisfies  Eq.  (42),  but  Eqs.  (43) 
and  (44)  can  be  satisfied  simultaneously  only  if  up/uj  = 1 and  then  yield  u2  = u^.  Thus,  a gas  jet  cannot 
be  rigorously  equivalent  to  a nonequilibrium  gas-particle  Jet.  Deviation  from  equivalence  can  be  expressed 
by  the  parameter 
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which  is  zero  for  up/ui  - 1 and  should  be  as  small  as  possible  if  equivalence  is  to  be  approximated.  (The 
conditions  T)  « 0 or  1\  » =°  which  also  lead  to  8 » 0 are  trivial.)  Substitution  for  11  from  Eq.  (42)  leads  to 
the  relationship  between  p and  up/uj  shown  in  Fig.  15  for  several  combinations  of  light  and  heavy  gases  and 
negligible  particle  volume  (aj  » Pi).  For  any  value  of  up/uj,  a maximum  of  p is  reached  when  p2/al  “ 1 + 
(uj/up).  The  locus  of  these  maxima  is  shown  as  the  broken  line.  Clearly,  the  particle  velocity  must  be 
greater  than  40  or  50%  of  the  gas  velocity  if  the  equivalence  condition  is  to  be  satisfied  within  a few 
percent.  Rudinger  (46)  injected  nitrogen  jets  carrying  glass  particles  of  33-p,m  average  size  into  a sub- 
sonic cross  flow  of  air  at  loading  ratios  between  2.3  and  25.  Penetration  of  the  centerline  of  the  Jet 
trajectory,  y,  was  determined  from  photographs.  Separate  experiments  with  a laser-Doppler  system  indi- 
cated that  up/ui  ranged  from  about  0.4  to  0.6  as  the  velocity  of  the  carrier  gas  was  varied  between  20 
and  90  m/s.  The  condition  for  equivalence  thus  was  approximately  satisfied,  and  the  momentum-flux  ratio, 
which  alone  determines  penetration  of  a gas  Jet,  was  computed  (with  subscript  s indicating  the  cross  flow)  as 
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Some  results  of  these  experiments  are  shown  in  Fig.  16,  where  the  filled-in  symbols  indicate  that  two 
experimental  points  were  too  close  together  to  be  shown  separately.  Ihis  figure  indicates  that  particle 
penetration  Increases  as  the  loading  ratio  decreases.  However,  schlieren  photographs  of  light  and  heavy 
pure-gas  Jets  yielded  a penetration  that  is  considerably  below  that  of  the  particles  (broken  line  in  Fig. 
16).  It  was  suspected  that  the  particles  separate  from  the  carrier  gas  and  behave  as  Isolated  particles. 
Indeed,  calculations  on  this  basis  gave  trajectories  which  could  have  been  brought  into  agreement  with  the 
experimental  observations  by  a small  shift  of  the  injection  point.  Such  a shift  should  be  expected  as  the 
result  of  the  initial  gas-particle  interaction  before  the  particles  separate  from  the  carrier  gas.  Con- 
firmation of  this  behavior  was  obtained  from  simultaneous  schlieren  and  regular  photographs  of  Jets  for 
which  helium  was  used  as  carrier  gas. 

It  is  interesting  to  compare  these  results  with  those  obtained  by  Salzman  (45)  and  Edelman  et  al.  (44) 
referred  to  in  Section  2.  Salzman  used  loading  ratios  between  about  11  and  22,  and  15-pm  particles.  These 
particles  probably  reached  a larger  fraction  of  the  gas  velocity  than  the  33-pm  particles  of  the  foregoing 
experiments.  He  assumed  velocity  equilibrium  and  derived  several  correlations  for  the  penetration,  one  of 
which  is  given  by 
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but  also  noted  a displacement  of  the  particle  centerline  relative  to  the  gas  centerline.  Edelman  et  al. 
used  1 to  j-pjn  particles  at  loading  ratios  below  0.5  and  observed  that  the  gas  and  particle  centerlines 
coincide.  Thus,  If  the  particles  are  sufficiently  small,  the  concept  of  an  equivalent  Jet  may  be  applied, 
and  existing  correlations  for  gas  jets  may  be  used  provided  the  momentum- flux  ratio  Is  properly  computed. 
Large  particles  separate  from  the  carrier  gas;  their  trajectory  can  be  computed  as  for  isolated  particles, 
but  a small  shift  of  the  injection  point  may  be  needed.  The  magnitude  of  this  shift  and  Its  dependence 
on  the  various  parameters,  such  as  the  relaxation  time  tv,  loading  ratio,  and  jet  diameter,  are  not  known 
at  this  time.  Intermediate  particle  sizes  lead  to  Intermediate  trajectories. 


9.  EXAMPLES 

A few  examples  are  Included  In  the  preceding  sections.  Several  more  are  given  In  the  following  to 
demonstrate  other  applications  of  the  analysis  to  specific  problems  or  to  Illustrate  some  aspect  by  ex- 
perimental observations.  A few  references  to  related  studies  also  are  given. 

Often  a flow  field  is  too  complicated  to  be  amenable  to  a mathematical  analysis,  for  example,  the 
particle  trajectories  in  the  flow  through  turbine  cascades.  Tabakoff  et  al  (18,19,109)  injected  particles 
into  the  flow  through  a two-dimensional  model  and  obtained  high-speed  motion  pictures.  A few  trajectories 
are  shown  In  Fig.  17  which  Indicates  that  small  displacements  can  lead  to  widely  different  trajectories. 
Impact  velocities  derived  from  such  data  are  useful  for  erosion  studies,  as  discussed  In  Section  2 (see 
also  A. 4,  B.3.5,  and  C.7  and  7). 

Particle  motion  In  a corner  expansion  of  a supersonic  flow  - Prandtl-Meyer  flow  - represents  another 
example  of  a two-dimensional  flow.  Assuming  that  the  particles  do  not  affect  the  flow,  Marble  (78)  derived 
the  equation  for  the  particle  trajectory  and  solved  It  numerically.  He  showed  that  the  departure  of  the 
trajectory  from  the  gas  streamlines  is  due  more  to  lag  In  the  tangential  than  in  the  radial  direction  (see 
also  B.2.2). 

The  development  of  chemical  lasers  requires  understanding  of  the  flow  in  small  mixing  nozzles  where 
measurements  are  difficult  because  of  the  small  dimensions  and  high  flow  velocities  and  tonperatures.  To 
assess  the  suitability  of  the  laser-Doppler  technique,  calculations  were  performed  for  a flow  of  nitrogen 
through  a typical  convergent-divergent  laser  nozzle  in  which  M = 4 is  reached  in  a distance  of  about  60  nm. 
These  calculations  were  based  on  the  drag  coefficient  given  by  Eq.  (10),  Because  of  the  low  densities  in 
the  expanding  gas,  the  correction  factor  of  Eq.  (11)  was  Included  and  reached  several  hundred  for  low  re- 
servoir pressures.  A few  results  are  shown  In  Fig.  18,  Moon  (110),  for  a reservoir  pressure  of  10  atm., 
a reservoir  temperature  of  1330  K and  calcium  fluoride  particles  of  various  sizes.  The  velocity  ratio 
up/u  is  plotted  as  a function  of  x/L,  where x Is  measured  from  the  throat  and  L is  the  total  length  of  the 
nozzle.  C' 

Continuous 
loading  ra 

ence  of  this  loading  on  the  flow,  an  equilibrium  flow  was  assumed  and  the  velocity  at  the  nozzle  exit  with 
and  without  particles  was  calculated.  The  difference  was  found  to  be  less  than  0.57.  as  long  as  the  load- 
ing ratio  did  not  exceed  0.01.  Thus,  0.1-pm  particles  would  not  affect  the  flow  significantly  (see  also 
A. 3,  C.7  and  D.5). 

To  demonstrate  how  the  nonlinearity  of  the  equation  of  motion  affects  two-dimensional  trajectories  of 
single  particles  (Section  3.1),  consider  a particle  injected  with  the  velocity  components  upQ  and  vpo  in 
the  direction  and  at  right  angles  to  the  direction  of  a constant  gas  flow  having  a velocity  u.  Using  the 
drag  law  of  Eq.  (10),  Rudlnger  (111)  derived  an  analytical  solution  for  the  maximum  particle  penetration, 

Y.  The  ratio  of  this  penetration  to  the  penetration  without  cross  flow  is  shown  In  Fig.  19  as  a function 
of  a Reynolds  number,  defined  here  as  Re  - pDvPo/p,  for  several  values  of  (u  - Up0)/vpo.  Clearly,  pene- 
tration Is  significantly  reduced  unless  u - upD  is  smaller  than  vpo.  Numerical  solutions  for  lateral  In- 
jection of  single  particles  were  obtained  also  by  Brandt  and  Perlni  (112),  who  used  various  drag  coeffi- 
cients and  Mach  number  effects  in  their  calculations.  Otterman  (113)  analyzed  the  motion  of  particles  in- 
jected into  the  boundary  layer  flow  around  the  stagnation  point  of  a body;  he  treated  the  flow  as  incom- 
pressible but  Included  lateral  forces  on  the  particle?  that  result  from  particle  rotation  and  boundary- 
layer  shear.  Hie  effect  of  four  assumptions  for  the  drag  coefficient  on  the  particle  trajectory  In  dif- 
ferent high-speed  flow  fields  was  explored  by  Korkan  et  al.  (66).  These  authors  found  that  the  results 
are  significantly  affected  for  particles  Injected  at  45°  into  a supersonic  flow  and  for  particle  motion 
In  a Prandtl-Meyer  expansion.  In  contrast,  particle  trajectories  behind  an  oblique  shock  are  insensitive 
to  the  assumption  for  the  drag  coefficient. 

Hie  foregoing  examples  deal  with  single  particles  or  particle  concentrations  low  enough  not  to  affect 
the  gas  flow,  Liu  (114)  analyzed  laminar  incompressible  gas-particle  flow  along  an  infinite  flat  plate 
which  is  Impulsively  accelerated  to  a constant  velocity  U.  He  assumed  negligible  particle  volume  and  Stokes 
drag  and  derived  t.he  variations  of  the  friction  coefficient  with  time.  Subsequently,  Hamed  and  Tabakoff 
(115,116)  analyze  the  same  problem  but  extended  it  by  allowing  for  a finite  particle  volume  according  to 
Eq.  (14)  and  by  including  In  the  analysis  the  lateral  lift  forces  caused  by  particle  slip  and  rotation  In 
the  shear  layer  near  the  plate.  Ihey  defined  a Reynolds  number  Re  * U2Tvp/p,  which  Is  thus  a measure  of 
the  plate  velocity,  and  Introduced  a transformed  friction  coefficient  Cj*  - CfU(tp/p)%  where  the  wall 
shear  stress  is  given  CfpU^/2.  Some  of  their  results,  for  Re  ■ 25  and  several  values  of  the  density 
ratio  pp/p  are  shown  In  Fig.  20.  Hie  friction  coefficient  reaches  a maximum  shortly  after  the  start  of 
the  motion  and  then  decreases  asymptotically  toward  Its  equilibrium  value.  Liu's  limiting  values  for 
frozen  and  equilibrium  flow  are  Included  In  the  figure.  For  frozen  flow  at  t - 0,  which  Is  equivalent  to 
the  absence  of  particles,  both  theories  give  the  same  value  CfQ*.  For  equilibrium  flow,  Liu's  theory  yields 
Cf  * « (1  + 7))’  which  agrees  well  with  the  theory  of  Hamed  and  Tabakoff,  It  also  indicates  a maximum  of  the 
friction  coefficient  near  t/%  » 1.5,  but  the  maximum  appears  to  be  slightly  lower.  The  good  agreement 


early,  particles  much  larger  than  U.l  pm  are  not  suitable  tracers  under  such  extreme  conditions, 
operation  of  the  laser-Doppler  system  requires  up  to  about  105  partlcles/mm2.  The  corresponding 
lo  Is  about  6 x 10*5  for  o.l-pm  particles  and  6 x 10*2  for  1-pm  particles.  To  check  the  Influ- 
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between  these  theories  indicates  that  the  additional  effects  Included  by  Hamed  and  Tabakoff  have  only  a 
minor  effect,  at  least,  at  the  low  velocities  considered.  Hie  lateral  forces  produce  a particle-free  layer 
of  thickness  yv  near  the  wall.  The  ratio  yv/D  is  shown  in  Fig.  21  as  a function  of  the  dimensionless  time 
t/rv  for  «0  ■ 0.001,  pp/p  “ 450,  and  several  values  of  a Reynolds  number  defined  as  before.  The  thickness 
of  the  particle-free  layer  is  only  a few  particle  diameters,  but  the  largest  indicated  Reynolds  number  of 
10C  corresponds  to  a plate  velocity  of  only  about  3 m/s  for  10-;*m  particles. 

The  problem  of  an  oscillating  Infinite  flat  plate  was  analyzed  by  Liu  (117),  and  steady  flow  over  a 
semi-infinite  flat  plate  was  treated  by  Marble  (29,78)  and  Singleton  (118).  Boundary-layer  flows  of  gas- 
particle  mixtures  were  studied  also  by  Soo  (2),  Otterman  and  Lee  (11)  and  Tabakoff  and  Hamed  (22). 

According  to  these  analyses,  the  pressure  gradient  in  a laminar  flow  increases  with  increasing  particle 
loading.  In  turbulent  flow  through  pipes  the  opposite  behavior  may  be  observed  under  some  conditions.  This 
effect  has  been  kncwn  for  some  time,  and  various  explanations  have  been  suggested,  as  Indicated  in  Section  2. 
Typical  results  by  Kane  and  Pfeffer  (50)  are  shown  in  Fig.  22.  These  data  were  obtained  in  a vertical  flow 
of  air  with  36-^m  glass  particles.  The  friction  coefficient  with  and  without  particles  is  defined  in  terms 
of  the  gas  variables  for  consistent  interpretation  of  the  data.  Thus,  the  pressure  gradient  for  the  flow  with 
solids  is  given  by  fspu2/2d  and  that  for  the  clean  gas  flow  by  f„pu2/2d,  where  d is  the  duct  diameter.  The 
ratio  of  the  friction  factors  fs/fg  is  shown  in  the  figure  as  a function  of  the  loading  ratio  for  several  Rey- 
nolds numbers,  defined  here  as  Re„  = pdu/^.  A marked  reduction  of  the  friction  factor  can  be  seen  to  occur 
for  some  conditions.  Drag  reduction  was  achieved  for  all  particle  sizes  (15-55  |*m)  in  vertical  tubes  but  only 
for  the  smaller  sizes  in  horizontal  tubes  Indicating  a detrimental  effect  of  sedimentation.  The  effect  is  only 
a weak  function  of  Reg.  As  pointed  out  in  Section  2 it  appears  that  this  phenomenon  can  be  expla  lned  by  the 
thickening  of  the  laminar  sublayer  caused  by  the  particles. 


10.  UNRESOLVED  PROBLEMS 

The  preceding  sections  indicate  the  variety  of  gas-particle  flows  that  are  of  practical  Importance 
and  methods  available  for  their  analysis.  Unresolved  problems  still  remain,  and  some  of  these  are  dis- 
cussed in  the  following. 

The  choice  of  an  appropriate  drag  coefficient  must  be  faced  for  almost  every  analysis.  This  problem 
is  discussed  at  some  length  in  Section  3.2.  The  customary  approach  has  been  either  to  ignore  this  problem 
and  use  the  standard  drag  coefficient,  or  even  Stokes  drag,  or  to  use  a drag  coefficient  derived  from  ex- 
perimental observations  in  a different  flow  system.  Since  determination  of  the  drag  coefficient  for  the 
system  of  interest  is  generally  not  practical,  a more  detailed  analysis  of  the  processes  that  affect  the 
particle  motion  would  be  desirable  but  extremely  difficult.  In  numerical  evaluations  of  specific  flows, 
at  least  the  sensitivity  of  the  results  to  various  assumptions  for  the  drag  coefficient  should  be  assessed 
(see,  for  instance.  Fig.  6).  A similar  uncertainty  exists  with  respect  to  the  heat-transfer  coefficient. 

This  parameter  affects  mainly  the  gas  and  particle  temperatures  and  has  only  a small  influence  on  the  par- 
ticle motion,  as  indicated  by  Figs.  6 and  7.  It  should  be  of  concern  for  calculations  that  Involve  mass 
transfer  between  the  gas  and  the  particles  or  heat  transfer  to  and  from  the  walls  of  the  duct  which  were 
not  considered  here. 

In  many  flows  of  practical  Importance,  the  gas  cannot  be  considered  as  perfect.  Work  is  currently  in 
progress  at  AECD  to  study  particle  trajectories  in  real-gas  flow  through  Prandtl -Meyer  expansions,  oblique 
shocks  and  nozzles,  using  the  latest  drag-coefficient  data. 

Throughout  this  entire  review,  the  particles  are  considered  to  be  spherical,  a condition  that  is  rarely 
satisfied  lr  engineering  applications.  Several  measures  to  indicate  the  size  of  Irregularly  shaped  particles 
can  be  defined,  Fuchs  (1),  Lapple  (27),  such  as  sieve  size  or  the  diameter  of  a sphere  having  the  same  vol- 
ume as  the  particle.  A frequently  used  measure  is  the  Stokes  diameter  which  is  the  diameter  of  a sphere 

having  the  same  settling  velocity  and  therefore  the  same  drag  coefficient.  This  measure  appears  to  be  the 
most  reasonable  one  to  use  for  calculations  of  non-equilibrium  flows.  Some  shapes,  such  as  disks  or  cylin- 
ders tend  to  settle  in  a preferred  orientation  and  may  also  perform  an  oscillatory  motion  during  settling, 
Marchildon,  et  al.  (120,121)  Christiansen  and  Barker  (122);  it  seems  unlikely  that  this  orientation  Is 
maintained  In  the  turbulent  environment  of  a flowing  suspension  with  a resultant  uncertainty  of  the  drag 
coefficient.  Relationships  between  drag  coefficient  and  particle  shape  and  concentration  remain  to  be 
established. 

The  lack  of  information  on  gas-particle  jets  is  emphasized  in  Sections  2 and  8.  If  such  Jets  are  in- 
jected into  a cross  flow,  sufficiently  small  particles  remain  substantially  immersed  in  the  gas  Jet,  and 
their  penetration  does  not  depend  on  the  particle  size.  Sufficiently  large  particles  should  be  expected 
to  separate  from  the  gas  Jet  and  move  along  a trajectory  that  depends  on  the  relaxation  time  Tv;  penetration 
then  Increases  approximately  with  the  square  of  the  particle  size.  The  extent  of  the  different  regimes  and 

its  dependence  on  the  flow  variables  is  largely  unkncwn  at  this  time, and  considerably  more  work  is  needed 

before  these  flows  are  fully  understood.  Injection  Into  a cross  flow  through  a longitudinal  slot  is  men- 
tioned in  Section  2 in  connection  with  water-augmented  air  cushion  vehicles.  No  information  on  such  flows 
appears  to  be  available  even  for  pure-gas  Jets. 


In  Section  9,  the  different  behavior  of  the  friction  coefficient  for  laminar  and  turbulent  flow  is 
discussed,  and  it  is  clear  that  simply  replacing  the  viscosity  by  a constant  eddy  viscosity  could  not  account 
for  the  observed  variations  of  the  friction  coefficient  in  turbulent  flow.  As  in  the  case  of  turbulent  flow 
of  a pure  gas,  the  basic  equations  should  Include  terms  corresponding  to  the  Reynolds  stresses.  Such  equa- 
tions were  derived  by  Panton  (123),  but  no  models  have  been  developed  so  far  to  evaluate  these  additional 
terms. 

There  are  many  still  unresolved  problems  also  in  the  fields  of  instrumentation,  design  of  engineering 
equipment,  and  analysis  of  flows  that  Involve  processes  specifically  excluded  in  the  preceding  sections. 
Their  discussion  would  go  beyond  the  scope  of  the  present  review. 
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Figure  1 Standard  drag  coefficient  for  sphere^ 
after  Schlichting  (58).  (Reprinted  with 
permission  of  McGraw-Hill  Book  Co.) 
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Figure  2 Sunmary  of  sphere  drag  coefficients. 
Ballistic-range  data  of  Bailey  and  Hyatt 
(65).  (Reprinted  with  permission  of  AIAA) 
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Figure  3 Drag  coefficient.*  obtained  by  several 
investigators  for  various  flows,  after 
Soo  (2).  (Reprinted  with  permission  of 
Ginn  and  Co.,  Xerox  Corp.) 
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Equi librium  speed  of  sound  for  gas 
y*  icie  mixtures,  after  Rudinger  (82) 
(Reprinted  with  permission  of  Marcel 
Dekker,  Inc.) 


Figure  5 Relaxation  zone  behind  a shock  front 

for  10->m  glass  spheres  in  air  at  standard 
temperature  and  pressure,  tor  T]  » 0.2  and 
Mt  ■ 1.30,  corresponding  to  6 ■ 1.125,  \M  ■ 
1.20,  aeo/aQ  « 0.88  and  Me  * 1.70;  based  on 
Ingebo's  drag  Eq.  (13),  and  Nu  * 2,  after 
Rudinger  (53).  (Reprinted  with  permission 
of  Marcel  Dekker,  Inc.) 
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Figure  6 Effect  of  different  assumptions  for  drag  and  heat 

transfer  on  pressure  and  particle  velocity  for  the  same 
mixture  and  shock  strength  as  for  Figure  5,  Rudlnger 
(53).  (Reprinted  with  permission  of  Marcel  Dekker,  Inc.) 


Figure  10  Dispersed  shock  wave  computed  for  the  same  conditions  as  for  Fig.  5 
but  for  Mf  " 0.95  corresponding  to  Me  = 1.08,  after  Rudinger  (53). 
(Reprinted  with  permission  of  Marcel  Dekker,  Inc.) 


Figure  11  Attenuation  of  sound  waves  by  water  droplets  in  air;  6 * 4.17,  Pr  * 
0.17,  v “ 1.4,  a'ter  Temkin  and  Dobbins,  (8).  (Reprinted  with  permission 
of  the  American  Institute  of  Physics) 
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Figure  12  Dispersion  of  sound  waves  by  water  droplets  in  air  for  the  same  conditions 
as  for  Fig.  11,  after  Temkin  and  Dobbins  (8).  (Reprinted  with  permission  of 
the  American  Institute  of  Physics) 
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Figure  14  Comparison  of  the  throat  conditions 
based  on  "exact"  and  isotnerroal  flow, 
after  Rudinger  (53).  (Reprinted  with 
permission  of  Marcel  Dekker,  Inc.) 


Figure  15  Equivalence  between  gas  Jets  and 
gas-particle  Jets,  Rudinger  (46). 
(Reprinted  with  permission  of  ASME) 
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Figure  16  Correlation  for  the  pi 
gas-particle  Jets  into  a c: 
Rudinger  (46).  (Reprinted 
permission  of  ASME) 
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